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AUTHOR'S ABSTRACT 


Sbornik zadach po dinamike tochki v_ pole tsentral'nykh sil, 


[Collection of Problems on the Dynamics of a Point in a Central 
Force Field], Ye. N. Polyakhova, Leningrad, Leningrad University 
Press, 1974, pp. 1-145. 


The collection is a detailed selection of problems on the 
dynamics of the motion of a material point acted on by a central 
gravitational force of attraction, in particular, . the dynamics 
of space flight. As an exception, the book presents several prob- 
lems on the motion of a point acted on by central nongravitational 
forces. 


The collection is intended mainly for correspondence students, 
however it can also be used as a text in the course on theoretical 
mechanics for students in day and evening departments. It may 
also prove useful to instructors providing practical exercises in 
the course on theoretical mechanics. 37 illustrations, 6 tables, 

9 references. 
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FOREWORD 


This collection is a textbook for the course of theoretical 
mechanics ("Point Dynamics" section). It is intended mainly for 
students in correspondence departments of Leningrad State Univer- 
Sity and other higher educational institutions. It may also be 
used in part by students in the day and evening departments. 
Moreover, the problem book may prove useful for beginning instruc- 
tors in providing practical exercises in the course on theoretical 
mechanics, particularly, when they prepare modifications of 


test problems. 


Most of the collection is a detailed sampling of problems on 
the dynamics of a material point acted on by gravitational force, 
in particular, problems on the elementary dynamics of space flight. 
Several problems on the motion of a point acted on by central non- 
gravitational forces are presented. Altogether, the collection 
includes about 200 problems of varying degrees of difficulties, 
with solutions. 


The basis for Sbornik zadach was material from lectures and 
practical exercises in the course on theoretical mechanics given 
by the author for a number of years in the mathematics and mechanics 
division Jof Leningrad State University. The conditions of the 
problems were set up by the author or else were taken from various 
domestic and foreign texts, indicated in the bibliography. Solutions 
to all problems, as well as related computations, have been provided 
anew by the author or have been carefully verified. 


It should be noted that the collection presents solutions of 
most of the problems proposed in the new chapter "Dynamics of 
Space Flight" from Sbornik zadach po teoreticheskoy mekhanike 
[Collection of Problems on Theoretical Mechanics] by I. V. Meshcher- 
skiy. 


The proposed problems, as related to their subject matter, 
are grouped into 11 sections. At the beginning of each section 
(not including the last) a brief theoretical background and 
essential formulas on the dynamics of a point are presented. The 
section "Miscellaneous Problems" groups problems in whose solution 
information from various parts of the collection is required. 


A distinguishing feature of this Sbornik zadach is the expo- 
sition of a number of theoretical questions on the dynamics of a 
point in a central force field in the form'of problems, so that 
quite often references to the numbers of formulas are derived 
directly in the solution are encountered in the text. The author 


LV 


hopes that the presentation of several elements of theory in the 
form of problems will promote an easier grasp of the material 
when presented in the correspondence teaching form. 


The author is obliged to thank the head of the celestial 
mechanics faculty of the mathematics and mechanics department of 
Leningrad State University, Professor K. V. Kholshevnikov and the 
docent of the theoretical mechanics department of this same faculty, 
oS. A. Zegzhda, for attentively reading the manuscript, for valuable 
counsel and comments that did much to promote improvement in the 
book. The author’ will aiso be appreciative of all who wish to 
report errors found or to express their critical comments and 
Suggestions. 
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CHAPTER ONE 


CENTRAL FORCES. FORCE OF GRAVITY AND ITS 
DYNAMIC CHARACTERISTICS 


In mechanics a force F whose line of action extends through 
a point called the force center (center of attraction) is called 
&a central force. The central force can be expressed by the 
formula 


Fe Fir, | (1.1) 


where r- is the unit vector of the radial direction connecting 
the force center with the point of its application. In attrac- 
tion the magnitude of force F is negative, and in repulsion -- 
-- positive. 


Central forces are subdivided into forces of attraction 
(directed toward the center) and forees of repulsion (acting from 
the center). Among the central forces most widespread in nature 
we should note first of all mutual attraction, that is, the gravi- 
tational force of attraction. Its law of variation, formulated by 
Newton, is known as the law of universal gravity. 


The law of universal gravity states that the masses m and M 
mutually attract each other with a force directly proportional to 
the product of these masses and.inversely proportional to the 
square of the distance r between them: 


[Fl= fam. (1.2) 


re 


The coefficient of proportionality f = 6.673°107° om>/g* see" is 
called the universal gravitational constant. 


* Numbers in the margin indicate pagination in the foreign text. 


From Newton's second law it follows that masses and accelera- 
tions ina gravity field ‘are associated by the relation 


Wt. We | (7.3) 


According to Eq. (1.3) two problems in the dynamics of a point 
in a central gravity force field ;are distinguished: 


1) the motion of a point with small mass m in the gravity | 
field (of a point with large mass M (m << M), where the ratio 
of masses is such that acceleration Wy can be neglected and there- 


fore we can assume mass M to be fixed (the limited problem of two 
bodies, or the problem of the motion of a "nonattracting" point 
in the gravity field of ja "attracting" point); and 


2) the motion of a point with mass m in the gravity field 
‘of a. ‘point with mass M (m < M), when acceleration Wy cannot 


be neglected (the problem of two bodies, or the problem of the 

motion of one "attracting" point jin the gravity field of | 

another). /6 
Problems: of the first type include all problems on the motion 

of artificial satellites and spacecraft in. central gravity fields 

of planets for the Sun, that is, on the dynamics of space flight. 

Their solution then makes up the bulk of the collection. Problems 

of the second type are examined only as exceptions. 


If not specifically stipulated, we will not differentiate 
from the dynamic point of view the concept of the gravitationally 
attracted material point and the gravitationally attracting body: 
in the class of these problems we will assume bodies whose dimen- 
sions cannot be neglected to be homogeneous bodies, and the attrac- 
tion of a homogeneous sphere, as we know, is equal to the attraction 
of a material point coincident with the center of the sphere at 
which the entire mass of this sphere is concentrated. 


Problem 1.1. Set up an equation of the limited two-body 
problem, that is, equations of motion of a point with mass m in 
a field of central attractive force with fixed point mass M (m << M). 


Solution. Let us place at the point of mass M the origin of 
an absolute inertial coordinate system x, y, z- The position of 
the point with mass m relative to this origin we will characterize 
as the radius-vector r= rr’. The central attractive force acts 
on mass m. According to Eqs. (1.1) and (1.2) this force is of the 
form: 


Farr, Fe-fmft,. (1.4) 


pe 


‘and the equations of motion of a point acted on by this force are 


2 =. ar 2= _ 
mi E=-fmM pe) op mo b=-fme z. 


By dividing both parts of the equality by m and denoting fM = HL; 
we get finally 


d°r + ps =Q.: (1.5) 


Eq. (1.5) is called the vector /equation of motion for a limited 
two-body problem. From this equation it follows that the motion 
(acceleration of a point depends only on the gravitationally — 
attracted point mass M concentrated at the origin of coordinates. 

The coefficient u = fM, called the gravitational parameter of the /7 
eentral body, characterizes the intensity of the gravitational 

field induced by mass M. 


The vector Jequation (1.5) is equivalent to three equations 
expressed in coordinate form: 


- 


ax , x = 0 d?y dz é 7 
diz * Ph ps => oH + #520, deat Mh ye7 0. (1.5"') 
Problem 1.2. Derive the formula of action of a central 


attractive force of fixed mass M as a point with mass m moves in 
the field of this force. 


Solution. To derive the formula of the elementary work done 
by the attractive force, let us use the expression for force (1.4): 


1 This mass is also called the central mass, or the mass of the 
central body (if the body is a homogeneous sphere), since the line 
of action of the force continuously passes through this point mass, 
also called the force center. The term "central mass" by no means 
signifies that this mass is at the center of the trajectory of 
motion, although in a particular case this is possible. 


SA-Fedpa— fmM (F.gp)-- fa ggz)- FmMg (1.6) 


By integrating (1.6), we get expressions for the work done by the 
attractive force as a point mass m at the surface of the imaginary 
sphere with radius ry is moved to a new position on the surface of 


a sphere with radius Py! 


: ramet). | (1.6") 


A- it fut dr = ffl 


If a point with mass m approaches mass M(r., < ri); the work 
done by the attractive force is positive. If point m separates 
from this mass (Pr, > ri); the work is negative. 


Problem 1.3. Show that the attractive force is a potential 
force and determine the form of the potential function. 


solution. The function of coordinates U whose differential 
is equal to the elementary work done is called a potential function, 
or a force function, while the force and force field for which this 
function exists are potential force and force fields. 

Based on Eq. (1.6) we conclude that dU = 6A = F-dr = ~ —= dr. 


Thus, the potential function exists and is of the form 


We can determine the constant-C from the relation 
C=lim U=U,, : (1.7) 


that is, C is the value of the potential function at infinity, 


and it can be set equal to zero: C= Ue = 0. By thus fixing /8 


the arbitrary constant, we get the expression: 


,; M _ fmm 
Us fim + U, = r ? 


(1.7') 


which we will call the potential of the attractive force (the 
potential function with fixed arbitrary constant is usually called 
the potential). 


The set of services of the potential level for different r 
values is a set of concentric spheres whose common center coincides 
With the center of attraction where mass M is concentrated. 
Equalities (1.7) and (1.7') enable us to write 


fra sU-Us= Aeon s | (1.7") 


so that the gravitational potential at a given surface of a level, 
that is, at the surface of a sphere with given radius is equal to 
the magnitude of the work done by the attractive force (potential 
difference) as point mass m is moved from infinity to this surface 
of the level. The attractive force can be expressed here in terms 
of the potential: 


F = U/Ar =- fmMIr?, 


Problem 1.4. Show that the gravitational potential of a 
homogeneous sphere with mass M is equal to the potential of a 
point mass m equal to the mass of the sphere concentrated at its 
center. 


Solution. Knowing the expression for the gravitation poten- 
tial of point mass (1.7'), we can compute the potential for point 
mass m attracted by a homogeneous sphere with mass M and radius R. 
Let .us subdivide the elementary volume of the sphere dt and calculate 
it in a system of spherical coordinates p, 98, and (Fig. 1) bound 
with a fixed direction to point mass m at distance r from the center 
of the sphere. We can easily see that the elementary volume dt = 
= do°pdé°p sinédd, therefore the potential dU of this volume for 
point m at distance h:-from element dt will be a= fmdtt, . fmydr | 5 


where hoa r?4p?-2rp cos 9; v is the density of the sphere. By 


integrating in p from 0 to R, in 6 from 0 to mw, and in ¢ from 0 to 
en, for the entire sphere we get 


Pa \\ ind Bde _ 19 
Yrs p?- 2 rp cos B 


000 


-oatms |(dpseb ds | 
J ) r + p*- arp cosé ! 


We set W=U(O) = 724 p2)-2rp cos§ , | so that au -2 rp sin TA 


Therefore 
___psnod® = (du _[1 (3, 92 2rpcese)*] ~ 
aera r \Se af (v's pe 2epeest) - 
< [r+ p) -[r-pl]= 22 | 
r r 
SO 


R 
U= afm \ tap =(4 0%) fm. ! 
0 


But the mass of the homogeneous 
sphere M = 4g R?y, ; « SO the 


gravitational potential of a homo- 
geneous sphere with mass M for 

| point m at distance r from the 
center is 


| u-taM, | (1.8) 


that is, it is equal to the poten- 
tial of point mass M concentrated 
at the center of the sphere. Here 
the force of gravitational attrac- 
tion of the sphere is also equal 
to the central force of attraction 


of the point mass: 


We note that the result is valid not only for a homogeneous sphere, /10 
but also for a sphere with spherical distribution of densities, 

that is, when the densities are equal at points equidistant from 

the center. In this case y.4n -9(p) p*d p. 


Based on the foregoing we can conclude that two spheres 
(homogeneous or with spherical distribution of densities) are 
attracted as corresponding material points. This important property 
of the potential of the sphere enables us within the framework of 
the problems considered in this collection to assume the gravitation 
of planets and the Sun to be gravitation of point masses, considering 
that the distribution of the densities of these bodies can be assumed 
spherical, with adequate accuracy. 


Problem 1.5. An artificial Earth satellite with mass m is 
moving uniformly in a circular radius at given altitude H acted 
on by the Earth's gravitational force. Write the law of action 
of the force as a function of angular and linear velocities of 
motion. Determine how many revolutions per day are made by the 
"Zero", that is, fictive, satellite around the Earth, moving along 
the surface of its spherical Earth in the circular orbit. 


Solution. In uniform motion along the circle, the accelera- 
tion w of a point consists only of the radial component w= wl = 


r 
= -w°(R, + H}. Here Re is the radius of the karth, and wis the 
O - 


angular velocity of the AES. Then in accordance with the equations 
ofmotion we can write 


F=F. =-mw"(R, +H), (1.9) 
} 
or 
af «- mv ; | ' 
- Ry +H (1.9') 


where the angular velocity w and the linear velocity v are associated 


by the relation ye(R,+H)w.' Eqs. (1.9) and (1.9') characterize 


the Gesired law of action of the force. 


We know that the "zero" satellite (r = R,) will move with 


linear velocity v = 7.9 km/sec (first escape velocity). Let us 
find its angular velocity: 


- 0. _70kmSec, 86400 _ 79-8,64-10* _ 
O= "R6370 km 628 6, 3710°-6,a8 ~ 17.1 revolutions. 


Thus, actual Earth satellites cannot make more than 17.1 revolutions 
per day. 


Problem 1:6. An artificial. Earth satellite with mass m moves 
uniformly in a circular orbit at given altitude H acted on by the /11 
Earth's gravitational force. Knowing that the mutual attraction 
of Earth and AES [artificial Earth satellite] obeys Newton's law, 
express the force of attraction in terms of the acceleration due 
to gravity at the Earth's surface and at altitude H. 


‘Solution. For this problem Newton's law is of the form 


\F | = fmm, / (Ry + H)°. We know that at the Earth's surface _{H = 0) 


the attractive force is numerically equal to the weight of the 
body: |F| = D = mg (we will assume the Earth to be fixed, and 
thus we will neglect the inertial forces). Hence fmm _, 
or RE MGs | 


FM, = yR?. (1.10) 


The attractive force F can be expressed here in terms of the 
acceleration due to gravity g: 


mg Re . R? 
[FI = “RE? F<- ~-2oe Fo (1.11) 


Note also that from Eq. (1.10) there derives a formula for computing 
the acceleration due to gravity at the Earth's surface: 


RE (1.11') 


_ £mM _— mgRe | 
|= qe = TRAE = Oo 


is the acceleration due to gravity at altitude H. Finally, we get 


For arbitrary. altitude H, where By 


(1.11") 


. ; 2 | 
MMe GRy=9 (RytH). 9 = yr Ate | 


Problem 1.7. To which altitude must a circular-orbit Earth 
satellite moving in the plane of the equator be inserted into 
orbit in order to be continuously over the same point on the 
Earth's surface (Ry = 6370 km). 


Solution. To fulfill this condition for stationary status 
of the "Satellite track" point it is necessary that the angular 
velocity of the AES in orbit be equal to the angular velocity of 
the Barth's rotation. “8 which can be readily calculated by the 


2K -5 -1 
formula ®, SR aR 7°10 ~ sec (P, is the period of the 


Farth's rotation about its axis). This AES is called a diurnal 
(24-hour) stationary [geostationary] satellite. 

To solve the problem, let us use Eq. (1.9), which for our 
case Must be rewritten as {F{=mw} (R,+4),' and Eq. (1.11), 


enabling us to express the force in terms of the acceleration 


due to gravity. By equating the force expresstons, we find the 
3 
formula for the altitude H = yh - R. By substituting in it 


= 9.81 m/sec“, ae 6370 km, we get H % 35,800 km. 


Problem 1.8. Determine the velocity required for rectilinear 
vertical climb of a missile with mass m from the Earth's surface 
to altitude H in the Earth's gravitational field. 


Solution. In the ascent of a missile from the Earth's sur- 
face to an assigned altitude H there is a change in the kinetic 


energy of the missile equal to the work done by the Earth's 
gravitational force (see problem 1.2). In this case we have 


2 
mph — mph = Frit, (gor te) = eae 
At altitude H the final velocity is equal to zero (the missile 


comes to a halt), and the initial velocity v= Vata. 


or after replacing oM. by gR¢ according to Eq. (1.10), we have 


O 
v= OgR, .! 
a” R,+F (1.12). 


Let us examine the two limiting cases: 


1) for H << Ry (low altitudes) vj-P2gh|io k= Vig | (Gali- 


leo's formula), and 

2) for H + » (separation from the Earth's gravitational 
field), lim Y. = Tim, V 2gR, [H+ Rs - /29Ry - 11.19 km/sec (second escape 
velocity). | 

Problem 1.9. Determine the velocity of fall of a point with 
mass m to the Earth's surface if this point is dropped at altitude 


H without initial velocity. 


Solution. The problem is the inverse problem 1.8. In this 


— Ge, 
case there is a change in the kinetic energy of the missile: = —) 
frm MeH _ . 
mk = fms ($- +)? GH’ Then vy = 0, but the velocity 


™N 


10 


of landing is 


~ 2 MsH . 2g 4 . f 
U = YR (CR, +H) ~Y Ry+H ( (1.12°') 


We can see that the velocity of landing found from Eq. (1.12') is 
the same as that caleulated by the formula of the launch velocity 
for a climb to altitude H (1.12). For H << R,, we again have 


5 
v % V2eH, and for H > » (arrival from "infinity") v%® Y2geR, = 11.19 


km/sec. Hence, in particular, it follows that meteorites falling 
on the Earth "from infinity" in a parabolic trajectory can have 
this velocity. 

Problem 1.10. Determine the velocity Vo that must be imparted 


along a vertical directed upwards to a body at the surface of the /1 
Farth in order for the body to ascend to an altitude equal to the 
Barth's radius (H = R, = 6370 km). Here we will assume only the 


Rarth's gravitational force. The acceleration due to gravity at 
the Earth's surface g = 981 em/sec* = 9.81°107> km/sec’. 


Solution. Substituting into Eq. (1.12) the value H = Ryo let 


us find v,=y29R,'/2R, = VgR, = ¥ Q,81-10°- 6570 ~' 7.9 km/sec. 
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KEPLER'S LAWS 


Kepler's laws on planetary motion are usually formulated as 
follows: 


1. Hach planet moves in an ellipse, at one focus of which 
is the Sun. 


2. The area of the sector described by the planet's radius- 
vector changes in proportion to time. 


3. Squares of the periods of revolution of the planets 
relate as cubes of the semi-major axes of their orbits. 


Problem 2.1. Show that the coordinates of’ a planet moving 
around the Sun according to Kepler's laws can be expressed as a 
function of time. 


Solution. Let us turn to Kepler's first law and consider an 
ellipse (Fig. 2) described by point M (planet) whose center C we 
will take as the origin of coordinates. We will use as the Cx axis 
the direction of the semi-major axis Cn = a@ extending through 
focus S at which the Sun is, and as the Cy axis -- the direction 
of the semi-minor axis CB = b. Let us replace the equation of an 
ellipse in this coordinate system 


2 2 
By + at (2.1) 


by parametric equations, selecting the parameter as follows. From 
point M (x, y) representing the position of a planet, let us drop /14 
a perpendicular MN to the X-axis. Extending this perpendicular 
upwards until it intersects a circle whose diameter is the major 

axis am, we get point M'. Specifying each of the points M and M'! 
uniquely defines another point. But the position of point M' can 

be characterized by angle E for the center of an ellipse measured 

from the semi-major axis Caw along the line CM' in the direction 


11 


of planetary motion. This angle 
is called the eccentric anomaly 
of the planet. Eq. (2.1) is 
equivalent to the parametric 
equations of an ellipse xX = a cos B£, 
and y = b sin E, the first of which 
| is obvious from geometrical consid- 
erations (x = CN), and we obtain 
| the second after substituting the 
| first into the canonical equation 
(2.1). The ratio e = CS:Cn < 1 


a“ 
\ 


\ S defining the shape of:the ellipse 
~\ rs ® e * 
~ Pa is its eccentricity. Obviously, 
——~ M CS = ae, CB = b = av 1-e%. 
Fig. 2 Not let us introduce a system 


of rectangular orbital coordinates 
S&n, whose axes are parallel to 
the system axes Cxy, and whose origin is at focus 8S. Then 


“§ aCN-CS = @ cose ~ we, p=WM =u V1~0? sink, | (2.2) 


Let us denote by r and ¢ the polar orbital coordinates corresponding 
to € and n. The angle $¢ measured from the radius vector Sq oriented 
at the perihelion 17 is called the true anomaly of the planet. Since 


ea rcosg, T=TsN Gs | (2.3) 


then, by comparing equalities (2.2) and (2.3) we get the formulas 
rsing =a yvi-e* sink, reos p= a(vosE-e), (2.4) 


which serve in calculating the polar orbital coordinates r and ¢ 
for a given E. 
From Eqs. (2.4) we can easily determine that 
r=a({-ecosE), | (2.5) 
Subtracting from (2.5) the second of the equalities (2.4)..and then 


adding them, we find r (1-cas Y) = a ({+e)(4- cos Eg) , rn (44008 Q¢) =| 
ea({-e)({+cosE): » OF | 


VF sin} =¥aGre) sin» YP cos y=7 x(i-2) (0s a ( 2. 6 ) 


where in extracting the root the sign is determined uniquely, 
Since the angles ¢/2 and E/2 are always in the same quadrant 
(E = 180° corresponds to = 180°). From (2.6) there also follows /15 


gers we. (2.7) 


Thus, to calculate r and 9» for a given E we can use, instead of 
(2.4), Eqs. (2.6) or Eqs. (2.5) and (2.7). 


Note that by canceling out E from Eqs. (2.4) and (2.5), we 
get the equation of an ellipse in the polar coordinates: 


_ Bp | 
= Trecs@ ° | (2.8) 


where p-a@(f-e*)=b¥i-e7 4s the focal parameter of ‘an ellipse, 


that is, the ordinate of a point for which ¢ = 90°. Eq. (2.8) 
is also the general equation of conic sections (when e = 1 it is 
a parabola, and when e > 1 it is a branch of a hyperbola concave 
with respeet to focus S). 


It remains to. show that angle E can be expressed as a func- 
tion of time. To do this we turn to Kepler's second law. 


If we let vo and t refer to the moments of time at which the 


planet is at the perihelion and at arbitrary point M, and if we 
let T represent the planetary period of revolution, based on the 
second law we have 


Q _ tr-éx ; 
Rab6|UT!”Ct”C (2.9). 


where nab is the area of the ellipse, Q@ is the area of the focal 
sector SM equal to the difference between the areas of the curvi- 
linear trapezium NMr, and a triangle NMS, that is, 


Q is area NMn -~ area ANMS. (2.10) 


To calculate the area of the trapezium NMn let us use the 
ratio of the areas of the two curvilinear (elliptical and circular) 
trapez1a: 


area NMn . oo 


area NM'n #4 (2.11) 


13 
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The area of trapezium NM'a in turn is half the area of the circular 
segment M'nM" and can be determined from the familiar formula of 
the circular segment with aperture angle 2E: area NM'aw = 1/2: 


“or /2 (2E-sin 2E) = 1/2 a2 (E-sin E cos E), whence by means of 
Eq. (2.11) we get area NMa = 1/2 ab (E-sin E cos E). The area 
ANMS appearing in Eq. (2.10) can be determined by means of Eqs. 


(2.3) and (2.4): area ANMS = $9 C8, )=—H SING cose =-4 ab sin E(cosE-e),' 
The "minus" sign in front of &, corresponds to the "minus" sign f16 


in Eq. (2.10), from which we finally have Q = 1/2 ab (E - e sin E), 
owing to which Eq. (2.9) can be written as 


E-e sin E=M, (2.12) 


OT 
where M=n (t - t_)> n= a- 

The quantity n, the mean rate of change of angle E, that is, 
the mean angular velocity of the planet,. is called the mean motion 
of the planet, M is the mean’anomaly, and Eq. (2.12) is Kepler's 
equation (see also Section [Chapter] Six). 


Kepler's equation enables us to completely solve the problem 
of determining the angle E as a function of time. Thus, based 
on the results obtained above it can be stated that the coordi- 
nates of a planet moving around the Sun according to Kepler's 
laws can be expressed as a function of time. 


Problem 2.2. Using the formula r = a(l = e cos E), derive 
the formulas of the aphelion ry and perihelion re distances of a 


planet from the Sun and the orbital eccentricity (Fig. 2). 


Solution. For the perihelion st and the aphelion a, angle E 
is 0 and 180°, respectively, so that 


rx a(i-e), nate), a-$ (ratte). (2.13) 


Hence the eccentricity is 


e--2 23, 
— Ty + Tx ' (2.13°) 
and the ratio of distances is 
7x . {-e , (2.13") 
Va T+e 


We can easily see that the mean (arithmetic mean) distance of a 
planet from the Sun is equal to the semi-major axis; 


_ 1 _ 
ray 7 5 (Pe + r a. 


Problem 2.3. The semi-major axis of the Earth's orbit as it 
moves around the Sun is 149.6°10° km, and the orbital eccentricity 
of the Earth e = 0.01678. Caiculate the largest re and the 


smallest re distances from Earth to Sun. 


Solution. From Eqs. (2.13) we have r= 149.6*10° (1-0.01678)* 


% 147.1°10° km, re 149.6°10° (14+0.01678) = 152.1°10° km. 
Problem 2.4. A satellite moves around a planet with radius 


R in an elliptical orbit, whose eccentricity is e. Find the semi- 
major axis of the orbit if the ratio of the pericenter altitude 


to the apocenter altitude H O/H, =y<l. 


~. 
fH 
=~] 


Solution. For an elliptical orbit we can write 2a = r_ t+ 
+ ri Hs + He + QR = HA (Ll + H/H,) + 2R = rick) 1 + y) + QR, 
but vr. = o(1 + e), therefore 20-2R = a(l + e})(1 + y-R(1i+ y), 
aLl - y-e (1 + yi) = R, whence we have « = R (1 - y)/l - y-e 
(1+ y). 


Problem 2.5. Determine the mean radius-vector in time [rJ], 


if a is the semi-major axis of the orbit and e is its eccentricity. 
Consider the values [r] for averaging with respect to the other 
variables of elliptical motion. 


Solution. By a time-averaged radius-vector of a point moving 


in an elliptical orbit we mean the quantity [r],= 4 J rae 


where T is the period of revolution. Let us use the formulas of 
elliptical motion obtained in problem 2.1, 


‘paa(1-e cose), E~e cose = n(t-T)=M,° 


and let us replace the variable of integration, ‘changing from t to E 
in the formulas 


dE -ecosEdE =ndt =dM, dt = +-200s£ gf,, 


Let us rewrite the formula of the time averaging in terms of a 
new variable, by changing the limits of integration for one 
revolution of the satellite: 


LS 
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—~2e cos€ + e7cos?F dé = be (-2esinEs | 


2x { 


- at, 2. — e e 
= £5 z(i+te } 


i 2 e* .. ) 
tz e-F— er a singE 


Considering that the period of revolution T = 21/n, we finally 
get 


Ir], = a (4+ +4 e*).' 


We can easily show that the same exact result ‘can be obtained by 
applying the formula for averaging r with respect to mean anomaly 
M. The latter is an analog of time in formulas of elliptical 
motion. The period of its variation is 217. In this case 


2x 


{ 
47 BK ya 


~ 
pa 
Co 


where 


| 


r=a(t-ecos£), AM= ({-ecosé )dé, | 


so that 


au . 
(rl, = | | (1- ecose)* dé = a(itgse*), 
5 t 


While averaging over the mean anomaly is equivalent to aver- 
aging over time, averaging over the eccentric anomaly E gives a 
different result: 


2% ot : . _ 
(rl, =o \rae = H& \(!- ees &) aE ae(& -esinE)| a 
0 | 0 


that 18, we again obtain the arithmetic mean Tay =o = mn Tas 


(see problem 2.2): averaging over E leads to a "loss" of ellip- 
ticity; ellipticity is manifested when eccentricity is present 
in the formulas derived above for [rd, = irJy 


Let us estimate roy and [rl, = [ry for the Earth's orbit 
(e = 0.01678). a =r, = 1.4960000°10!% em is the exact value 
of the astronomic unit (Earth-Sun distance). Estimate irl, 
= [r], =. (1 + § e%). = 1,0001408 = 1.4962106"10"> em. Within 


the limits of accuracy needed to solve problems in this collection, 
we can neglect the difference between these quantities and assume 


a oa. ..,0 
that ir], “0 Poy FOF L49.6°10° km. 


Averaging over the polar angle » (true anomaly) reduces to 
computing the integral on 


[r], = a | ra@. 


0 


This integral can be easily obtained by means of the subtitution 


tag -/ ree Wf 


and when this is differentiated we have 


2 i+e ,_o@F _, 
(1+ y°$) de - {-€ cos?— 


Replacing tg £ by its value from the substitution, we get 
tte . sin? = ! 
f+ tg $= 4+ 72S tye = 14449 ws ok? 


which enables us to express add in terms of dE: 


ag = tte , @eE REE 
{-e + {+2 sin’ 


Using the formulas 


cos = t-(1+cos&), sin? = 4 (4 ~cos E), | 


1 
2 
let us find 


1-e 4-ecos— {-ecose r 


dy = ise (j-e)d& _ yi-e* dé = fe dE 
Y 


if 
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so that we finally get 


™ 
aad 
\O 
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CHAPTER THREE fi 
INTEGRAL OF AREAS 
The theorem of the variation in the kinetic moment K (angular 


momentum) of a point with mass m moving under the effect of a 
central force F is of the form 


aK _. a (Fx mi) -FxF=0.! 


Hence follows the law of conservation of kinetic moment when the 
point is moving under the influence of central forces, usually 
called the integral of areas: 


rFxv=t, (3.1) 


where the vector constant c is the constant of the areas. 


To determine the geometrical meaning of c, let us introduce 
the vector Aco, whose modulus is the area AOMM' (Fig. 3): 


bG= + (F x ar). 


Dividing both parts of the equality by at and letting at tend to 
zero, we get the differential formula 


© ._ wry @F ~~ 
Sat = TX ge = TxD, (3.2) 


where do is a vector whose modulus is equal to the elementary area 
swept by radius vector r in time dt. A comparison of 8.23 with B.D 
enables us to write 


(3.2') 
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” The quantity do/dt is called the 
sectorial velocity of a point so 
that the constant of areas is 
equal to double the sectorial 
velocity of a point. The modulus 

| of the sectorial velocity is 

| equal to the velocity of variation 

| in the area ao swept by radius vector 

| r. Its dimension is length?/time. 


The vector of sectorial 
velocity is perpendicular to the 
plane containing the vectors r 
and v, that is, the areas of the 
trajectory of a point. From the integral of areas it follows that! 


Fig. 3 


when a point moves under the influence of central forces the /20 


sectorial velocity is constant in magnitude and direction just 
like theconstant of areas. In the particular case, this assertion 
is valid also for the force of attraction. 


Probiem 3.1. Show that computing the integral of areas (3.1) 
is sufficient and necessary for the motion of a point to occur in 
a single plane passing through the center of attraction. 


Solution. Let us place at the attracting center O a rectan- 
gular system of axes xyz of arbitrary orientation. Let the com- 
ponents C in these axes be denoted by ec c._, and Co. that is, 


x? 
C = ec .{c,, Cus c,}. Then the vector equality (3.1) can be written 
as three scalar quantities: yz - zy = Cy» 2% - XZ = Cys and 


xXy - yx = Che Multiplying trese equalities by x, y, zZ, respectively, 
and adding the results, we arrive at the equation 


C,2+CyyY +€,2 = C:F=Q0, (3.3) 


which is the equation of a plane extending through the origin of 
coordinates, that is, through the attracting center, and perpendi- 
cular to C. For a central force of attraction these results mean } 
that the motion of a planet occurs in an unchanged plane extending ~ 
through the center of the Sun. This factor is reflected in Kepler's 
first law (see Chapter Two). 


Problem 3.2. Using the vector formula of sectorial velocity 
3.2), find the law of variation of area o swept by a radius vector, 
with time. 


Solution. From vector equality (3.2) follows the scalar 


equality 2S ste , and integrating it gives us o = s ct + oy: 


This law of linear increase in area o is written for an 
arbitrary central force. For the particular case of an attrac- 
tive force, this result corresponds to Kepler's second law (see 
Chapter Two). 


Problem 3.3. Find the sectorial velocity do/dt and constant 


of areas ec of the elliptical motion of a point acted on by a central 


force. Determine the values of the same quantities for the motion 
of the Earth around the Sun, assuming the Earth's orbit to be 
circular. 


Solution. When a point moves in an ellipse acted on by an 
arbitrary force its radius-vector sweeps out a total area of an 
ellipse nab in one period of revolution T, so that 


4g. ab , c= 2hab (3.4) 


In particular, this result is valid also for the case when a point 
moves in an ellipse acted onby agravitational force of attraction, 
when the attracting center is at .a focus of the ellipse. Taking 
the Earth's orbit: as circular, we have a = b = 150°10© km, T = 365 


14 2 


so that do/dt = 1.94°10 astronomical unit?/ 


/day. 


um /days = 0.86°10- 


Problem 3.4. Using the integral of areas determine the 
relation between the planetary velocities at the perihelion and 
aphelion (Fig. 4). 


solution. The vector equality 


rx v= eis satisfied for any 

point on the orbit, including for 
the apsides (perihelion «+ and 
aphelion oa): | 


PROU=I, KUy = Ta % Ug = Cc , 


where the modulus of the constant 
Fig. 4 of areas c can be determined by 
the formula 


c=rusin(F, 3) >0, | (3.5) 


-~-AQ L — — 
where (7,0) is the angle between vectors r and v at an arbitrary 
point M (Fig. 4). The quantity ¢c is always positive, so that the 


J2l1 


es 
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angle CP y does not exceed_180°. At the apsides the vectors 
r and v are orthogonal (sin (r 7 v) = 1), therefore the following 
relation is valid 


Ux, Ta 


Ya Tx (3.6) 


which means that the velocity of a point at. the apsides is inversely 
proportional to the distance between them and the ellipse focus § 
(that is, from the Sun). At the perihelion the planetary velocity 
Will be in the greatest, and at the aphelion -- the smallest. 


Remark. In generalizing these formulas to the case of any 
conic section and considering the orthogonality r and v at the 
pericenter of this section, we can write: for the ellipse c = 


= Tain “max ~ “max “min? and for the parabola and hyperbola 


c=r "Vy : 
min max 


For an arbitrary conic section the following formula is valid: 


sin (Fe v) - a = x7 max ; 


Problem 3.5. If a space rocket at an altitude 230 km over {22 
the Earth's surface is given a velocity 10.00 km/sec parallel to 
the Barth's surface, its orbital apogee will be roughly 370,000 
km from the Earth's center (near the orbit of the Moon). What 
velocity will the rocket have at its apogee? 


Solution. The computations can be made directly by Eq. (3.6). 
Assuming the mean radius of the Earth R, = 6370 km, let us find 


é 


= ge _6370+230_ __ 66,000 _ 
Y= nw = 10,00 F3%0+3r0000 = Sresro ~ 9-18 km/sec. 


Problem 3.6. Two meteorites describe the same ellipse, at 
focus S of which the Sun is situated. The distance between them 


is so small that the are M,M. of the ellipse can be assumed to be 


a segment of a straight line. We know that the distance MMs is 


« when its midsection is at the perihelion 1 (Fig. 5). Assuming 
that the meteorites will move at equal sectorial velocities, 
Getermine the distance MM. when its midsection will pass through 


the aphelion. The distances re and Yr, are known. 


Solution. Let us denote the 
symmetric positions of meteorites 


near the aphelion by My and M4. 


The condition of smallness of dis- 
tance between the meteorites 
enabling us to approximate the 

- are of the ellipse with a chord 
lets ‘us ; use, instead of the formula of the area of an ellip- 
tical sector, the formula of the area of a.triangle. In the time 
the meteorite M, transits the arc MMi; its radius-vector sweeps 
out the area M,SM}. On the condition that the sectorial velocities 
are equal to each other, this area is equal to the area MoSM5 


swept. out by the radius vector of meteorite M, as the latter 


transits the arc M,M}. From the equality of these areas there 


follows the equality of the areas of the triangles M,SM, and 


MjSM3, that is, LearedL(Mmay’ , from whence we have m/Mf ea je.! 


Since r,| <r, close to the aphelion the meteorites prove 


to be closer to each other than close to the perihelion (we have 

in mind the symmetric positions). Here their distance from each 

other near the aphelion proves to be at a minimum of all possible 
distances. 


Problem 3.7. Describe in polar coordinates the equations 
of motion of a point acted on by a central field and obtain the /23 
first integral of the equations of motion in the scalar form. 
Show that the resulting integral is an integral of areas. 


Solution. The equations of motion of a point acted on by a 
central force F = FY. in the polar coordinates r and ¢@ are two 


equations in projections in the radial and transversal directions: 

fxm (For @°)y Fy =—Peeae- (7? @)=0. (3.7) 

From the second equation at’ once there follows the first integral 
r°o =¢ 

se (3.8) 


To show that it is identical with the integral of areas presented 
above, let us show that the constant c from (3.8) is in fact the 
constant of areas. Actually, by turning again to Fig. 3, for the 
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area OMM' swept by r in time dt we get area OMM' = do = 5 r“do, 
from whence there follows r°4'= 2 Sts ec. The identity has been 


proven. 


Problem 3.8. Using the notation for the integral of areas 
in polar coordinates (3.8), derive a formula for the radial vn 


and transversal v5 projections of velocity v of a point as it 


moves in a conical section (Fig. 4). 


| Solution. When a point moves in a conical section r = 
=p/(l tecos 9), this problem consists of determining the quantities 


Uv, = r oe »iand vayue+ve e. 


ar ar. ad 
Y= gE "ag ae’ 


We obtain the derivative drfde = Faso from the equation of 


the section, and the derivative dofdt-efr*, from the integral of 
areas. Finally, we have 


V= S esiny, Vo “5 (i+e cosy), v= pF (1+ e*+2e cosy) ‘ (3.9) 


From (3.9) there follow the formulas of elliptical motion for 
velocities at the apsides: 


= ¥% (920) = 2 (1+e),v,=0, ((=180")= £-(1-e). (3.10) 


Knowing the~ratio of velocities and radius-vectors at the apsides 


(see problem 3.4), let us find 
‘Us Ts . {re | 
Vs Tx {-e (3.11) 


Substituting into Eg. (3.10) the expression p = a (1 - e*) and 
the constant of areas from (3.5) c= roy, = Pyvy> we again get 


the familiar expressions (see problem 2.2) ro =a (1 - e), ro 
: . a TF 
=q ({(1+e), ande= a 
a 4 
Problem 3.9. A point M with mass m moves around a fixed 
center O under the effect of a central force that depends only 
on distance OM = r. Knowing that the velocity of the point 
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will change according to the law v = a/r, where a_is a known con- 
stant, find the magnitude and direction of force F and the trajec- 
tory of this point. 


Solution. -To determine the force F let us use the first of 
the formulas (3.7) F=F, =m (rf - r¢2), in. which we must substi- 


tute the values of the quantities r_and ro2, To determine them, 
let us use the integral of areas r* = ¢ and the formula for the 


. ° ° ® y* = pry yr? 5? — 2 | 
expansion of velocity in polar coordinates rg = (a/r) - | 


2 | 


o 48 . ° 22 2 -2. ater | se 8 C%-@ 
Writing the relations g=ta > r@ = ty; p?= Saf, and F= | 


and substituting them into the formula for force F=m [Sag - £2). m 2 
= 

let us determine that the force F proves to be the attractive 

force inversely proportional to the cube of the distance to the 

points. Obviously, this force is not the force of Newtonian 

attraction. 


We can find the trajectory of the point from the relations 


~  fetae® | rar | ab ac. riage a dt, 

Ps ra *| - Yateer dt, r®°@=c, and oc 2 = d " from which 
. dr_ © dw. . 

after canceling out dt there follows arg = ¢ oe! By inte- 


erating the left and right parts of the equation, we get 


c 


ahr “9 +0 ,| where C =—patag Ine. | 


Hence it follows that the equation of the family of logarithmic 


spirals rene | is the equation of the set of trajectories. 


25 


26 


CHAPTER FOUR 
BINET'S FORMULAS FOR CENTRAL FORCES 
Velocity formula. Suppose a point is moving under the 


influence of central forces. We know that in the polar coordi- 
nates r and‘? the velocity of the point is expressed by the 


formula 
vis ult vy =($5) + r(agy. 


Transforming the expressions for v,, and v, by means of the integral 


2 de ° 
of areas r° gf 2¢,, we find 


and 


Let us introduce the new variable u = 1/r. Then considering that 


al du. _ 
gate gs we have ¥,=-¢ Te’ and = “w= eu. | 


Finally we get the formula of the square of the velocity of the 
point moving in the central force field, 


vis c( gu)” + ut], | (4.1) 


which is usually called Binet's first formula. 


Acceleration formula. Writing out the theorem on the change 
in kinetic energy of a point moving under the influence of a 
eentral force and dividing both parts of the equation by d¢, we get 


m , d(v*) _ Zr, 
2 ae rag 


™ 
in) 
I= 


! 


™ 
N 
Wh 


| 


Let us replace ve by its value from Binet's first formula: 


me*  @ (2 


o- f[ (44) +e? ]-6-$. 


Using the formula 


let us write the expression 


2 a 
me"fo du d?u du] _ 
a [2 de dg? + eu gu] ap 


w, = --c7u* (Sty +u). | (4.2) 


This expression of the acceleration of a material point moving 
under the influence of a central force is called Binet's second 
formula. Here the vector of acceleration is defined thusly: 


bd 


G = O, = w,Fo = EP, D, =O. | (4.3) 


Problem 4.1. Using Binet's formulas, find the law of action 
of a central force in which a point is moving along a circle with 
radius r = R. Determine the dependence of force on velocity. 


Solution. From Binet's second formula (4.2) there follows 
the expression of the force 


2,,8 _° ce 
fo a-Me*Uu -- US = const , 


that is. constant in magnitude. From Binet's first formula (4.1) 
we find 


A > | | 
so that 
mew 
fo- B = const, 


Thus, the motion of a point along a circle originates under the 
action of a magnitude-constant attractive force F. at velocity 


™ 
Nh 
ON 
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v, also constant in magnitude. This force is called the centri- 
petal force. The acceleration’ corrésponding to it is also constant 
in Magnitude and coincides completely with the centripetal component. 


Problem 4.2. <A central force causes a point to move in a 
logarithmic spiral r = ry eA? (4 = ctg a} (Fig. 6). Find the 
law of action of a force by using Binet's formulas. 

Solution. Based on the 


equation of the motion of a point, 
let us set up auxiliary equations: 


u- + e4 
! Bu =— 2 9 - hu, -Gity = 
: =? OO ay 
Fig. 6 which must be substituted into 


Binet's second formula: 
2 | 
Fp m= me7u> (A741) — —meoi(ctg ast) Q, | 


so that the force F_ is an attractive force directed toward the 
asymptotic point of spiral O, which is a dynamic force center. 
This force is proportional to 1/r3 and is among the forces of 
the same type as those examined in problem 3.9. Forces of this 
type are encountered, for example, in the theory of motion of 
microparticles. 


The velocity of a point defined by Binet's first formula 
varies inversely proportional to the distance of the point from 
the center 


v=cy/(qu)"s ur =cuV r+] = Lyf ctgaatt, 


so that the velocity increases with approach to the center and 
vice versa (the direction of velocity can be arbitrary). 


Based on the assigned r and v, we can determine the constant 
of areas 


¢ = —-|2— = rv sina aconst. . 
a+] | 
¥dgta+) (4.4) 


This same relation follows directly from the integral of areas 


CalFxd| = rv sin (F.3)=const, : 


Problem 4.3. Using Binet's formulas, find the exact law fet 
of action of a central force directed toward the focus of an 
ellipse, by using the equation of an ellipse in the polar coor- 
dinates r = p/(l1 + e cos ¢), where p = a(1 - e¢). (Direct problem 
of dynamics for a central force.) 


Solution. We know that if a point moves according to Kepler's 
laws, the force causing this motion is a central force directed at 
the focus of an ellipse. Based on the equation of the trajectory 
of the point, let us set up auxiliary expressions: 


. A+ @ cos ( @ i 
u Op, Be ae Sing, tues cos @, 


which must be substituted into Binet's second formula to find the 
law of action of the force 
Boo= --me?u2 1 _ me’ || 
r p pre *. 


inversely proportional to the square of the distance from the 
force center. The "minus" sign once’again confirms that the force 
considered is an attractive force. Usually the constant c¢/p is 
denoted by uw, so that 


re (465) 


The coefficient of proportionality u can be readily determined 
from the expression p = a (1 - e*) and. Eq. (3.4) of the elliptical 


sectorial velocity in which %=-=aYt-e?: 


(4.6) 


27 _ gt 27h? 47a? - | 
f p pT =—“yE =const. 


Using Kepler's third law (see Chapter Two), we can easily estab- 
lish the physical significance of the constant u = e2/p. Actually, 


if a°/PT* is constant not only for one body moving in an orbit with 
given o and T, but is also constant and identical for all bodies 
moving around the same central body, then it is apparent that the 
constant characterizes the intensity of the gravitational field 
produced by this body. The constant u is called the gravitational 
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parameter or Gauss' constant. Each body has its own gravitational 
parameter independent of the gravitational parameters of the other 
bodies. The relation between the gravitational parameter and the 
universal gravitational constant f will be determined below, in 
problem 4.6 


Problem 4.4. Determine the acceleration of a point moving 
under the law of universal gravity at the moments it transits the 
pericenter and apocenter. 


Solution. Based on the law of acbF0g 
force determined in problem 4.4, F 
acceleration we have 


i {Zhe gravitational 
/; pro), so that for 


4d | ; 
° =~ cr | woo - — = : 
We W, = -S. M =- 357, and a =~ ppE 


‘Knowing that 


pza(i~e*)= a({+e)(1-e), Tr, =e (1-@), and r, =a(i+e),. 


we get 


wy» Sse)? 4, _ otUl-e2) we — fises?® | 
Pp? 9 oa “ps? LD on (43 e, 


Problem 4.5. Based on the law of action of a central force 
(4.5) (see problem 4.4), derive the law of universal gravity and 
find the relation between the universal constant of gravity f with 
the gravitational parameter u and the mass of central body M. 


Solution. Let us examine any two bodies, for example, 
Sun and the Earth, with masses Mo and M,. Each of these bodies 
has its own gravitational parameter Ue and u,. The force with 


which the sun attracts the Earth can be stated using Eq. (4.5): 


Fl>-tee (4.7) 


Tr 2 


and the force with which theEarth attracts the Sun in accordance 
with Newton's third law can be written as 


= TE (4.77) 


™ 
NO 


From the law of equality of action and reaction, we get 


LOM, jr? “Bala Jr, whence 
He _ ps - _- pl _, 
M, Ms eee Mp1 const ’ 


(4.8) 


where Hy My is the ratio of the gravitational parameter of any 


planet to its mass. Therefore, the ratio of the gravitational 
parameter of any body to its mass is a constant. It is called 
the universal constant of gravity, or the constant of universal 
gravity. Let us denote the constant of gpravity by f, and then 


Ho PS (4.9) 


These formulas establish the relation between the universal gravi- 
tational constant with the gravitational parameter and with the 

mass of the body. Substituting (4.9) into (4.7) and (4.7') {2 
enables us to write the law of universal gravitation for the Sun 


and Farth _ 
- _fMoM 
[Ff] = Elo, 
or for any other two masses m and M 


\F I = Mm. 


r2 


The gravitational constants are these: in the SI system 


f = 6.673°10711 m°/ke*sec®, and in the SGS system f = 6.673°107° 


em?/g'sec®. 


Problem 4.6. Knowing that the mass of the Sun Mo = jl. 97°10°° 


g, the mass of the Earth M, = 6. 10°! ge, and the mass of the Moon 
Ma = 1/81.5 Earth mass, ae?ermine the gravitational parameters of 
these celestial bodies using the universal gravitational constant 
f= 6.67°107° em>/g sec”. 

solution. The gravitational parameters of the Sun, Earth, 
and Moon, calculated by the formula uw = fM, are, respectively, 


Un = 1327°10° km>/sec“, He = 398 ,600 Km?/sec* , and p< |= 4900 km?/sec 
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Problem 4.7. Knowing the radius of a celestial body R and 
the acceleration due to gravity at its surface g, determine the 
gravitational parameter pp of the celestial body and calculate it 
for the Earth (R, = 6370 km, g = 9.81 m/sec?). 

O 


Solution. Compare the formula » = fM with Eq. (1.10) from 
problem 1.6. As a result we get 


PgR, (4.10) 


that is, this. formula lets us determine the gravitational parameter 
of any celestial body from the values of g and R. From this 
formula we can calculate the gravitational parameter of the Earth 
already calculated by us in problem 4.6: 


BPs = Ms = RZ = 398,600 km>/see* 


Problem 4.8. From 
" 
parameter yp of a celestial body and the acceleration g due to 
gravity at its surface if we know the ratios of its mass M and 
radius R to the mass Ms and radius Re of the Earth. Compute 
these quantities for the Moon, Venus, Mars, and Jupiter. The 
corresponding ratios are given in Table l. 


and g, determine the gravitational 


TABLE 1 


ae 


Mars 
Jupiter 


Solution. From Eq. u = fM = eRo it follows that the gravi- /30 
tational parameters of different bodies relate to each other as 


the masses of these bodies: H/H = M/M,, whence uw = 398,600: 


é 


*M/M km3/see“, that is, from the mass ratios given in the problem 


& 
conditions we can caleulate the parameters of all these celestial 
bodies. Also, from this same formula it follows that the accelera- 
tions due to gravity at the surfaces of these celestial bodies 
relate as follows: 


£- OTe 


Whence g = 9.81 M/M 4 (R/R4)* m/sec’. By performing the computa- 


tions, we get the following results (see Table 2). 


TABLE 2 


Problem 4.9. Assuming that the central force is a gravita- 
tional force Fi= — um/re, determine the trajectory of a point under 
arbitrary initial conditions by using Binet's second formula 
(Newton's problem, or the inverse problem of dynamics for a gravi- 
tational force). 


Solution. We know that Binet's second formula (4.2) in 
conjunction with the law of areas (4.4) gives a system of differen- 
tial equations from which, by knowing F = P (r, ¢, r, ¢, t), we 


can determine the law of motion of a point acted on by a central 
force. For the case when Fe does not explicitly depend on time. 


(the gravitational force is specifically in this class of forces), 
Binet's second formula is a differential equation of the trajec- 
tory of the point. From this formula, the second-order inhomo- 
geneous differential equation 


d?u/ dg +u =—F. [me*u*=p fe? , 


Gerives for the gravitational force. The solution to this equation 
must be sought for in the form 


w= Ts A cos(@+e), | (4.11) 


f 


where A and ce are the constants of integration (we can easily 
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verify that the proposed solution satisfies the equation). Con- 
verting in (4.11) from u to r = l/u, we get the formula 


c® 


= ——_—_##______, 
} + oe cos (Q +8) (4.12) 


which must be compared with the equation of the conic section 


rea deee (4.12") 


From Eq. (4.6) (see problem 4.3) it follows that p = o-/u, where 
wis the gravitational parameter, so that the numerators of both 


formulas under comparison are equal to each other. If we set the 


constant of A/y equal to eccentricity e, and the initial phase 
of e« equal to zero, the formulas can be assumed to be identical. 


The zero-equality of the initial phase means that measuring 
the polar angle 9 (at the focus of the conical section), usually 
called the true anomaly, proceeds from the position of the radius- 
vector corresponding to the angle 9 = 0. Since in Eq. (4.12') 
¢ = 0 corresponds to r = Pain? then therefore measuring the angles 
proceeds from the direction toward '!the pericenter, that is, to the 
point with polar coordinates (rs = Train? oq = 0) 

Thus, the trajectory of a point moving under the influence 
of the gravitational force, under arbitrary initial conditions, 
is a conic section, whose shape (ellipse, parabola, or hyperbola) 
depends on the initial conditions imposed on the motion. 


Problem 4.10. Using’ the equation of a conic section and 
Binet's first formula, establish the relation between the orbital 
eccentricity and initial kinematic characteristics of motion By 


and Vo: 


Solution. The equation of a conic section (4.12') allows 

us to write 
_ 4+ e00s - _¢ Git __ . 
us tee , ecosy- up-t, ae = Bp SING. 


(4.23) 


Expressing the derivative du/d¢d by means of Binet's first formula, 
for the initial conditions we can write 


-£ cng - (22) ~ " 
D sing, = (GH) = e- uG 


Replacing p by its value from (4.6), we get 


+ C* Va _—,,* 
esa a2 (BR a, (4.18) 


where the signs correspond to positive and negative values of /3e 
sin 6) where $) = « (Q,7) and by = € (wT, On). 


From the second equation (4.13) it follows that 


= aC” — pe 
PCOS, FE (4.15) 


(the sign of cos. is regulated by the numerator). 


0 


Taking the squares of (4.14) and (4.15) and adding them, we 
find the desired relation: 


e=y/1+ Gr (v5 - Dupe). (4.16) 


As follows from (4.16), eccentricity does not depend explicitly 
on the polar angle 6 (the initial value of the true anomaly), 


however for a given constant of areas c, the relation of b0 with 
Vo and ry = 1/u, is expressed by a formula deriving directly 
from (4.14) and (4.15): 


usc fe | (4.17) 


The law of sign selection here is the same as in (4.14). We note 


that the constant of areas c = vup of conic sections for an ellipse 
a | 

c= fue (1-e*) (e<1),. for parabola csy2pr (e-1),, 

and for hyperbola c=yYpa(e*-4) (e>1), whence the 


positive sign of c follows once again. 


35 


36 


™~ 
WN) 
NO 


CHAPTER FIVE 


ENERGY BALANCE AND VELOCITY ALONG A SPACE TRAJECTORY 


In examining problems on the motion of artificial celestial 
bodies or spacecraft, we must bear in mind that the shape and 
linear dimensions of trajectories (conic sections) determined by 
eccentricity e and by the radius-vector of the pericenter re 


depend exclusively on the initial launch conditions. In problem 


4.10, the formula relating e with Po = 1/U, and Vo was obtained, 


and this formula can be rewritten thusly: 
(5.1) 


Analysis of the energy significance of (5.1), in particular, the 
energy significance of the expression appearing in the parentheses, 
enables us to establish the dependence of the kind of motion on /33 
these conditions. 


Problem 5.1. A point with mass m is attracted to a fixed 
center under the law of universal gravity. Write the integral of 
the energy of the point. 


Solution. To write the integral of the energy of a point, 
let us employ the formula of the potential of the gravitational 
force (1.7) (see problem 1.3), writing it in the form U = um/r, 
where uw is the gravitational parameter. By selecting the integral 
of energy in the form T - U = h, let us find mv?/2 - um/r = h, 
where Hh is the constant of the energy of the point, or the total 
energy. Dividing both parts of the equation by Mm, we get 


p?- 3h = 2h = kh. = const, (5.2) 


where h denotes the doubled total energy of unit mass. 


Customarily, Eq. (5.2) is called the integral of the energy 
of a point moving in a gravitational force field. It enables us 


to find the energy: significance of vé - 2u/r, from Eq. (5.1). 


Actually, since the integral of energy is satisfied for any 
position of the point in question in orbit, that is, 


2 2 2h fF . | 
v?- 2p vg- =f =A =const, (5.2°) 


then obviously the difference ve 


of unit mass (very frequently it is precisely this quantity that is 
called the energy constant). Now Eq. (5.1) can be rewritten as 


- 2u/Po is twice the total energy 


@ = f+ rh . (5.3) 


Problem 5.2. Based on Eqs. (5.2) and (5.3) derive the formulas 
for the initial launch velocities needed for insertion of a point 
into orbit with conic section. 

Solution. By specifying in Eq. (5.3) the eccentricity e 
corresponding to different kinds of conic sections, by means of 
(5.2') we get the velocities determining each kind of section. 


We can readily obtain the formula of local circular velocity 
as follows. Setting e = 0 in (5.1), let us write 


Ue ~ 2p fr, =—p7/c’, ; 


whence 


vp= (peje?) + 2p] ro- 
To determine yp/e? , let us use the equations 
B =¢%p/ and rela = PA + € COS Y= P » ; 
whence 
Biet=$'/p =4fr,) and > HIT ae | 


In the particular case of the launch of a body from the { 
Earth's surface (r 9 = Re = 6370 km), the local circular and 
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local parabolic velocities are called the first and second escape 


velocities, respectively:* 


When yp 398,600 km>/sec* (see problem 4.6), v. = 7.91 km/sec, 


I 


} 
and v 


tl 11.19 km/sec. 


When calculating Vy and Viz we can use other formulas. Actually, 


by comparing Eqs. (1.10) and (4.9) we conclude that u, = f£M, = gR,, 
whence g 5 5 


v,-VGR; . U=V2 R, - 
hee OES (5.5) 


In the case when a body is launched at an arbitrary distance 
from the center (Py > Ry) we have 


Tie | fabs ,! 
Vivei - (=| “i.par © | (5.6) 


By combining Eqs. (5.6). and (5.5), we can set up other formulas 
relating velocities for the launch of a body at altitude H above 
the Earth's surface: 


| 
afte . [gR = 4. / Re _ / Ra a 
Ya.ci Vr RH VR HO NV RR | 


: (5.7) 
foe, __fegRd _, [Re - ia: ! 
Ya par *We ~ Ror = Uy RH 7 N13 Rae it *| 
Using Eq. (5.6), we can write /35 


); 


6 = + EO fy 2 

pe (Us 1. par 
eC The terms local circular and local parabolic velocities mean that 
by these formulas we can calculate the theoretical values of these 
velocities for any given r. The actual velocity of a point may or 
may not coincide with these values. In the case when these velo- 
cities are the actual velocities in circular or parabolic trajec- 
tories, we will use the terms circular or parabolic velocity. In 
place of the term second escape velocity we can use the term escape 
velocity. 


on the basis of which, allowing for all the above-—presented 
relations, we can set up the following table (see Table 3). 


Table 3 


Problem 5.3. The second Soviet space rocket had a velocity 
of 2.31 km/sec at the distance 320,000 km from the Earth's center. 
Considering that motion occurred in a conic section, determine the 
Shape of the trajectory. What velocity did the rocket have at the 
altitude H= 230 km from the Earth's surface? 


Solution. Let us use Eq. (5.2) determining the invariancy 
of the following quantities for this trajectory in the calculations: 


E(B) (oe) 


Pe 320000 r°¢600_ . 


Let us determine the constant of energy from the velocity: 


2 


~ . km 
h = (v' Ee) = 2,84 —— >0, 
T 2 


whence it follows that the trajectory is a hyperbola. Further, 
we can determine the velocity of the rocket at the altitude 
H = 230 km: 


, 2} 3 | = ) km yy 
i / ue RetH N,l2 Sec - 


For comparison, we can calculate the local parabolic velocity at 


this altitude: 
v =f 2ns/R*H) = 10.99 SE. 
l.par ° sec 


39 


40 


The slight difference hetween the actual velocity and the theo- 
retical parabolic velocity enables us to state that the hyperbolic 
trajectory is nearly the same as a parabola. 


Problem 5.4. Determine the first Vy and second Viz escape 


velocities for the Moon, Venus, Mars, and Jupiter, by using the 
gravitational parameters and radii of bodies listed in Table 4. 


TABLE 4 


Solution. By calculating velocities using the formulas 


wos. 


v, - fz. | and %y > Y2uf/R = Yo 2, 


we get the following velocities (Table 5). 


TABLE 5 


Problem 5.5. To the first approximation, the orbit of the 
Moon can be assumed to be a circle with radius r = 384,400 km = 
> 60.4 R,. Determine the local circular V1 oq and local parabolic 


Vi.par 


velocities of a body relative to the Earth at this distance. 
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Solution. From Eqs. (5.7) we have 


W427 uv, V s0eRe = 7,91 / tox = 1,02 km/sec, 


= eo) ° = 1,44 km/sec. 
v1 par v2 Sei 


™N 
tad 
ze 


Problem 5.6. Determine the local circular and local para- 
bolic velocities of a point moving around the Sun in the Earth's 


circular orbit (r =a = 149.6°10° km). 


Solution. The velocities of the solar satellite can be 
calculated from expressions analogous to Eqs. (5.6): 


Po _ {397 -{0° . 
Moi Vr 7 Y teae-10F = 2978 kn/sec: 


| 


-~ YD v... =42,4 km/sec. 
V1 par v2 Lei ° 


Problem 5.7. Knowing the expressions for the radius-vector 
of a point in elliptical motion about an attracting center, 


F = pF° /(i+ecos ¢) =a (4 -ecosé) F°, 


where the unit vector po is the unit vector in the radial direction, 
E is the eccentric anomaly, and » is the true anomaly, find the 
expressions for the vector of orbital velocity of this point in 
orbital and inertial coordinate systems. 


Solution I (for the orbital system of axes). Any system of 
axes one of whose planes coincides with the orbital plane is 
called an orbital system of axes. In this case we can use as the 
orbital system a moving system of polar axes r’and 9% (Fig. 

7 a), all_the more so because the components of the orbital 
velocity v in the radial and transversal; directions are already 


determined (see problem 3.8). Using Eqs. (3.9), that is, 
y= Sesing, v= (i+ € 00s), | and c=ypp. | 


let us write the expression for the vector of orbital velocity in 
the orbital system of axes: 


2,05, VE esingr +()+ eco] 
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Solution II (for an inertial system of axes). The same system 


of polar axes pe and ge placed at a point which for this problem 

we Can assume to be fixed can be considered an inertial system of /38 
axes. It is convenient to select the pericenter of the orbit as 

this point (Fig. 7 a). 


a Let us write the formulas 
for transforming the axes of the 
moving system r® and $9 into the 
axes of the fixed system r® ~ and 
oo: 

re = Fy cosy + @ sin wy, | 


Go = OF cos@-F2 sing, | 


where » is the true anomaly of 
point M, that is, at the origin of 
the moving system. By substituting 
the formulas of transformation into 
the formula velocity obtained above, 
we get the expression for the 
orbital velocity framed in the 
inertial system of axes: 


° “5 [e sing (Fz cosy + G2 sing) +(1+ ecose)x 
(@2 cose - Fz sing) =YF- [ry singe@z(c ~ se) 


Key: 1. Vv 


Yi.ei To transform, use the relations 
of elliptical motion (2.4) and 
(2.5), whence 


, Vi-e’ sink ! 
sing =—~tre sine | 
a 1 -@ cos& | 
_ _fOSE -e 
Gos @ d-ecose * 


({-e*%) cos F 
€ + go = we eee 
* 8 @ {- ecosé 


so that 


™ 
WW 


- | ¥ {-e# sinE (1-e?) cos E zs | 
B= Ye ~£ = es = 


4~ecsf * 4-ecosE 


=f VA ef sine #2 + V4 ~e? cose Gt]. 


P {-ecosE 


When p = @ (1 - e“), we can write finally 
v «| ae [- on E Fr, +¥1-e* cosE ha 
vs jh. Vsin®E + (4-e *)cos?E _ fe fi+ecose 
a 1- ecose 4~eCOSE 


For the particular cases of the pericenter (E = Q) and apocenter 
(E = 180°), we have 


ven PAE Pas fe VS . 


Note that the latter formulas can be obtained directly from relations 
(3.10) (see next problem). 


Problem 5.8. Derive formulas relating velocities at the 


ae 


apsides v_ and v_; energy constant h, and total orbital velocity 


v with semi-major axis a of the elliptical orbit. 


Solution. In problem 3.8, based on the integral of areas 
the following formulas were derived: 


Uy = c(i+e)/p, vy=c(l-e)/ p, \ ) 


and in problem 4.3 it was established that p= c%/p.| so that c=ypp 
and c/p-v¥uj/p, whence 


v= /E JE | 
R 5, (4+e), v4 = (t-e). (5.8) 
In addition, for elliptical motion p=a(i-e),| therefore 

v= qe 5 ys fe ts», "| (5.9) 


and the formula relating velocities at the apsides is as follows: 


Va = vw, te, v, =u tte , ; (5.10) /40 
Vmin ® Unae Toe 3 Umax = Umin {-e - (5.10!) 
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From problem 5,1 it follows that the constant of energy h 
is defined by the formula 


ow 2 9 — a2 2 _ 
ha vi- SE = ug - Af = const. (5.11) 


Substituting instead of Vv. the expression from (5.9) corresponding 


to it, and instead of rie its value r_ = a(l - e), we get 


h Ete} 2 _ HU-e)  -. 
hs a(i-e) a(f-e) a(t#-e) «a | (5.11') 


From Eqs. (5.11) and (5.11') there follows the formula relating 
voand a: | 


teks Bean (G a): (5.12) 


Formula (5.12) is also customarily called the integral of energy. 


Problem 5.9. Express the velocity at any point of an ellip- 
tical orbit in terms of the eccentric anomaly. 


Solution. Let us substitute into Eq. (5.12) the formula 
for the radius-vector expressed in terms of the eccentric anomaly 
(see problem 2.1), r =a (l-ecos FE). As a result we get 


v= fe e-+) = /E PS - 


TaticemsE)_ fE fireest) 
a(i-e cof @ ~e@cosé | 


The same result was derived in problem 5.7 in another.way. 


Problem 5.10. Determine the total orbital energy of a point 
moving under the influence of a gravitational force as a function 
of the semi-major axis of elliptical orbit a. 


Solution. Let us write the equation of the energy balance 
of the point: 


= = _mv_ em omy 24\ mz | 
h=T+V2=+T-Ue= 5} +r - (vo fame. | (5.13) 


By substituting into (5.13) instead of the constant h, its value 
expressed in terms of a, that is, Eq. (5.11'), we find yy 


A= Meh =~ Eb = const. (5.13') 


Thus, the total energy of a point’ moving in an elliptical orbit 
With semi-axis does not depend on the radius-vector. At the peri- 


center the potential enerby V = - pm/r is at a minimum (the radius- 
vector r is at a minimum, and the potential U = um/r is at a 
maximum). Therefore the kinetic energy T = mv?/2, and this means 


also the orbital velocity v is at a maximum. But at the apocenter 
the potential energy is at a maximum and the kinetic energy and 
velocity. are at a minimum. 


If two satellites move around a planet in circular orbits, 
the one that is farthest from the planet has the greatest orbital 
energy. 


Problem’ 5.11. Set up a table of energy balance for four 
characteristic points of an ellipse by using the integral of 
energy. 


Solution. The following six points (Fig. 7 a) are called 
the characteristic points of an ellipse: the points of inter- 
section of the ellipse with its semi-major axis, or apsides, that 
is, pericenter1(¢ = E = 90°) and apocenter ao (¢ = E = 180°), 
points A (@ = 90°) and At (o = 180°), whose radius-vector is the 
focal parameter of the ellipse Pr, = Pp, =D =o (1 - e cos Ey) = 


=a (1 = e“), and the eccentric anomaly En = Bat = arccos e, and 


the points of intersection of the ellipse with its semi-minor axis 
B and B', whose radius-vectors are equal to the semi-major axis 
based. on the relations 


' ne intdszwuth_oe#l) 
Tr, * Ty 3 yY (ae)? +) j =¥are* +a (-e*) aa,: 


For the two last points the correspondin olar angle (true anomaly) 
¢ = arccos (-e), since sin ¢ = b/a = mere cos ¢ = -e, and the 


eccentric anomalies BE, = 90°, En. = 270°, 


Using the energy relations (5.12), (5.13), and (5.13'). for 
four (of six) characteristic points of the ellipse, let us deter- 
mine the quantities 


pa BRAS ye Be km 


and let us set up the following table. 
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TABLE 6 


Ne ere we ee ep ewe eee, 


~~ 


8 arecas-ay 90" Pa —! ~ + ~ 


9 4 
a. 189° 180 aci+e) Da(ire) ~Z+e) “Za. 


Problem 5.12. Prove that when a satellite moving in an 
elliptical orbit transits the end of its semi-minor axis B its 
velocity is equal in absolute magnitude to the local circular 
velocity V4 4° 

Solution. It was proved in problem 5.11 that the radius-vector 
of the end point of the semi-minor axis of an ellipse is equal to 
the semi-major axis (r, =o) (Fig. 7 a), therefore the orbital 


velocity of this point can be calculated by Eq. (5.12), 


vy Vp (2/m - Ya VELE 9} 


and this expression coincides with the expression for the Local 
circular velocity at distance r = a from the attracting center. 


Fig. % b'shows the elliptical and circular orbits corresponding 
to this case. These orbits have the same total energy, since the 


energy constant hy n-pfan-p/ry| is the same for them. Note 


that the equality of velocities can derive directly from the 
equality of the kinetic energies of these two orbits, since total 
orbital energy h and the potential energy V = - um/r for them are 
also the same. From Fig. 7 b it also follows that the velocity 
vectors Vi od and Vp form an equilateral vector triangle. 


Problem 5.13. Find in an elliptical orbit the points where 
the velocity is equal to the mean geometrical velocity at the 
pericenter and at the apocenter. 


Solution. In problems 5.7 and 5.8 we found the formulas for 
velocities at the apsides: 


ite =¥f vise: 
Ye =f t-e *, and a= i {+e 


™ 
= 
(ee) 


The geometric mean of these velocities. .is of the form 


V* os, vv, - Vila. 


However, in problem 5.12 it was shown that this velocity is 
attained when a satellite moving in an elliptical orbit transits 
the end point of its semi-minor axis where Pa = a. thus, both 


points desired lie at the end points of the semi-minor axis of 
the orbit (E = 90° and 270°). 


Problem 5.14. When a satellite moving in an elliptical orbit 
transits the end point of its semi-minor axis B (Fig. 7 b), its 
velocity Vp is equal in absolute magnitude to the local circular 
velocity Vi ai? $9 that the vectors Vp and Vi oq form an equi- 
lateral vector triangle (see problem 5.12). Determine the relative 
velocity Veer = Yi.ci 7 YB and express it in terms of the orbital 


eccentricity. 


solution. The vector of relative velocity is the closing 
leg of the equilateral vector triangle, so that 


2 C 2 ° 


Yrer = Yicea * YB @ Y2,c1%B COS a2” 


_ _ o _ _ 
L.ci Vp) cos (360° .- y) cos y 


= cos (¢, - 90°) = cos (90° - $,) = sin $,- 


where v = Vp and cos (va 


ei? 


It was established in problem 5.11 that the true anomaly of the 
end point of the semi-minor axis 4, ~ arccos(-e), so that 


t 


sin &, ~ sin[arccos (-e)| sy¥{-e, 


Substituting this expression into the velocity formula, we have 


= 72 V1 o4 U-G-et)*]" 


rel rel 


2 2 ) 
v = 2 vo a4 (1-v-e*),I r 


Expanding in a series, we get an approximate formula of relative 
velocity: 


i 


Yrel 7 Yucca e(t+$). . 
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The angle at the base of the vector triangle is determined by 
the formula 


y* g, - 90" = G,-=,- 


Problem 5.15. At some point of an orbit at distance r from 
the attracting planet, its satellite has the velocity v. The 


local parabolic velocity at this point is V1 par: Determine the 


semi-major axis of the elliptical orbit a. /44 
Solution. The integral of energy can be written in terms 
of a local parabolic velocity as follows: 


he yr - — Pi = ye _ y? 3 
r l.par 
where h = - Hy ts so that 
U 
go = —  — Pi 
ye 2 
l.par ~ ¥ 
If Ho is unknown, it can be determined from the known r and 


_1 2 
Yi.par’'p1l ~ > Yi.par *° 

Problem 5.16. Upon transiting across the pericenter, the 

satellite of a planet has the velocity Vee Its local circular 


velocity at the pericenter is v Determine the eccentricity 


1.ci a 
of the elliptical orbit. Solve this same problem for the case of 
a satellite transiting the apocenter. 


| Solution. By writing the integral of energy ve = u(2/r - 1/a) 
and using the formula of the radius-vector of the pericenter, we 
can state 


ve oh Be SR BOS) Kye) =. 2 


x Vs l.ci (1 + e) 
Whence . 
Q | 
@ = = ~{. 
Viet ul 


Since for an elliptical orbit Q < e < 1, then we have 


UR 
Lame <2 6. Ve <v z= _ 
Oise Ss Te “ois <Ms ri UY part | 


Thus, to determine the eccentricity it is sufficient to know the 
Vo and v at the pericenter. 
Tt 1.ci 
We can easily show that the formula obtained above for the 
eccentricity can be derived directly from Eq. (5.3): 


aw a Vx) — 
e=ylt-ark = jo Abate (un -B) = 


Vz 


V-wa) +a) 


Since v_ > v the expression in the parentheses is negative, /45 


L.ci wr’ | 
and therefore To ensure e > Q we must select the formula with a 
"minus" sign, then we.have 


2 
a 2% ~{. 


v 
1.ci*® | 


When the satellite transits the apocenter, we have 


ae ee ea a 
whence 
ent- er 
so that 0 < - <1,. that is, 0< Va < Vicia: 


Problem 5.17. The greatest distance of Sputnik-3 from the 
surface of the Earth H = 1880 km, and the smallest h - 230 km 
(Fig. 8). Determine its velocities at the apogee v_ and at the 
perigee v_. “ 
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Solution. From Fig. 8 the geometrical relations ©. 


2a =AR, tH th, 


ask, +$(H+h) = 6370+ 
+ $ (1880 +230) = 7425 in. | 


are obvious. The eccentricity can be determined by Eq. (2.13°): 


| : 
| | v, = /He y4=8 = SBR / thet = 6.06 see 
: | 


Fig. 8 


@ = Ya -lx _ (Re tH)- (Ry +h). 
Paths "(Rath +(Reo th 


H- Wy {650 = 0,44. 
~ OR, +H+h 14.350 


| 
| 
| 
To compute the velocities at the 
apsiacs apply Eqs. (5.9) and 
> 


_m tte . ted _ km 
Dy = Vee = 6,56 0,89 8,13 see" 


Thus, the velocity at the perigee. 


is 25 percent greater than at the 
apogee. 


Problem 5.18. Calculate the velocity of a satellite whose 
orbit has a large eccentricity at the perigee vo and at the 


apogee Vy The largest distance of the satellite from the Earth's 


surface Hi and the shortest distance h are 42,450 and 252 kn, 


respectively. 


Solution. Calculate a and 


e, as in probiem 5.17: 


~&=Ro+ 5 (H+h) = 6370+ £ (42450 +259) -27701 km, 


C= OR THT *Behao = 676, : 


H 


vas PER YTS = YSPBGM JOE. 1,40 | km/sec, 


Ve a Ve 142 - {+ 0,76 
Ug = Va Gre = 140 4 Otae = 


In this case the velocity of the perigee is more than seven times 


the apogee velocity. 


= 10,27 ‘gn/sec. 


J 
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Problem 5.19. Find the largest and smallest velocities of 
the Earth around the Sun. Compare the velocities obtained with 
the Earth's mean heliocentric velocity. 


Solution. In problem 2.3 the larges (aphelion) r, and smallest 
(perihelion) r distances of Earth from Sun were caleulated: re 
= 147.1°10° km, and ry = 152.1°10° km, and it was also shown that 
the arithmetic mean of these distances is equal to the semi-major 
axis: ro = 1/2 (r_ + ry ») = 149.6°10° km. In problem 2.5 the 


mean integral value of the heliocentric distance of the ERarth 
[r], was determined. Using these values, we can calculate the 


quantities we seek. Thus, the orbital velocity of the Earth 
reaches its largest and smallest values at the perihelion and 
aphelion, respectively: 


‘ 2 4 $ 2 { 
YVR (R> &) = 1182710 aatriae ~ Tag) ~ 80428 km/sec, 
\ 
_ Qo t\_ Jan? 2 _ : \ 
Vay Po (Fr zx} = y 1327 10 (estar HEB TIE) = 2% S3EM/ sec. 


(we can find v, also by Eq. (3.6) from problem 3.4). The arith- /4 


(5.14) 


metic mean of these velocities is: 


~4(v,+0,)= + [pH (2-4) +p. (2. -)]= 28,00, km 


sec" 


<j 


av 


Now let us determine what we must adopt as the mean integral 
velocity [v]. As already pointed out in problem 2.5, averaging 
can be carried out with respect to any variable used in the theory 
of elliptical motion: with respect to time t, mean anomaly M, 
eccentric anomaly E, or true anomaly @ (see problem 2.1). 


Let us first look at time averaging, when 


Lv], -4+Jodt, 
0 


where T is the theory of the Earth's revolution around the Sun. 
Adopt as the new variable of integration the eccentric anomaly E, 
since in problem 5.9 we derived a convenient formula relating Vv .to 


E: 
=y ytrecost . 
U= + 4-ecose 


a1 
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From Kepler's equation (see problem 2.5) we have 
— -esin & = 
en(t-T)=M, d&-ewsEdE = ndt=aM, dt = 
: =(f-ecosE)dE/n,i 


so that we finally get the elliptical integral of the second kind 


[v), =e JE lyme ease yt} -e? cos?E de 1B \U-e cas* £ - 


2x 
- + e°cos" E+, dé =f JE (E-F ef -F e? sin2& + eeten} : 
0 


- 2% /E (1-Le*+ 3 e*4.,.), | 


Since T = 2n/n, we find 


[v},=/e (I- gers a e*4...). 


We can easily show that the same result can be obtained by 
averaging over the mean anomaly M, which is an analog of time in 
elliptical motion. Actually, since ndt = dM = (1 - e cos E) dE, 


2z 
{r], = ay \uan . she fe \Vi-eF os Ea = VE () - 
-te +e e*+,..). 


On the other hand, averaging over the eccentric anomaly E 
leads to a different result: 


ex 2m | 
{ 
(v), ag JUdl = de Ea dé | 


_ E ((1secme)o otenstes ak = VE (+ £ eP+..), 


Averaging with respect to the true anomaly 9, we have (see problem 


2.5, where do = a/rv¥l = e% dE) 


@x 


4 7 i+ecose Yi-e* d& | 
[v], dz \vdq= & E | (eRe ecos—E [{-ecosE) ~ 
on 


=a f= 1-e teeny cos 2 E+...) E = 


“VE (1- $ ete Ses.) VE (14 Bete). 


° 


Thus, it is difficult to give an exact determination of the 
Earth's mean velocity, however considering the smallness of the 
eccentricity of the Earth's orbit (e = 0.01678), with adequate 
accuracy we can adopt as the mean Earth velocity the first term 
in the expansions obtained above, that is, [v] = Yu _/o a! O« In other 


words, the mean Earth velocity (or the velocity of another planet ) 

can be assumed, with adequate accuracy, to be the velocity that it 
would have if it moved uniformly in a circle with a radius equal 

to the semi-major axis of its orbit. This definition of the mean 
velocity of the Earth coincides with the definition of its helio- 
centric circular velocity as the Earth transits the end point of /49 
its semi-minor axis at a distance a from the Sun (see problems 

5.12 and 5.13). Moreover, we note that the quantity Lv] = Yu fo 


is simultaneously the circular velocity of a point at distance 
= ir, = o@ from the Sun. On. the other hand, the mean velocity 


can be defined as the heliocentric circular velocity of a point at 
distance Lr], = Lr Jy: 


1-85, * Vatey Ve 


For small e, as before we can use the definition of the mean 
velocity presented above. 


¢ -t¢e* + ae. = 


Thus, by agreeing to call the quantity [v] = vu the mean 


Earth velocity, we note that it does not coincide with the arith- 
metic mean of the velocity v_ and Vos To derive the Formula | 


relating [Lv] with Vo and Vos “Let us substitute in (5.14) Iv}° 
y/o: 
© 


[uv]? = Spe _ v2 _ 2p, _y? | 
x Vn ~* (5.15) 


We-can set up simpler formulas of the relationship for small 
eccentricities. Actually, from (5.14) we have 
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x a Px: 1-eé 
e_ Me, 2a-n p, ry i-@ (5,16) 
Ma sng See Se be: 


Hence, for small e we have 
v. ® [v]}° (442 2 [y]? 
x e), %, = [vu] (1-2e), 


vy, = [v] (1+e), v, ~ [uv] (f{-e), | (5.17) 


[v] ~ ¥, 4-e) = v, (14e). | 


\ 


Problem 5.20. What minimum initial velocity must be important 
to a spacecraft moving parallel to the Earth's surface at altitude /50 
h = 230 km for it to reach the Moon at its apogee, by moving in an 
ellipse that is tangent to the Moon's orbit. The semi-major axis 
of the Moon's orbit @&j| = 384,400 km, and the eccentricity of the 


Moon's orbit ey '= 0.055. Solve the same problem for the case of 
arrival..at the Moon when it is at the apogee of the given spacecraft. 


Solution. Calculate the apogee and perigee distances of the 
Moon: Tae= Ee (1 ee,),= 405,500 km, and Mequ @(4-&¢) =!) 363,000 km. 
To reach the Moon situated at the apogee a of its orbit 


(Fig. 9), the spacecraft must move in ellipse I (or It) tangent 
to the orbit of the Moon at apogee ag, and with its perigee at the 


launch point A (under the condition the initial velocity is hori- 
zontal). Here the semi-major axis of orbit I is 


ae (A+ 
+0d,)=3(% +a.) , 
where r.2R, +h: = 6370 + 230 = 6600 km is the perigee distance 
of the spacecraft, so that a, = 5 (6600 + 405,500) = 206,000 km. 


Let us find the initial launch velocity: 


| _ | 
v, “VPs (% - a) - (398800 (oq - 000) = 


= 10, 90 ' km/sec. 


Using the equality fa: 7 Tae > we can determine the eccentricity 


of the transfer ellipse: 


_ Tar -%e _ Tae Ts _ _398900_ | 
CL Ta try tae Te BIO ~ FH9EB 


(the orbit is strongly elliptical). 


For the case of the spacecraft 
reaching the Moon situated at the 
perigee Tq of its orbit, we must 


consider the transfer ellipse Il 
| (or II’) (Ta, = reeds | In this 


case we have 


/51 
Fig. 9 a,=$(B0+ Om) = 4(n+4,,)=4 (66004565500) = 184950 km, 
Key: 1. Orbit of Moon | «7VPse-®- (oe ls = teal (0. “(i 
oe ASE: SE SER - 0 
By comparing the resulting values of Vn and oo with the 
corresponding values V1 par for Pa = — we find Va par = 


ja = 10,09 km/sec ; 


Problem 5.21. Two satellites with equal masses move in the 
same direction around an attracting center in coplanar orbits, 
one of which is circular with radius ry» the other is elliptical 
with radius-vectors of the pericenter and apocenter Po and Br ys 
respectively. Assuming that the satellites by direct docking 
connect with each other at the point of tangency of their orbits 
and move together in further motion, find the radius-vector of 
the apocenter of their new orbit. 


Solution. Since the semi-axis of the circular orbit I (Fig. 
10) a> = To» and the semi-axis of the elliptical orbit II a:.. = 


TI 
= 1/2 (To + 8Y 4) = 9/2 o> 


the velocities of both satellites at 
the point of tangency of their orbits, being simultaneously the 


25 
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pericenter w of orbit II and the new orbit III, are as follows; 


_,f = 2 “4 = {\_ 4 
v, VE, me Pr a * —(1-9)= s7#. 


Upon satellite docking, velocity change occurs. The velocity 
of a new satellite, which is the arithmetic mean of the velocity 
Vr and Viz (for the same masses of both satellites), is 


ny) 201+ SJE - S/F 


The semi-major axis of a new orbit is found from the integral of 
energy 


m, “hae 


whence 
= tH (o- 49).23 ~ 26, _ 413 ! 
ay Y% (2 2) 6° Fr, and Gag = lay rg. 
Knowing that 2a, “htt, ./ we find 
22 _ 49 
ig = (a3 71) % = Gee 


Fig. 10 


Problem 5.22. Derive the formula relating the angle aoe (FLD): 


to the angle of the true anomaly » for the parabolic orbit (Fig. 11). 


j 


_~ 


Solution. Let us use the relations of paraholic motion 
r =ip/ +cos g), v=] 


In addition, for.any conic section the formulas sina = ¢/TY,| 


and c=<V¥pp,; are valid, whence 


— 


sind _YHP (1 +c0sg) Vp 
_ ppap (i +cosg) 


_ fit cosg | 
Se = +c05 £. 


| 
: 


O16 


At the pericenter @=0, and a= 2, With separation from the 


pericenter » increases, and ack decreases, so that 
sina -cos(%-a)=+osf  a=t(n-g@). 
2. ’ 2 


This result is valid only for parabolic orbits. 


Problem 5.23. Establish the dynamic meaning of the integral 
of energy 
he-u*- op 


r 


b 


for parabolic (h. = 0) and hyperbolic (h > 0} motion. 
Solution. Let us rewrite the integral of energy in the form 
y= h+2pn/r, 
whence 


.. ~ P) ~ 
lim v?= v2 = lim (A + 4) = Rk = const, 
P= oo oo rT 03 m 


so that the constant of energy h (double the total energy of unit 
mass) can be interpreted as the square of the velocity "at infinity." 
Here the interval of energy can be written as 


° : (5.18) /5 


I a 
- ~4 ~ V e 
h- veiw v Ya par /53 


The quantity v_ is also called the excess hyperbolic velocity. 


af 
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In hyperbolic motion (h = 0), the quantity W-ve-0 means 


that the true velocity coincides wholly with the theoretical value 
of the local parabolic velocity at a given point of the orbit 

= Tt : ; , u : : 
(v = V3 par? and that "at infinity" a stoppage occurs ‘(y -¥% -0), 


In hyperbolic motion ( >0) | the inequality Ah =u= >0 | 


means that the true hyperbolic velocity exceeds the theoretical 


value of the local parabolic velocity (v > Vy par? and the point 


will move indefinitely long along the asymptote of the hyperbola 


at constant velocity Do. - vie=comst, the velocity Ven is the smallest 
possible velocity for the given hyperbola. Here the difference 
between Vion and v wili:tend to zero . (Jim of 4.20) and the 
theoretical velocity "at infinity" Von coincides with the actual 


hyperbolic velocity (lim v= v..) 


Problem 5.24. Determine the velocity at which a meteorite 
enters the Earth's atmosphere if its velocity "at infinity" 
Vv = 10.0 km/sec. 
CO 
Solution. Let us write the interval of energy in the form 
of expression (5.18), determining in advance that; 


. 2 2 


h = Vo = 100 km“/sec > 0: 


thus, a hyperbola will be the meteorite's trajectory. To determine 
the theoretical value of the local parabolic velocity upon entry 
into the atmosphere, we will assume that entry occurs at an altitude 
of about 200 km above the Earth's surface so that v = 11.0 km/ 
/see (see problem 5.3). Let us calculate the entry velocity 

of the meteorite into the atmosphere: 


y = v2 + Ur at - (224 =14.9 km/sec. 


Problem 5.25. Knowing that for a hyperbola e > 1 anda < Q, 
ana that the formulas of elliptical motion can be transformed into 
the formulas of hyperbolic motion by replacing a with - lol], let 
us write the formulas of the velocity of a point v and the constant 
of energy h for hyperbolic motion, bearing in mind the motion along 
the branch of the hyperbola whose focus is the attracting center 
C (Fig. 12). 


1.par 


Solution. Using the formulas of elliptical motion {54 
2 ~ 
ye=p(2 - +), h - Bb 


and replacing « in them by -|lal, 
for the hyperbola we find 


pe 


We note that the resulting 


Actually, 


v 5) ! 
v= p(F + FE) = mR ge tte 
Fig. 12 
whence, by taking the difference ve = ve - ve to be a con- 
, : cD l.par 
strained quantity, we again have limv=v_. 


P+ © 


Problem 5.26. For a spacecraft 500,000 km from the Earth's 


center, the theoretical velocity "at infinity" vo 7 4.00 km/sec. 


Determine the true velocity of the craft v and the length of the 
semi-major axis of its orbit |e], by first ascertaining in 
advance the shape of the orbit. 


solution. Determine the constant of energy: h = v4 = 16 


aD 
km“ /sec* > O. Now we can conclude that the craft trajectory is 
a hyperbola (Fig. 12). Let us find the semi-major axis of the 


hyperbola from the integral of energy 


so that 
laj=n0 =: 398,600/16 = 24,900 km. 


We can calculate the instantaneous hyperbolic velocity from the 


formula 


Oo _2 2 
vio= Vm Ft Vi.par? 


where 


V1 par =V2p,/r = ¥2-398600/500000 =1,26 m/sec, 


so that v = 4.19 km/sec. 


formulas (5.19) enable us again to 
obtain the results of problem 5.23. 


va p(S+rh) A= Be (5.19) 


™ 
WI 
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Problem. 5.27. Establish the formula relating the smallest 
velocity Vo and the largest velocity vo for hyperbolic motion 


(Pig. 12). 


Solution. Using the formula 
a a 
peevs- Shi ang wax =e 


let us obtain the relation 
v= utl- 2n/p.! 


Knowing that 
n/p = 1/(R+1) 


let us determine 


fe-t | fe-4 
U5 = Uy a ? | oT Onin = Vier a i (5.20) 


(compare with Eq. (5.10') for elliptical motion). 


Problem 5.28. A spacecraft moving along the asymptote of a 
planetocentric hyperbola approaches a planet with gravitational 
parameter u "from infinity" at velocity Vip’ The distance of a 


craft from a pericenter is ri Determine by what angle 2’ the 
vector of the planetocentric velocity will rotate after the craft 
transits the pericenter. 


Solution. Since CB = An for a hyperbola (Fig. 12), AO = CO, 
so that 


_mo., 1%) 2 
cosA = 5 ~ laity” 4 t+re/lal ’ 
Yr ; 
where laf= te, 
ve 


Finally, 


2h =2arccos (it vi /y y* | 


Here the radius-vector of the craft will rotate relative to the 
launch planet by the angle y = 180° +-2\. If, for example, it is 
not r, that is known, but the distance of the craft from the 

planet CB. (distance of planet from the asymptote of the hyperbola), 
by denoting CB = Aw = d, we can write 


2) = 2arctg 5% = darctg a . 


Let us illustrate this problem by the following example. 
Suppose a spacecraft is inserted into an Earth orbit at the alti- 


tude v/s = 230 km with a hyperbolic launch velocity V Sd = 12.00 


km/sec. The orbit of the craft is situated in the orbital plane /56 
of the Moon. The approach of the craft to the Moon at the peri- 


selenion is m,{ = 6300 km. In this case, to determine the angles 
2x and y first from the initial data we must find the elements of 
the geocentric hyperbola 


2 yo 
Win Ps Fy, + jay + bal =17170 km je = 724. +4= BOO t= 1,384, 


The theoretical velocity "at infinity" must be determined from the 
formula | 


Maye Vh= Vpn Thal 


or from the formula 


Yosy™ Vy p, ¥E- Ve t ),. 
so that 
Vp /A = 4,82 km/sec. 


Since the direction of the asymptote of the geocentric hyperbola 
coincides at the limit with the direction of the geocentric radius- 
vector of the craft, it can be assumed that the absolute geocentric 
velocity of the craft Yon /$ is perpendicular to the translational 
circular velocity of the moon Uy,» Which, as was determined in 


problem 5.5, is 1.02 km/sec. The relative selenocentric velocity 
of the craft v,, is, in accordance with the theorem\of the 


addition of velocities, 


Vooyg =V Voy, t Ya, = Y (4,82)? + (1,02)° = 4,92 km/sec. 
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Knowing that Hf = 4900 km/sec’, for the angles 2A and y we 
get 


25-1 | 
2 = 2 arccos I + 8300.0) | = 2arccos 0,033 ~{76", : 


{ =180°-2)~%4°' The parameters of: the selenocentric hyperbola 


are as follows: 


| 6300 
el = “2 Ggaye 7") SM? Ta _ | 


and the selenocentric craft velocity at the periselenion is 


e+! 35 
pj * Yong ¥ G-T = 4,92 59,5 °, 5.08 km/sec. 


From the velocities obtained Vise 12.00 km/sec, 4.82 


4 “wo /é ~ 
km/sec, and Veoje = 4.92 km/sec, we can determine the minimum angle 
of rotation of asymptote 2\Aithat is possible for this problem, 


corresponding to the case when the hyperbola is tangent to the 
surface of the Moon ( ", = R, = 1740 km): 


aye . | | 
eh... = 2 arccos (1 +“oRe | 's 2arccos 0,41 = 167" | 


This angle corresponds to the maximum angle of rotation of the 
geocentric radius-vector yy * 13°. 


We can easily show that the values of the angles 2A and y 
lie at the limit 


[ 


2arecos (1+ B=) ¢ 2d < 180°, O< ] < 180°-2A, | 


where R is the radius of the planet. This means that even when 
tangency is involved the trajectory will not become a hyperbola. 
So in the case of no values of rp can capture occur for this 


problem, that is, there is no possible case when the hyperbolic 
trajectory converts into an elliptical trajectory. 


™~ 
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CHAPTER SIX 


TIME OF MOTION IN A SPACE TRAJECTORY 


The time during which a satellite twice successively transits 
its pericenter is called the period of revolution (orbital period) 
T of the satellite around the attracting center. The number 

n= 2 is its mean angular velocity, or its mean velocity (see 
problem 2.1). In problem 4.3 the relation of period T and the 
semi-major axis a with the gravitational parameter of central body 
u (4.6) was established: 


ps = 4n%*a/T 4 
The relation between the periods of revolution of different 


satellites around the same attracting center is characterized by 
Kepler's third law, established empirically: 


Ts — aa (6.1) 


In the general case of motion in a conic section, the time 
of motion of a body to an arbitrary point on the trajectory is 


determined by using the integral of areas r-> = c (see problem 
3.7), from which in particular we can get the familiar Kepler's 
equation derived in problem 2.1 from geometrical ratios. 


Problem 6.1. Derive formulas for calculating the period of 
revolution of a "zero" satellite T (r = Ry) s and also for the 


zZer 
period of revolution of an arbitrary Earth satellite T expressed 
in terms of the period of the "zero" satellite T oer’ 


Solution. From Kepler's third law (6.1) it follows that 
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7? a’ a \7F 
oer Re ° r= Tor Re) sf (6.2) 


Let us calculate Ter from the formula for uniform circular 


motion: 
Bek, _ 206370 _ 7 
Ter 7 op ab || «84.40 min, / 

so that 


3/2 


T= 84,40 (gaz) /min. (6.3) 


Problem 6.2. The perigee altitude of a satellite He = 230 km, 


and the apogee altitude HO = 950 km. What is the period of its 
revolution? 


Solution. By Eq. (6.2) we can write 


- oa 3/2 
r= (28 these) = To 1+ eae) = 
Ze c QR, = oR, | 


wT (j+ duets) . 3, 2304050), , 
Taft Big) = 66,40 (102 ARP) 96.8 nin 


Problem 6.3. Find the period of revolution of a satellite T 
and its mean motion n if the semi-major axis ao and the gravitational 
parameter of the attracting center uw are known. 


solution. Let us use the result of problem 3.3, in which it 
was established that when a point moves in an ellipse its radius- 
vector sweeps out the total ellipse area nmab:in one period of 
revolution T, so that T = 2tab/c, where c is the constant of 
areas. Substituting in place of c its value Yup from problem 4.3, 
we get 


ro ZOO  _axab_ _ ona 
VRP YE Vea Vp (6.4) 


. an . . 
since nm=>>» The mean motion is 


‘ - ~3/2 
mate aM (6.4') 


Using Eq. (6.4), for example, we can determine the radius of the 
circular orbit of a 24-hour AES: 


3fxT®  3/398600(24- 3600)" 3 
Gers oe = SD -/75 5-0" = 42,170 km . 


The altitude of the orbit of a 24-hour AES is 35,800 km (see 
also problem 1.7). 


Problem 6.4. Show that the period of revolution of a 
satellite T depends only on its initial velocity and does not 
depend on its direction. 


Solution. Based on Eq. (6.4) we can conclude that the 
period T is uniquely determined only by the magnitude of the 
orbital semi-axis. In turn, the semi-major axis a is associated 


with the algebraic value of the initial velocity Vo by the integral 


; a 72 4a\! , 
of energies Vv, = we (e - a): , so that the period T does not 


depend on the direction of the initial velocity. This means that 
if at some point in space a launch is made with elliptical 
velocity Vo (Vo og < Vy < V1 par)? for any direction of this 


velocity ensuring the existence of the flyby trajectory, the 
satellite will be inserted into orbit with the same semi-major 
axes a and, therefore, with the same period. These orbits can 
have different eccentricities. In any case the satellite will 
return to its launch point after the same time interval. 


Problem 6.5. Find the relation between the periods of revo- 
lution of planets around the Sun qT, and the semi-major axes of 
their elliptical trajectories a, (give the mathematical substan- 
tiation of Kepler's third law). 


Solution. To justify Kepler's third law let us utilize the 
familiar relation, which is of the form 


Ho a 
“Gut ~ 7% * 


for heliocentric motion. Since the constant Ue =°1327°10° 
characterizing the Sun's gravitational field appears in the left- 
hand part of the equation, for all heliocentric orbits the 
following condition is satisfied: 


Te 0 TR ARE (6.5) 


km>/sec* 


™ 
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or 


(6.5") 


expressing Kepler's third law. 


Similarly, for geocentric orbits we can write 


: 3 3 3 
7,2 ~ Ti ~ eee T? gg? (ps 398,600 sec’ ). 


Problem 6.6. The orbital perihelion of the first Soviet 
space rocket that became a satellite of the Sun is at a distance 
from itr, = 146.4°10° km, and its aphelion is at a distance 


ri. 197.2°10° km. Determine the period of revolution T of the 


artificial planet / around the Sun. 


solution I. Find the semi-major axis of the orbit of the 
artificial planet (AP): 


} 
at ad (ry +m) = $ (146,44 19762)+10% = 171,810 em, 


To determine the period of revolution of the AP let us use Eq. 
(6.4) sinee the parameter Me = 1327°108 km3/sec? is known: 


3/2 . | 
_  2na” | far _ /(171.8-10°)° _ 
Typ V He = 2m po - 6.28 1327-108 7 | 


=6,28V58,21-10" = 38,82-10"c = 449,27 gays 


Solution II. Knowing the orbital semi-major axis a Of the 
AP, we can use Kepler's third law (6.5'): 


JAP . ZAP , Whence when T, = 365, 2b aysa, =149,6-{0°: km 


is WEY = 369.26) = | 


= 365,26 ¥1,52 = 449,27 days. 


Problem 6.7. Using the integral of areas in the equation of 
28 conic. section, set up the general integral formula of the time 
of motion in the conic section as a function of the true anomaly 
(polar angle) ¢. 


ON 


Solution. Using the integral of areas in the form r°4 = 


= p* @ : . $ = | 
r’ Ffe-e¢!] and the equation of a conic section ” TEESE + 


we get j 


2 


adt= ridge __P FeEXOLE 
~ Cc Cc {+ecos¢ 


Let us reduce the coefficients p“/e by means of the relation °*VpEp | 
(see problem 4.3) to the form p@/Vy , and then let us inte- 


grate the expression for dt from 0 to » (from-the pericenter to 
the arbitrary points): 


4 
- t-t, = -— de 
te t-t, = Fo \ tha - 


(6.6) 


where t. is the moment of transiting of the pericenter. 


Problem 6.8. From the integral formula (6.6) and the 
relations of elliptical motion, by using the substitution 


ig Z -l ite ts 
(see problem 2.1), derive Kepler's equation for an elliptical 
trajectory. 


Solution. Based on the formula pea(l-e) |, the integral 
formula can be written as: 


y 


7 V8 (4_ 92)??? a | 6.6! 
t ee | tier: 0-0") 


Then, by means of the substitution indicated in the condition of 
the problem, we can convert to the variable of integration E 
(eccentric anomaly). To do this, we must differentiate the 
expression obtained by substitution, after which we have 


2Q@ . _ . to ne £ | 
sec £ dy yite sec’ 
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2 {+e 
Further, using the relation sec? $ a+ tg $ [=e 


we find the expression 


{+e | 
e=yi-e cos’ 5 + 4+8 sin? 
Employing the formulas 
cos*£=+(I+cos—), and sin? £ = 4 (1-cosé) ' 


we can write 


= e (l-ed— _ yi-e® dé 
ap = fe 


j-ecost  {-ec0sé 


We can also éxpress the integrand by E, by employing the 
relations of elliptical motion 


r=a(t-ecos£) jand rcos@ =a (cosé -e) | 


{-e? 


{+ €COS@ = 7—ScosE . 


(see problem 2.1): | 


Finally, the integral becomes 


c= aAiner)t {t-emee ae - =f lE -esinE). (6.7) 


Using the formula neypa*? we get Kepler's equation that we 


already derived in problem 2.1 by another method: 


F-esinE = n(t-t,)=M.. 


This expression served in the determination of the time of motion 
in an elliptical trajectory. It can also be solved for E for 
known t. In this case Kepler's equation is solved by the method 
of iterations. 


Remark. We can similarly obtain Kepler's equation for a 
hyperbolic trajectory by using the analog of the eccentric anomaly 
H = iE, which has real values. The corresponding formulas when 
e > 1] are of the form 


gd -/2+} tho, p-=lal(e?-t), 


reosg¢ sjaj/(e-chH), re=laf(echH-), 
eshH -H= n(t-t,) = yp lal (t-t,). 


Problem 6.9. A spacecraft is flying to the Moon in an ellip- 
tical geocentric trajectory tangent to the lunar orbit at its 
apogee %|(Fig. 9). Determine the time of flight of the craft 63 


by the apogee of the lunar orbit a_| (orbit I). Solve the problem 


for the case when the flight is made at the perigee of the lunar 
orbit h,_; (orbit II).if the perigee and apogee distances of the 


Moon are xd = 363,300 km and 7*| = 405,500 km, respectively. 
Solution I. In problem 5.20 we determined the major semi- 


axes of the transfer orbits I and II, a, = 206,000 km and on 
= 185,000 km. 


Now let us write Kepler's. third law.(6.5) for geocentric 
motion: 


GF ay ff, — 398600 | . km? 
TE Te GR WG Gaye LOND TG 


whence the total periods of revolution in orbits I and II are 


a? _ ,/ (2,06-10°)° -y To. 
“Vion?  ~ V1on-107 6,69: 10 


=9,30-10°sec= 10,76days /Odys 18,2 hr 


-// ai = {1085-10°) «1/6 26-10" v | 
Ty 7 4,044 -40 7 {,044-10% : ac " | 


= 7,91+{05ec 9,16days= GdayO3,8nr 


The half-periods we seek are, respectively: Tr = 5.38 days = 


= 5 days 09.2 hours, Try = 4.58 days = 4 days 13.9 hours. 


solution II. The quantities Tr and Try can be obtained 


directly from Kepler's equation (6.7), by substituting instead 


69 


{0 


of E its value E = 7m: 


2 

3/2 Te yy . 

- —bh- (f-esinE) = 4A = : 
tig / Ms 


3/2 -2 3/2 
~ SelG Fur = 9 49810 ayy 


~ ¥ 398800 
9 3f2 0.498 - ig g 40+40° =, 
Uy = Q, 4Q8& -40 ay ” a . 
=4,66:{Osec= 5,38days= Sdays09,2 hr 
‘ \ 
- - { 
T, = 0.498 « 10°@,""= 0498-402 7,98 -10"= | 
= 3,97: (Osec= 4, S8days= 4dayst3.9hr 


Problem 6.10. A spacecraft moving in an Earth orbit is given 


an additional velocity impulse at point A bringing it into a trans- 


fer heliocentric elliptical orbit I tangent to the orbit of Mars 
at point B (Fig. 13). Determine the flight time to Mars’ orbit. 
Solve the problem for the case of flight to Venus! orbit (orbit II). 
The orbits of the planets are assumed to be circular with radil 


r= 150°10° km, r, = 228°10° km, and ry = 108°10° km. 


1 


| Solution. The point of tangency 

| of the orbits B is the aphelion 

| of orbit I, and the point of 

tangency C is the perihelion of 

orbit II, since at the points of 

| tangency the vector of the orbital 

| velocity is perpendicular to the 

| radius of the circle of the orbit. 

| Let us determine the semi-major 
axes of the transfer orbits: 


Fig. 13 . a, = 2 (r+r,) = 189-40° ih, 
Key: 1. Venus! orbit { g 
2. Earthts orbit a= 3 (r +7) = 129-40 km. 
3 Mars! orbit | 


Knowing that the period of revolution of the Earth around 
the Sun is 365.26 = 365.3 days, let us use Kepler’s third law: 


Zz, | 
fat 365,35 // oats) » whence 7, = 365,3f (ee) =545,4/ days, 
y (422 3 oO 


TF, = 365,3 /(128)* = 291.2 / days, : 


so that the flight time to the orbits of Mars and Venus in an 


elliptical trajectory is tr = 257.6 days and Try = 146.0 days. 
Problem 6.11. A spacecraft will move in a geocentric trajec-— 
tory transecting the Moon's orbit (Fig. 14). The insertion into 65 


orbit occurs at the altitude H = 230 km at a horizontal initial 
velocity Vy = 10.95 km/sec. Considering that the Moon will move 


in a circular orbit with radius "™/= 384,400 km, determine the 
Flight time to the Moon's orbit. 


Solution. A comparison of 
the initial launch velocity Vo = 


= 10.95 km/sec with a local para- 
bolic velocity at the launch point 


V1 par Vassy 7) 10.99 km/sec 


(see also problem 5.20) enables 

us to conclude that the flight 
occurs in a strongly elongated 
elliptical orbit, since the launch 
velocity does not exceed the para- 
bolic velocity, but slightly 
differs: from it. Let us determine 
the semi-major axis of the transfer orbit from the integral of 
energy written for the perigee (the launch is horizontal): 


Fig.14 
Key: 1. Moon's orbit 


Us = MRK - 4). ! 


Hence from the given vy = v, and H = H_, we find a = 450,000 km. 
Let us determine the orbital eccentricity by the formula r_ = 
= (1 -—- e), whence we have 


e=1-Ty/a =4-AS5 = 0,985. , 


™~ 
ON 
On 


Thus, the ship executes the flight in a strongly elongated geo- 
centric elliptical orbit I, whose center C, as we can easily 
determine, lies beyond the Moon's orbit (Fig. 14). The time of 
motion in the are of the trajectory 7A of interest to us (from 
the perigee to the point of intersection with the Moon's orbit) 
ean be found from Kepler's equation (6.7): 


‘Ya 
t= F— (E -esinE). 


f2 


We can determine the eccentric anomaly of the point of intersection 
E in advance from the equality r = a (1 - 3 cos E) = 384,400 kn, 
whence 


a-r _ 450,000 - 384400 | 
cosE = Ge = —450000-0,085. 0,148 ? 


sin FE =0,989 , F 24.4299 = 84°29". 


Substituting these quantities into Kepler's equation, we have 


3/2 | 
t a (1,422 -0,985-0,989 ) = 2,14+10°/ sec = 


=59.4hr = 2,45 days. : 


Remark. We-can consider as the point of intersection also 
the point B, for. which when cos E, = 0.146, E, = 278°24', and 


B B 
sin E, = -0.989, so that 


2 


af 
450,000 
- 4 (4,422 + 0,985°0,089)= N,47°10°/ sec = 


V398600 
i 


= 3418, hr = 13,3! days. 


Problem 6.12. Solve problem 6.11 by assuming that the /6 
initial velocity vo = 12.00 km/sec. 


~] 


Solution. By comparing the given initial velocity with the 
local parabolic velocity Vy par = 10.99 km/sec, let us conclude 


l.par? hyperbola II is the flight trajectory 


(Pig. 14). By writing the integral of energy for the hyperbola 
(see problem 5.5) 


that since Vo > V 


vo = p, (2/7, + t/lat), 


let us calculate the semi-major axis: 


|| = IZT70= 2,7 Ry. 


From the relation ma lal(e-1)), we get the eccentricity 


ry, _ 66900 
er +{= FPO +4 =1,384 


To determine the flight time to the Moon's orbit 1t, let us 
use the formulas of hyperbolic motion obtained in problem 6.8: 


re=faul(echH-1), ! 
e shH -H Vp, la) (t-t,) | 


(Keplerts equation), which in our case can be written as: 


r, s|a|(echH-1) = 384400 xm, 
t =t-t= 7 (e shH-H). 


From the first relation we find 


echH = ito = 22,40+4= 23,40, 


ch H = 23:40 216,91, sh =/sh™H~1 16,88, | 


eshH =23,36, AH=3,52 = 201°44', 


which, on being substituted into the second, allow us to determine 


(17470) _ 3. _ | 
wT, * V sah600 (23,36 ~ 3,52)= 3,56+40 19,84 


f 


= 7,06 40sec = 19,6 ‘hr. 


Problem 6.13. Derive the approximate formula for the flight /68 
time in an elliptical trajectory as a function of true anomaly 
(polar angle) ¢ for small eccentricities. 


Solution. Let us use Eq. (6.6') (see problem 6.8) derived 
for the elliptical trajectories: 


13 


¢ 
_ ae)" (de 
ve VE J (f+ecos@)? - 


Let us expand the integrand in a series in powers of eccentricity 
e. Limiting ourselves to the- first power of e, we get 


e | 


_a” 2 d -2e sin ). 
tT ~ ie SC e cos ay = Ee (4- e sin g (6.8) 
; Vie -3/a | 
Since the mean velocity m= 20/T = Vp a + | we have 
T= xh (y-2e sing). (6.9) 


Problem 6.14. A spacecraft satellite inserted into Earth 
orbit has a perigee altitude H_ = 180 km and an apogee altitude 


= 340 km: The moment the ship transits the perigee t. = 9 hr 


OO min. In the next revolution it is required to fire the retro- 
engine at the moment when the true anomaly of the craft 9 = 270°. 
Determine the semi-major axis:of the orbit of the spacecraft a, 
its eccentricity e, the period of revolution of the craft T, and 
the moment of firing the retro-engine t*. 


Solution. From the given altitudes H and H, let us find 
Oy Ts and e: 


© = Ry +-5(Hy+H,) = 6630 km, ry = Ry +Hy = 6590 km 


but Tr, = a({-e), whence e oa-7 = Bo, #=Q,012. 


From Kepler's third law written in the form of Eq. (6.3') (see 
problem 6.1), let us determine T: 


_ a ye 7 $850) 
P= 84,4 (S320) 84,4 e370) = 89,5 min} 


Considering that the eccentricity (e = 0.012) is quite small, 
let us assume that the approximate formula of the flight time (6.9) 
obtained in problem 6.13: /6 


el 


T=t-t, ~ g-(g-2e sing), | 


whence we get for the time of flight from the perigee to the instant 


the retro-engine is fired (¢ = 3/2 7): 


270°. ae EL + 0,024)=66,7 min.' 


Let us determine thé moment of firing the retro-engine: 


t% = tT + + = 9 hr OO min + 66.7 min = 10 hr 06.7 min. 


Problem 6.15. From the integral formula (6.6) (see problem 
6.8) and the relations of parabolic motion, derive the formula of 
flight time in a parabolic trajectory from the pericenter am To 
the point A with given polar coordinates r, and ¢, (Fig. 11). 


Solution. For a parabola e = 1, therefore the integral 
formula (6.6) can be rewritten in the simplified form: 


3/2 Na 


ES: le lag 


Using the relation 
{/cos* $ = sty Se 


let us write 


3/2 ‘at 
©) = sige J 


¢ 


. ‘2 
soll+ h tg 4). 


Thus, t can easily be expressed not only as a function of ae but 
also as a function of Pas for which let us use the equation of the 


Q 
~ Te = (l+ty 2). 


whence 


oO 5 
tg-$=y -! 
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We finally get 
t(m)= Pa fet B+) = nc ae (6.10) 


In this formula parameter p can be replaced by its value p = ers 
characterizing the parabolic trajectory, so that 


© (7%) =sypee V(r - 7a) (+27). (6.10! 


Problem 6.16. <A rocket at an altitude H_ = 230 km above the 


Barth's surface is given a horizontal parabolic velocity. Deter- 
mine the flight time over a distance equal to the mean radius of 


the Moon's orbit | = 384,400 km. 
Solution. Let us use Eq. (6.10'), by setting r, = r=! 
= 384,400 km, ‘%7% + Ag i = 6600 km: 


V/2(384400 - 6600) 3 


c= ERE) (334400 + 15200)= 1.82°10 sec = 50.6 hr. 


It is of interest to compare this result with the results of 
problem 6.11 and 6.12: flight in an elliptical trajectory takes 
59.4 hours, while it takes 19.6 hours in a hyperbolic trajectory. 
The solution for the problem of flight in the parabolic trajec- 
tory is simpler, since it does not require the integral of 
energy. For this reason, we need not use the parabolic velocity 
in calculations. The local parabolic velocity for the altitude 
H = 230 km was already determined in problem 5.20 (v = 10.99 
km/sec). 
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CHAPTER SEVEN 


CONDITIONS FOR THE EXISTENCE OF ELLIPTICAL TRAJECTORIES 


Among the elliptical trajectories, in which the main condi- 
tions for existence is the velocity constraints 


Vici s vs ‘1 .par? (7.1) 
we will distinguish two classes of trajectories -- flyby and 
ballistic. For a body launched from the Earth's surface (ry = R,) 


to describe about it a closed elliptical curve (0 < e < 1), that 
is, to enter a flyby trajectory, the flyby condition 


r> ry =R (7.2) 


} 
-must ‘tbe satisfied (the case of a "zero" satellite when r = Py = 


wien amen og = 


~ dy 
Rye 


it also can be included in the class of flyby trajectories). For 
a launch from altitude H above the Earth's surface these conditions 
(the circular orbit is included) become 


= Vrs and e = 0 is not considered by us here, although 


(7.1') 


V . < < 
L.ci v Yi .par? 


an 
WwW 
an) 


O° (7.2') 


In addition to these obvious and necessary flyby conditions, there 
are yet a number of other conditions that will be established in 
the following problems. 


An elliptical trajectory which often occurs within the Earth 
(r < Ry? and motion along it is possible only outside it, that is, 


Uf 
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where r 3 Ris ts called a ballistic trajectory. Below we will f{1 
establish the criteria that enable us to determine whether a4 
trajectory belongs to a particular class of elliptical trajectories. 


Problem 7.1. Show that for a flyby to be feasible, in addi- 
tion to satisfying condition (7.1) and (7.2), it is necessary that 
the launch point at the Earth's surface be the perigee of the orbit 
and the velocity at this point be directed horizontally. 


Solution. This condition follows directly from condition 
(7.2), since if the launch point is at the surface of the Earth, 
for a flyby trajectory tangent to the Earth at the launch point 
this point can be nowhere else except the perigee. In Section 3.5 
it was shown that the orbital velocity at the perigee: and apogee 
is perpendicular to the radius-vector, so that the launch velocity 
in this case must be horizontal. 


However, this statement can be proven more rigorously. 
Actually, the equation of a conic section pf tecos@)>p/fi+e cosy). 


corresponds to the flyby condition (7.2). When ¢ is varied from 
O to n7.this conditions means that 6 > 6° The latter is possible 


only when ¢ = 0 (the launch point coincides with the perigee). 


Now let us show that when the launch point coincides with 


a | _Awz,! 
the perigee Ay=E Zs Py = O,H - (Fig. 15), where oa),=(7%. %)>! 
and Po 3-8) is the angle of Vy with the local horizon. 


Let us use Eqs. (4.14) and (4.15) (see problem 4.10) written for 
O<QeRe, 


ce? 
b e Sir Yo= a5 ar - i, ? \ 


u 
e cos, = 


! a 
| Knowing that },=° 9, : 
(see problem 4.7) and 


t 


C=|R,x%, |= Rv, sina | 


Fig. 15 


™N 
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(see problem 3.4), we can write 


R? x? sin?a 
%, ‘ GR; Rs ve WANA, Rs gRs ( 
| 
a 2 .-.2 x. 8 | 
© 005 Y= RV, sin rs -{ = Vo Si Ay | 
GRE Gee 


‘When ¢) = 0, Eqs. (7.3) can be reduced to the following expressions: 


vy sina, ' 
sinPaa ctgdy-D, € = — ep tS 
SLTL 0 g 0 q gF 3 (7.4) 


For both equations (7.4) to be satisfied, ctg a) Must be 
equal to zero, since sin a6 cannot tend to zero by virtue of the 


condition 0 <e <1. Therefore, a = + 1/2, and 8, = 0, 7, and 


Q 
this means that the launch is horizontal. 


Thus, for a flyby to be possible from a launch from the Earth's 
surface simultaneously three conditions must be satisfied: 


v V9 < Vy pans GorOt Ot, pe Oya! (7.5) 


1.ci 


If a # w/2, for any initial velocity the body launched from 
the Earth's surface cannot be inserted into flyby trajectory. In 
this case only a ballistic trajectory can obtain. 


When a body is launched from altitude H above the Earth's 
surface, satisfying condition (7.5) as before ensures a flyby, 
however they are not mandatory for its execution. This means that 
the insertion into orbit can be made at an arbitrary point of the 
orbit, which is neither at the perigee or the apogee. The latter 
fact presupposes violation of the condition that the velocity Vo 
is horizontal. 


Problem 7.2. A satellite is inserted into elliptical orbit 


at altitude H above the Earth's surface with velocity Vo directed 


horizontally (along the tangent to the circle drawn through the 

launch point from the Earth's center). For the satellite to move /73 
around the Earth for an extended time it must not approach its 

surface more closely than the distance h(h < H). What must the 


launch velocity Vo be to ensure this condition? 
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Solution I. From the condition of the problem it follows 
that if the launch velocity is directed horizontally, the launch 
point at altitude H can only be the orbital apogee 


(Mnox® Ty ORg tH), | 


since the satellite can attain a lower altitude (h < H). The 
constraing h < H means that the perigee distance r_ must obey the 
condition " 


Prin > x >; +h, 


Let us write the formula of the launch velocity Vo = Vy 


analogously to (5.16) (see problem 5.19): 


= mes ees 8 


rs apts (1 RgtH 


» ~ 
Ry+H+rry RetH RetH + 


Let us evaluate the expression in the parentheses, beginning with 


| 
the constraints rz > Ry + Ry and Rg tH +Tn 2 2Ra +Hth ‘| 


Roth 5 __ Rg th | 
R, +H+r, ~ eRe +H +h *° 


Therefore, 


2 Opes _ Re+H 
Us PR (1 BR ash} 
and 
ops Rye t+h | 


UVa > VSR p+Hth  RgtH * 


The problem always has a solution, since the radicand is 
always positive. 


solution II. The problem does have a simpler solution. 


Let us write the integral of energy for the apogee and the 
perigee: 


pi_ Hs. yt ~ che, 
“ Ra +H od Re +h 


whence 74 
= _ H-h 
Uy Vg —~ Ap, (R,tH)(R, +R) | 
» =< Ry tH y 
Let us replace ve by its value Rg th OA deriving 


from the integral of areas, after which we get 


- Oy.(H-h)Rath) -4/2 Rath ts, | 
Me YR SAR SAR thy] YOM Reh OR ATR 


This then is the expression of velocity Vo = Vy required for the 


perigee altitude to exactly equal h. In order for this altitude 
not to be less than the given value h, it is necessary that 


_ Resth sd 
Us= Uy 2) “Pep a OR, +H+th °* (7.6) 


Problem 7.3. A rocket is 
launched at an arbitrary angie O65 


to the vertical at the Earth's 
surface. Find the minimum initial 
launch velocity Vo required for 


the maximum separation of the 
rocket from. the Earth's surface 
ry to be equal to the assigned 


distance d (Fig. 16). Determine 
the semi-major axis a and the 
eccentricity e of the corresponding 
elliptical orbit as a function of 
angle Oy > 


Fig. 16 


Solution. As follows from the solution to problem 7.1, the 
trajectory cannot be a flyby trajectory since when the rocket is 
launched from the Earth's surface at point A the flyby is feasible 
only in the particular case when a = + n/2. In this case motion 
occurs in a ballistic trajectory, more accurately, along the are 
AaB, (the are BrA lies within the Earth). 


8] 
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Let us write the integral of energy for the launch point A 


and the apogee a: 
ye? — 2h 2 ye 2h 


Rs ° a 
Simultaneously, from the integral of energies we have /15 
CHIR XO [= [Fx |= Rv, sinae do, 
v, = Ss vy, sina,. 


Substituting the value vy, we have obtained into the integral of 


energy, we get 


uf (1 AL site.) 2m - 2). 


so that the desired launch velocity is 
_,f apg ad(d-Re 
ed NCEE COR (7.7) 


Let us determine the semi-major axis of the ballistic trajectory 


ao from the integral of energy vs Hs Re - +), so that 


sR d*- R? sina 


Qy,-vik, ald -R, sin®d,) ° 7.7") 


We can find the eccentricity from the apogee distance Tyad ~a(iee), | 
from whence , 


Qa = 


a-a _ ol d(d - in® 
e= SS = -f- 2 rt ae (= EaRy sinPag(2d-Ry) | 


ad? Ri sin? a , j 


We ean easily see that from the resulting formulas for a = + m/2, 
the formulas we already obtained in problem 5.19 for tthe flyby 
trajectory follow as a particular case: 


a=t(d+Rs)-4 (hats), 


yee yt «—2bsd MoT Ma | To 


~ R,(a +Rs) (Ta + 7x) a 


a 


Problem 7.4. A material point moving according to the law of /76 


universal gravity is at the position Mg at the initial moment 


(Fig. 16), at the distance ry = Re + H from the Earth's center 


and has the velocity Vo: The angle between the velocity center 
Vo and the line of the local horizon is Bo» and the polar angle 
of the point M, (true anomaly) is ¢5- Determine the eccentricity 


of the elliptical orbit and establish the dependence of angles 
b5 and By: Is it possible from these problem conditions to 


determine whether the trajectory is a flyby or ballistic trajectory? 


Solution. To find the eccentricity, let us use the general 
formula 


pe (7.8) 


where 


° _ WL 2 b 
C27, VU, SIN(T,, Uy) = Ty Vo £05 Po, h=u9- Fe, (7.8') 


so that from the assigned Po» Vo» and Bo we can uniquely determine 
c, h, and é. 
To establish the relation between the angles 6 and By» let 


us use relations (7.3), writing them in our case as 


2 
e sing r2 y? cos f, v2 , _ uecos®s, tgp, | 
0G, Te VERE cost p, 7 G70 


r2 y? cos Po Us COS*By 
€ cos@ = -{= —{ 


3 ’ 
Var 9, "0 
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where g,, is the acceleration due to gravity at altitude H (see 
problem 1.6). Dividing the first equation by the second, we get 


| 
Wye. = ve COS Bo 9 By _ tgp 


Vy COS"By J To 1-9 7, / V9 COS* By ; 


Transforming the denominator of the resulting fraction by means 
of the relations 


2 a 2 
ve cos'p, =~ , Ps 2a 9,17 eee prs (7.9) /T7 
we find 
top, _ t a 
1." 7a ei ay SC! (7.10) 


The orbital parameter p - Mu, can be assumed known since 
from (7.8') we know the constant of areas. Since the eccentricity 
then can be assumed known, then consequently we also know the semi- 
axis a. If angle b4 had not been given in the conditions, Eq. 


(7.10) could serve for its determination from the known By» To» 
a, and e. 


We can: at once from the problem conditions establish whether 
the trajectory is a flyby or ballistic trajectory. To do this, 


it is sufficient to find the perogee distance Pain = °_ 7-4 (7 - e) 


and to compare it with the Earth's radius, by verifying condition 
(7.1). 


Problem 7.5. At which point of the ballistic trajectory 
does the angle of trajectory inclination to the local axis (plane 
perpendicular. to the radius-vector) reach its largest value? 

Solution. Let us use the results of problem 7.4. Let Mo 


(Fig. 16) be an arbitrary point on an elliptical orbit. We will 
denote angle By by 8, bearing in mind the arbitrary point on the 


orbit. Eq. (7.10) is valid for any point on the orbit: 


9 Le 
, ry. esting 
Up = tg y (I - p) ~ F+ecosg’ 


Let us set up derivative ae and equate it equal to zero: 


ae 
dp _ : _ ~e?- ecos@ <0 
do” at.fesn@¢@ \" fi+ecosomy -? 
© T+ Goat) Oreese) 
whence we get cos ¢ = - e fore # 0. Corresponding to this condi- 


Tion is the radius-vector 4 i 
rep/i+ecosy =a(i-e*)/l-e=a.. 


However, in problem 5.11 it was shown that the conditions cos ¢$¢ = 
=-eandr =a are satisfied only for two points on the ellipse 
B and Bt situated at the end points of the semi-minor axis 

(see Fig. 7), that is, when E = «+ 7/2. 


Problem 7.6. At the moment the spacecraft separates from 
the last rocket stage, it is at point M, (Fig. 16) situated at 


altitude H = 200 km from the Earth's surface, and has an initial 


velocity Vy = 8.50 km/sec, where the vector Vy) makes an angle of 


B a = 10°.00' with the line of the local horizon. Calculate the 


constant of areas with the craft trajectory c, constant of energy 
h, eccentricity e, trajectory parameter p, polar angle of launch 
point $), semi-major axis a, and perigee and apogee distances 


rs and Pye Determine from. the data obtained whether the trajec- 


tory is a flyby or ballistic trajectory. 


Solution. Let us use formulas from problem 7.4 for the 
calculations. 


Let us determine the constant of areas by the formula 


C= T% Vz COS Bo = (Ry tH) Vo lOS By = | 
=(6370 + 200)- 8,50: 0,985 =55,000 km“/sec . 


We calculate the constant of energy h by the formula 


. . 398.6 ~ 2 
ft = v2 -5E* 2 (8,50)? - ger = 72,25 121,34 =-49,4/ Ki 


Sec 


We find the eccentricity from (7.8) 


—. oo | 
e “i+ qb = /1-49.4 355500) = ¥1{-0,934 = 0,263. | 
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Parameter p, according to (7.9) is 


_¢c?_ §8000 7_ 
P* ys ~ 398600 = 7590 km. 


Determine the polar angle by Eq. (7.10): 


| tgp, _ 0476 _ 0,476 _ _ eno ng! 
Wy foe 7 _ BBO yas ~ 6,134 = 4,310, @ = 52 38 y. 
? 


where we can verify $6 by using Eq. (4.17): 


pe LE mT 
G9 =a eo 8 
me TB £53000)" — 398,600 


Let us compute. the semi-major axis in terms of parameter p: 


__P  ___ 7590 «7500 | 
a = TleF = 70,060 ~ 70,93) = 8450, km, 


and let us find the distances re and re from a and e: 


ry = 4(4-e) = 8450(1-0,263) = 6000° km, 
T= 2(4 +e) = 8150 (4+ 0,263) = 10,300 pp. 


Since the r_ proves to be smaller than R,, obviously the 


} 
ballistic is.a trajectory one. The flyby is impossible for the 
givenjinitial conditions. 


Problem 7.7. <A rocket at a distance Po from the Earth's 


center is given (Tr, $ R,) some initial velocity Vg directed at 
O 


an angle By to the horizon. What must this velocity be for the 
radius~vector of the perigee to be equal to the given value 


LP) 
re (mr <ry)t 


solution. As follows from the problem condition, the launch 


point M, (Fig. 16) is neither the perigee nor the apogee of the 


N 


orbit. Let us write the integral of areas for the launch point 
My and the perigee fw: 


C=[Fxd| = ToUglOS Py = Fy Uge, 


Here 


v_ -—L2Me00S pe 
XK Ld 


Writing the integral of areas 


Uy - 2 ps /T eu -2 p/n, 


for these same two points, in it we substitute the value of vo: 


2 { { 2 M7 = 
Uy = Ux 2), (4 -) = MEH rere 


_ Je Ue COS’Bo _ 9 Tx 
rz 3 Tove ? 


whence 
a . 

vi [1-(B)eos'p|=-2n, Boz 

and 
/: 2s (%-Tx) 

= $ Q x 

rer (yeast tT 78 | (7.11) 
For the problem to have a solution, the inequality FE cosp, 1, /80 


must be satisfied, whence 


*05 Pe? Fe 5 (7.12) 


so that we can determine from the given By» ro. and Po whether 


the assigned distance re has been attained. 


Eqs. (7.11) and (7.12) enable us to investigate both classes 
of elliptical trajectories for which r 2 Ry (flyby) or res Re 
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(ballistic). Upon inspection of the flyby trajectories we can 
assign not Py and Poy but the launch altitude and the perigee 


altitude H and h_. In this case (7.11) becomes 


{ 


y= 2 Bs (H -he) ; 
o” . RatH \2 2 _{ ; 
(Re+H)(Ry* hs) | (GE ) cos, 4] | (7.13) 
where COS P, >(R, + Ay V/Ry +H), so..that actually we determine what 


velocity must be imparted to the rocket at altitude H for the 
altitude of its orbital perigee to be equal to the given altitude 
hi. 

We can easily see that from (7.13) as a particular example 
Eq. (7.6) follows, obtained for the velocity in horizontal launch. 
Actually, when Bo = 0 we have again 


. -/ 2 pt, (H—-hts) - yn Rsths\ | 
va Me (Re +H) [(Rg+H)*~ (Re ths?) BR, +H thy ao) 


f 


Problem 7.8. Under which initial. conditions will the trajec= 
tory of a spacecraft launched at altitude H from the surface of a 
planet with radius R not intersect its surface? 


Solution. We can consider this problem as a particular prob- 
lem of the preceding one, with r= R (ha = 0) (the tangency is 


included in the case of flyby). By Eq. (7.11) we could easily 
establish the limiting velocity at which intersection will not 
occur, that is-, in which flyby is observed: 


Ee RN | eee! | 
‘lim r_ R] (4a) cos", 4] 7 (R+H)[(R+H) cos" p,-R*] * (7.13') 


so that when V, z Va4m> we get the flyby trajectory (the case /82 


Vv gives the case of tangency). This condition can be 


0°” “lim 
represented as: 


a 
Yq ? Yi.par (R+H)* cosa, R* | “lect ¥(R+H)?cos*B-R*  °, ~.13 


Y¥2 V 


where v OO 
1.ci 


Va ed and v 


1 par * are the 


local eircular and local parabolic velocities at the launch point. 


Now let us look at the conditions imposed on angle By, (with 
local horizon) or angle a, (with the local vertical) (Fig. 16). 
Obviously, the problem of determining velocity has a solution if 

a 2a 2 
the condition (R+H) cos pp-R >9s/ 45 observed, where (R+H)cosp,>R,| 
or _ 


(R+H) sind, >R.. (7.14) 


From the geometrical point of view, this condition means 

that the velocity vector Vo must be directed outside a cone. 
described’ around a planet and 
with its vertex at the launch 


point My (Fig. 17), so that the 


angle of the semi-aperture of the 
cone 6 = 180° - 6 is a function 


: of the launch point altitude above 

the surface of the planet H. 

| Actually, from the right triangle 
MCA we have 


R=(R+H)sin (180°- | 


Fig. 17 -O,,)=(R +H) sind, = 


(7.15) 
= (Rei) 605 Py 


therefore, if the velocity vector Vo, is directed along the cone 


generatrix, the problem has no solution. For any other vector 
Voo directed outside the cone, from the relation of the sides of 
AM, CB it is clear that the leg CB = (R + H) sin Goo = (R + H) 
cos B yo is larger than the radius of planet R, so that the 


condition is satisfied. Condition (7.14) and (7.15) are identical 
to condition (7.12). 
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CHAPTER EIGHT 


| 


TRANSFER FROM ORBIT TO ORBIT 


If at a point with orbital velocity v an addition velocity 
impulse u of arbitrary direction is applied, the point acquires 
the velocity v + u and passes into a new orbit. Here the impulse 
u must be specified in magnitude and direction. We will relate 
the direction of the radial impulse u with the positive direction 
of the unit radius-vector r° and assume that u = ur®. Here we 
will regard the impulse u as positive if u > 0, and negative if 
u< 0. We will associate the tangential impulse u with the posi- 
tive direction of the tangent to the orbit T° (directéd ‘toward 
the side of motion) and we will assume that u = ur® (rule of 
selection of the sign of u-is the same as for the radial impulse). 
In the case of a tangential impulse, the geometrical addition of 
velocities is replaced by the algebraic. 


Problem 8.1. A satellite is moving in a circular Earth orbit 
Po with circular velocity Voie Determine the radial velocity 
impulse u as a result of whose application the satellite passes 
into an elliptical orbit with assigned perigee distance re 


Solution. Let us look at the positive and negative impulses 
(Figs...18 and 19}. In both cases the applied impulse u is perpen- 
dicular to Vo4q therefore the resultant velocities, even for their /8 


——wee 


different directions, will be the same: v = Wve5 + u*. Therefore, 
in both cases the direction of velocity is positive (/ve, + u > 


> Vo? and the satellite passes into an elliptical orbit with 
semi-major axis a > Py: The orientation of the orbital line of 
apsides depends entirely on the direction of u. 

To determine the velocity at the perigee of the new orbit 7, 


let us write the integral of areas for two points on this orbit -- 
--— perigee an and the launch point My? 


™~ 
ee) 
N 


Fig. 18 Fig. 19 
Key: Voa Key: 1. Voar 
V Ve 
par ei 
Vv. Vv 
ci ci 


5 
| 
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~ 
e 
wy, 
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From vector geometry it follows that sin ae = Viq/V> so that we 
have 


O (8.1) 


We can reach the same relation by assuming that the launch point 
belongs simultaneously to both orbits, as a result of which the 
constants of the areas of these orbits are 


€ 
s 9 and c= | r, Uy = rs | Vv ° . 
Cl Cl 


Com xX Une TX D = XV 


4 (8.2) 
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Canceling out H/o by means of the relation 
O 


(8.3) 


we get for u 


so that 


Or 


fetal 2B +(S)] HB ay! | 


Yo rs omen 


Finally, for the desired impulse we have 


UT SY ci a (= -1%). (8.4) 
or | 
u=+tv, (4-1). | (8.41) 


- pa. 2 1: # ° ° . ° , 
where Vi ei a [| is the familiar local circular velocity at, 
the perigee of the new orbit. 


Since the expression in the parentheses is positive (Ty > ri)s 


the sign in front of it corresponds to the sign of impulse u. 


Thus, the given distance rr can be attained for both direc- 
tions of the radial impulse vector. The orientation of the 


ae 
Y 


elliptical orbit here differs. We can easily see that in the 


limiting case for the impulse |u| = Voz (%Qq > 0) parabolas A 
3 - ¥y2. + 2 = = 
and B are formed (v Vv" + u ¥2 Vos Voar with mutually 


opposite direction of focal axes. 


Interestingly, according to Eq. (8.3) the length of the semi- 
major axis of the new elliptical orbit does not depend. on the 
direction of the impulse (for its fixed value). Thus, both new 
elliptical orbits (see Figs. 18 and 19) have identical semi-axes 
and eccentricities (for equal ri), although the orientations 


differ. Correspondingly, the periods of revolution of the 
satellites in these orbits are also the same. (see problem 8.2). 


Problem 8.2. The satellite moves in a near-Earth circular 
orbit with radius Typ making one revolution in the time T. As a 


result of a radial velocity impulse u being applied, it passes 
over into an elliptical orbit. Determine the period of revolution 
of the satellite Tin the elliptical orbit. 


Solution. In problem 8.1 it was shown that the length of the 
semi-major axis of an elliptical orbit’aq formed due to the appli- 
cation of the radial velocity impulse and, therefore, the period 


of revolution Ty > do not depend on the direction of impulse u: 


for a fixed impulse u these quantities prove to be equal both for 
the positive and negative impulses (see Figs. 18 and 19). The value 
of ga can be determined from the integral of energy using the second 
formula of (8.2): 


; { 
vo =v, tut--ES +ut-p,(2 - 4). [85 
ci 
whence 
ao -—hs"s = r 
Peo Teus Wore ue/pe 8 (8.5) 


Assuming a to be known, let us use Kepler's third law: 


2 73 43 
Fate 2 ({- 72 y 
rp wae= (tu) 
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whence 


| ye 
T,<-TU- gee) - (8.6) 


For the satellite of a planet or star whose gravitational para- 
meter pis unknown, by means of the relations of circular motion 
we can cancel out, P/U from Eq. (8.6), then we have 


yt af = (HE)! mo T\? 
ci eo 7 /? ee laar) » 


so that 


net[t- G4t)']. ! (8.7) 


Problem 8.3. The problem conditions are the same as in 
problems 8.1 and 8.2. Determine how variation in impulse u with 
any sign influences the orientation of the focal axis of an 
ellipse relative to the fixed direction at the launch point M, 
that is, establish a relation between u and the true anomaly 
(polar angle) 9 of the launch point. 


Solution. As follows from geometrical relations (Figs. 18 
and 19), variation in impulse u influences the direction of v, 
and therefore, the change in the angle ay = arcctg U/V 4° In 


problem 7.4 we established the formula (7.10) relating angle % 6 


with the angle 6, = {2 Ty (angle with the horizon): 
_ to pe _-_ttga, _ uL ; 
bo = Trp 1- rof/P Ua -PefP) 


As we can readily see, the parameter of the elliptical orbit is 
equal to the parameter of the initial circular orbit, since 
parameter p is defined by the formula p = c?/yu, and the constant 
of areas c, as indicated in. problem 8.1, is the same for both 
orbits, namely: 


. 4 ga * 
Ce MUG MU SIN dy = ro VUi. + wu" SITLL Ag, 
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where 


so that 


pe c? r? vci_ r, 
Hs Hs 


This value of parameter P could have been obtained at once by 
the formula p = a.(1 - e*) for the circular orbit when a = Po: 


and e = 0. The value of the parameter p = ry = a (1 - e*) is 


retained also for the entire family of elliptical trajectories 
corresponding to the set of values |u| < V.4- In this case the 


orbital elements a and e of each ellipse differ. 


Thus, assuming p to be known, we get tg ae > @m, that is, 
b4 in general does not. depend on u, and thus, does not depend on 
v, and here we have a..fixed value dQ = t/f2, or oy = 3n/e2. The 
first value of by (see Fig. 18) corresponds to u > 0, and the 


second (see Fig. 19) -- u < 0. This means that the perigee of 


any elliptical orbit corresponding to a radial impulse with specific 


Sign is situated at a fixed angular distance from Mo: 


of u is changed, the direction Earth center - perigee is reversed. 


The directions of the focal axes of the orbits coincide here for 
u with any sign. 


A variation in the elements of the orbital family can be 
determined by Eqs. (8.5) for a and from p = a (1 - e*) for e: 


oe 
~ f-au® if hy ’ 
1- 7," | (ve) 


{1-2 Vita -/E\e _ 
e= j-2= {-—-= (y= Po. 


Hence it is clear that a and e depend not on the sign of u, since 


ue or |u| appears in the formulas, but only on its value. As u 


is increased, we observe a and e to increase. As a result, there 
is a lowering of the perigee, which is evident from the formula 


When the sign 


~~ 
Co 


Two parabolas formed when |u| = v_. so that v=v.. ¥e2v 
ei ci par 


are the limiting case of elliptical trajectories (except .:cfor the 
case of initial circle). The formulas for o« and e in the particular 


cases. (u = 0 and |u| = Vag? correspond to a circle (a = Py and 
e = 0.) and a parabola (ec > » , e > 1). The perigee distances of 
both parabolas are the same: r=bpetn. Thus, by specifying 


a set of radial impulses lu | < VoG? we observe a gradual evolution 


of the trajectories from the parabola A to parabola B, that is, 
an evolution of parabola A -~- ellipse - circle - ellipse - parabola B. 


Problem 8.4. A satellite 
moves in an elliptical Barth orbit 
I (Fig. 20) with specified radius- 
vectors of perigee re and apogee 


Dee Determine the tangential 


velocity impulse u plight of the 
perigee for which the apogee alti- 
tude is increased by a given 

' amount H. 


solution. The apogee alti- 
Fig. 20 ‘tude can be increased only by 
applying at the perigee a positive 
velocity impulse u_= ur? (u > 0), | 
during which the orbital velocity v = vi. tu (v = v_ + u). Here 
the :perigee of orbit I is the perigee also of the new orbit II, 
and the semi-major axis Op = 1/2 (ro tir.) is increased to Ory = 


= 1/2 (rv, +r + H). This occurs due to an increase in the orbital 
velocity and follows directly from the integral of energy ve = u 
(2/r_ - l/ary)- 


O+ 


From the integrals of energy written for the perigees of both 
orbits we have 


2p Qfly 2p, Hs Ta 88 
Se ed £28 


[Ls 2 hts Aes (T+H) ; 
pov, +uay/ Eo — Qe” fe TytYath ~ VIR(%etPatH) » 


(8.8) 


whence 


PEE (te Ve) (8.9) 


When the apogee altitude is reduced by the altitude H (see 
orbit. III in Fig. 20), we must bear in mind that the negative 
impulse u (u.< 0) is directed opposite to Vi. SO that there is a 


reduction in the perigee velocity from ve to v = vo tu<v_. 


The perigees of orbits I and III coincide, and the value: of the 
semi-axis a, = 1/2 (r_ tr) drops to the value rr = 1/2 (r_ + 


tr - H). We get the formula of the corresponding impulse by 
replacing H in Eq. (8.9) by (-H): 


_ Ks y/ Ta -H -y en | 
= Tx (a re): } (8.9') 


Problem 8.5. A spacecraft is moving in an Earth orbit 
(r = 150°10° km) with circular heliocentric velocity Voq = 29.78 


km/sec. What tangential velocity impulse must be imparted to this 
craft for its new orbit at its aphelion to be tangent to the 


circular orbit of Mars (r5 = 228°10° km). Determine the same for 
the case of light to Venus (ry = 108°10° km) (Fig. 13). 


Solution I. Let us use the formulas of the preceding problem, 


by writing it in general form for u = MG and r = ror: 


ee eo ( are 7 - 7x) =¥2 0. Wart - 75) 


where ¥2 Yor = “par = 42,11 km/sec. For the case of flight to 


the orbit of Mars, 


mor, u(l~>0, H=en-r, tue oa (oa 


and for the case of flight to the orbit of Venus 


rn“, u(fi)<0, H=7-h, uU= =42,0 (fe -+)- 
fe 1) _ joe ) __ 
42,0 (Po - to) = 4eutt( 253 Vi, 2.53 km/sec. 


~ 
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After tangency to the orbit of 'Mars or Venus,. the spacecraft will 
continue to move in an elliptical trajectory, returning to the 
point of impulse application. The transfer to the orbit of Mars 
or Venus is possible only when an additional impulse is applied 
at the point of tangency. 


Solution II. bet us write a series of formulas characterizing 
the transfer of the spacecraft from one elliptical orbit into 
another. The integral of energy is of the form 


a 2 { 
v “Pe (# - a)°0", (2-4). (8.10) 
If it is assumed that the craft, traveling at velocity Vo> 
acquires this velocity by an additional tangential impulse 
u (Vv = Vq + u) being imparted to it, then 
u= U-V,)= Vat a-F-u,, (8.11) 


where a is the semi-major axis of the new orbit. In particular, 
if the initial orbit were circular, Vo - Vay and the corresponding 
impulse would be 


u=v,(p2-2 -4), vVie. (8.11') 


If velocity v is assumed to be equal to the velocity of the craft 
before application of the impulse as the result of which it 
acquires the final velocity v = v+u, then we have 


fin 
| 
= —_ = -_ r 
u Vesn V Vegn Voy V2-Be (8.12) 
where a is the semi-major axis of the initial orbit. For the /90 
final circular orbit with v_..=v_., we have 
fin om - 
- ea | = He | ? 
U Veg (1-¥2 a). Vag r (8.12') 


For Mars, r/a = 3/3.78 = 0.794, u = 29.78 (v2 - 0.794 - [= 
= 2.95 km/sec. For Venus, r/a = 3/2.58 = 0.860, u = 29.78 x (I - 
-~ ¥2 —- 0.860) = - 2.53 km/sec. 


Problem 8.6. A satellite is moving in a circular orbit with 
radius r, making one revolution in the time T. When a tangential 
velocity impulse u is applied, it passes into an elliptical orbit. 
Determine the period of revolution Ty - 


Solution. Let us derive the formula for the period of 
revolution in the elliptical orbit for a tangential impulse in 


and direction. Based on Kepler's third law T°/TS = 7 Jar and 


and the integral of energy, we have 


2 2 
y 2 ; Vv 3 
a =2-¢ and ri = Ty -(2- “3 )~. 
Vi, " V., 
ei r "el 


With the tangential impulse ve =(v.. + uw)? (u with any sign) 
where vod = 2e7mr/T, so that 


ToT [1-22 - ayy 


Calculate the flight time to the orbits of Mars and Venus 
(Fig. 13). In problem 8.5 it was determined that a tangential 
impulse u = 2.95 km/sec > 0, was needed for a flight to the orbit 
of Mars, while an impulse u = —- 2.53 km/sec < 0 was needed for 
a flight to the orbit of Venus. Adopting for the Earth's orbit 


T = 365.3 days = 3.15710! sec, andr = 150°10° km, we get that 
for a flight to the orbit of Mars 
2,95: 3,15-107\ 4 2,95+3,t5 40" \2] 
T, = 569,95 (ecacno ) - z( 3,14 1,5 -10° ) a 


-3/2 1 \3 
= 365,3[1- 0,197 - 0,001] = 365,3 fara) = | 


= 365,3° 4,411 = 515,1 days, t,-47, = 257,6days. | POL 


™ 
-~ 


Under the same conditions, in the case of a flight to the 
orbit of Venus we have 


-3fa 


> 365, 2453 +5445 -107 2453-3115 10)" 
ne 369,3 [1 -( 3.14 + 1,5 - 10° ) 4 So 508 = 


100 


-3/a Vite) 
- 365.3[1+ 0,168- 001) = 365.3/(4te5-) 


‘ = 365,3 -0,798 = 291, 9days,t,-4 T, > 146.0 nays 


The same values of t were obtained in problem 6.10 by another 
method. 


Problem 8.7. A spacecraft moving in a circular satellite 
orbit around the Earth must be launched from it, receiving a 
tangential velocity impulse, and must be inserted into hyperbolic 
orbit with given velocity v.. At which radius r of \the initial 


circular orbit will the impulse be smallest? 


Solution. Insertion into hyperbolic orbit can be effected 
only by applying a positive tangential impulse u (Fig. 21), when 
the circular velocity of the craft rises to the hyperbolic 


velocity v = Vad + u > Voar = Ve Vis From the integral of energy 
2 2 2 


Cc 
for hyperbolic motion ve - ve =v -evi, =v =h it follows 
par ci co 


that (v_. + un)? _ ove, = v, so that u = v2v2, + v2 - v_.. Since 
ci om e0 ci oo el 


dv.,/dr = au /r # 0)/dr, the condition du/dr' = 0 corresponding 
O 


to the desired minimum of. impulse u 
eat 2 ( 2Ucej -1) dvci _ 0,. 


c.r 2Ue + US ar 
can be satisfied only when ve = ove =.y° = vey,/r 
| oma par to 


This corresponds to the impulse 


u= Vow hyperbolic velocity 
Vv = 2v_,, and radius-vector 
om 5 ) 
| r= 2u /vo = 2y./h. 
6” $ 


Problem 8.8. A transfer, 
which is called a transfer in a 
? : Hohmann ellipse, is effected from 
an initial circular orbit J with 


1 
eircular coplanar orbit F with 


radius r.. A A of orbi 
Fig. 21 Ss Yr, tT point of orbit 


radius r. (Fig. 22 a) to the final 
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—_— Ea 


— 


J the spacecraft receives a positive tangential impulse Uys As a 


result, motion continues along the ellipse T (a > v1) to point B 
(orbital apocenter T), where a second positive impulse U, inserts 
the craft into its final orbit F (r.5 > a). Determine the total 


incremental velocity u = u, + Us in transferring from one orbit 


lL 


to the other and determine at which ratios of ry TO rs, will the 


ratio u to the initial velocity Voq.l be at a maximum (the problem 


of transfer via a Hohmann ellipse). 


Solution I. Let us look at the initial craft velocity at 


point A as it moves in the circular orbit Vo4 yevn/r, - The value 


Vol + Uy at this same 


point, which is the pericenter of the transfer orbit T with semi- 
axis a = l/e (ry + r,)s can be determined from the integral of 


areas: 
ap (2 | / 2 e ys 
Vv = = eee = —_—_—s), = = e 
K rae a p a) " M+r, 


of the velocity after the increment v= 


so that 
= _ a jt rn 
Vo Voi. 1 te “FUF,° 
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Thus, for transfer to orbit T, the craft must be given the tangen- 


tial impulse 
— ~~ = pop Pe —_ _ ‘ 
MV 7 Voila Vn Fer TV 78> 


When impulse u. is applied at the apocenter of the transfer 


2 
orbit, that is, at point B, the velocity after the increment 
V =v +u,, where 
ci.2 a O 
Vot.2 ~ rs ° 
2 2). 72h, 2 y 
Va pe - Far) * ro My+le Om m” 


so that 


= yy yf te nf ee. oO. 
Un = Voy 2 -a VU ee 7 Oe 


Let us set up the formula for the overall velocity increment 
U = uy + U5, by representing it as 


-fFe ary + foe | 
ra roi+Ta) Pytls ° 


Introducing the notation R = ro/ry > 1, we can write 


Toa “V+ -Vateem + VOR “te 
- (FE (1-f+fh-1>o. | 


(8,13) 


Let us determine the desired maximum of this quantity by the formula 


to (Uy ld PRY +7 R-{ (PRY hn -0 
R Vad 1 7 a Tok R(1+R)* oR) ORPva™ Ss 


whence there follows the cubic equation R? = 15R¢ -~ QR-1= 0, 94 


which has a real root R ax = ro/ry = 15.58176. When R > Raax? 


W/V As 1 decreases. A plot of this quantity as a function of R 


is shown in Fig. 23. From this figure it follows that the energy 
outlays necessary, for example, to reach the orbit of Pluto from 
the Barth's orbit, are greater than to reach the orbit of Mars, 
but smaller than to reach Jupiter's orbit. The energy outlays 
for heliocentric Earth-Saturn and Earth-Uranus flights are about 
the same as for flights to Jupiter and Neptune. 


Solution II. To determine 


uy and Uy, we can use Eas. (8.11') 


and (8.12'):; 


6 5 0 0 2 26 30 3 Wy 
R= 


ome 


Ts 


Fig. 23 after which again get 


saturn . 
Jupiter vu Cau) Yh _ far (1- t) fig 


Mars - 14R 
Uranus ci.l | vB 


Neptune 
Pluta | 
Vv. This formula is valid for any R, 
cist and when R < 1, u < 0, and when 
R > 1 (as in this case), u > 0. 


Key: 


SIO FW 
<i 


Problem 8.9. Using the formulas in problem 8.8, calculate 
the total velocity impulse u = Us + u., needed for a transfer via 
a Hohmann ellipse from the circular Earth orbit (Vea i= 29.78 
km/sec , r) = 150°10° km) to the circular orbit of Mars (Tr, = 


= 228°10 km). Characterize the velocity change as the craft 
transfers from orbit to the other. /9 


Solution. Let us use Eq. (8.13) for the calculations, where 


7% _ 228 _ | 2 
R FT =“A5t = {, 52. . oO that V . 1 = mee ( -*) t 
| | ci.l_. 


+ + = VERE - 0.88) +¥055 -1 «0,976 + 0,811-1 


=0,187, and uw = 0,187-v ,, = 0,187+29,78 =5,57km/sec. | 
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The distribution of velocity impulses as the craft transfers 
from the Earth's orbit to the orbit of Mars is carried out as 
follows. On being launched from the Earth's orbit (Vos 1° 29.78 


km/sec), the craft given a positive impulse u, = 295 km/sec (see 


problem 8.5) sufficient for the transfer ellipse T at the aphelion 
to touch the orbit of Mars. The elliptical transfer is charac- 
terized by these velocities (see Fig. 22): 


U, = Uy = Ui + U,= 29,78 +2,95 = 32,73 km/sec, | 


Uz = Uy = Uy oa = 32,73 BS = 21,53 km/sec a 


The second positive impulse Uy =u-u, = 5.57-2.95 = 2.62 


km/sec is important to the craft at the aphelion of the transfer 
ellipse B and inserts the craft into the circular Mars orbit with 
circular velocity v,, 5 = Vp t u = 21.53 + 2.62 = 24.15 km/sec. 


This circular heliocentric velocity in the orbit of Mars can be 
obtained by the formula 


Voa.o WE = V Seecoe = 583,20 = 24,15 km/sec., 


Problem 8.10. Determine the total velocity increment in 
transfer via a Hohmann ellipse from initial circular orbit with 
radius ry to a final circular orbit with smaller radius Ly (Fig. 


22 b). 


Solution. For point A, which after application of a negative 
velocity impulse U4 becomes the apocenter of the Hohmann transfer 


ellipse T, we have 


| 
ea 


= _ = kK er Hef, | 
44 Vo Yoi.1 J ™% Mtr ry 


When the impulse is again applied at pericenter B, we have 


=e eu (2 - d).,/eh. 2 | 
Yoyo Pare ee VAG a VE ae, 


™ 
\O 
On 


Thus, we again obtain Eq. (8.13): 


7 < = fee (1-DAei<o, | 


ci.l 


where Re Za £4. , 
4 H 


Problem 8.11. Using the formulas in problem 8.10, determine 


the total velocity impulse u = uy + Us needed to transfer a space- 


craft from the circular Earth orbit (Vos 1° 29.78 km/See r, = 
= 150°10° km) to the circular Venus orbit (r., = 108°10° km). 


Solution. setting 


in Eq. (8-13), we get 


feed, “pee 
U/V 5 4 “Vitveo 1,389) +V7,589 -4 0,176 , | 


and the total impulse u = 0.176 vo = — 0.176°29.78 = - 5.24 


km/sec. tet 
The distribution of impulses occurs as follows. During 

a launch from Earth orbit (v_, 1 29.78 km/sec), the craft is 

given a negative impulse u, £°+2753 km/sec (see Problem 8.5) suf- 


ficient for the transfer etlipse T to touch the orbit of Venus. 


Further, we have v, = v. =v + u 


A oO ci.l Hi 
ve v/v = 27.25°150/108 = 37.77 km/sec. 


= 27.25 km/sec, Va =v = 


The second negative impulse Uy =U- UF - 5.24 + 2.53 = 


= 2.71 km/sec is imparted to the craft at the perihelion B and 
inserts it into the Venus orbit: Voq.o = Yp + u, = 37.¢f - 2.fl 


= 35.06 km/sec. It is precisely this circular velocity along the 


106 


Venus orbit that we find from the formula 
-i/ fe 4427 +10 ! 
“Vr, > Votoa-toe = FeeS =! 35,06 km/sec. 


Problem 8.12. Using the results of problem 818, determine 
the total velocity increment for the limiting case of transfer 
to an “infinitely distant" final circular orbit (r, >> r,) (Fig. 
24 a). 


° t 
Solution. Determine the limit of Eq. (8.18) as R=jt—oo: 
Jing ge» Je, [Ve (1. +p -1} = v2 


Thus, the total velocity increment tends to the limit 


“= uw, + wu = Yo -4 | caer 0,44 SO | 
(ee = ( Von \Yord " 


that is, from the standpoint of energy outlays this transfer is 
equivalent to the departure of the body "to infinity" in a para- 
bolic trajectory whose pericenter " at the point of application 


of the unit velocity impulse u = (/2-1) Voi.1 Voar — Youll’ 


Solving the inverse problem -- on return "from infinity" to the 
initial orbit J, which is equivalent to the problem of the depar- 
ture of the point "So infinity" with final orbit F, we obtain the 
formula for the total negative impulse u = u, + u, = - (/72-1) 


_ 1 C 
Vago 77 0.41 Voi.2 * 0. 


Problem 8.13. Determine the overall velocity increment 
needed to transfer a craft along a three-impulse bielliptical 
trajectory (Fig. 24 b) when its motion originates initially along 
ellipse tT to point B', where it is given a second positive 


impulse u, (see condition of problem 8.8), and then along ellipse 
T, to point C, where this craft receives a third, negative impulse 
U2 inserting it into the final orbit F. 


o> and U2 > let us /98 
use successively Eqs. (8.11'), (8.11), and (8.12'): 


u -PE((Py-), R=Ft. ae, \ 
uel (2-0 VE Vara Seta 


Solution. To determine the impulses Uj» U 


™ 
\O 
\O 


wf Eb kee + 


| since according to the theorem of 
| areas we have 


a (VE +m) Payee 


Finally, 


es 21 a Viet! + 
| + PCO acisay *VR ~VRCGTR * (8.14) 


We can easily see that from Eq. 
( (8.14) when R = Q (r., = Poe) Eq. 


| (8.13) follows as a particular 
; case. 


Problem 8.14. Determine the 
energy needed to rotate the plane 


Fig. 24 


Key: Yéi.l of the circular orbit of a craft 
Vou. moving along it at velocity 
Vis os ee ~~! ° 
ci. V4 Vag y-E. by 90°. 


Solution. Suppose a spacecraft is moving along a circular 
orbit with velocity Vy and at some moment of time receives the 


increment Av, perpendicular to the radius-vector r at point M. 
Then the resultant velocity V5 = Vy + Av, will also be perpendicular 
to r, and if lv. | > Iv, |; the point M of application of the impulse 


will become the pericenter of the new elliptical orbit (Fig. 25). 
Here the plane of the new orbit will be inclined to the plane of 
the initial orbit at an angle i (the angle between the vector is 
vy and V5 ). 


From the vector triangle we have by the theorem of cosines 
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‘ 
| 


¥ v2 +y? a 
AU=¥ 0 +Ul-24,¥, Cost .«: 


To conserve the circular velocity 


| Vi» the craft must be given the 
| impulse AU, =2%, sin Ef 
Fig. 25 In this case the inclination of 


the new orbit to the plane of 
the initial orbit i 1s arbitrary. 
In the limiting case (when i = 90° and Vv, = V5) we obtain /100 


nv, =073 -¥2E, | (8.15) 


Let us perform the calculation. Let the spacecraft move in a 
circular orbit of the Earth with circular heliocentric velocity 


v, > 29.78 kn/sec (artificial planet). For the craft to be 


inserted into a circular heliocentric orbit whose radius, is equal 
to the radius of the Earth's orbit and which is inclined to the 
plane of the Earth's orbit by i = 10° there is. ‘required the impulse 


{ 


Av, = 22, sin =2-29,78-0,087- 5,18; km/sec, 


and for insertion into a polar heliocentric orbit (i = 90°), the 
impulse | 
A= YU V2 =V2 -29,78 = 42,41 km/sec 


is required. Thus, the energy necessary to rotate the plane of 
the circular orbit 90°, m(Av)?/2, is equal to the energy needed 


to provide the circular velocity /Ze at the point of applica- 
tion of the tangential velocity impulse, that is, to transfer 


the craft. from a circular orbit to a parabolic orbit lying in 
the same plane. 


CHAPTER NINE /100 


SPHERE OF ACTION. PROBLEMS OF THIRD ESCAPE VELOCITY 


The sphere of action is an arbitrary geometrical concept 
employed to "delimit" two gravity fields. The term sphere of 
action of a lesser gravitating mass relative to a’ greater refers 
to an imaginary.surface within which it is useful to assume the 
lesser mass mas the central body, and the greater mass M as the 
perturbing body. The relation of the accelerations of a point 
with mass m relative to a point with mass M (m < M) at the boundary 
of the sphere of action is as follows: 


Perturbing acceleration due to M | _ Perturbing acceleration due to m 
Acceleration due to central mass m Acceleration due to central mass 
Mw 


If we neglect the perturbing acceleration, we can conditionally 
assume the sphere of action to be a sphere of equal zero accelera- 
tion delimiting the two gravity fields. The radius of the sphere 


2/s | 
of action must be determined from the formula p=a(m/M)’, 


where A is the distance between the masses. From this formula let 
us find the radius of the sphere of action of the Earth relative 
to the Sun, assuming m = M, = 6: ior! g, M = M, = 1.97°1072 g, and 


s = 149.6°10° xm: ° 
ane? \ ef 
p = 149,6-10°( S10.) "= 929,000= 930,000 km. /101 


Similarly, we can determine, for example, the radius of the Moon's 
sphere of action relative to the Earth (p = 66,000 km) or the 
radius of the sphere of action of the Sun relative to the center 


of our Galaxy (p = g"107* km). 
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In the problems presented below, we will allow for the dimen- 
sions of the sphere of action of the Earth only for geocentric 
motion of a point, For heliocentric motion we can neglect the 
dimensions of the sphere of action and assume that a point at the 
boundary of the Earthts sphere of action has the same heliocentric 
circular velocity as the Earth: 


_f_Pe [Ee _ 
Y1.ci/o Ys V740,6-10° ~ 1 149,6-10 E0678 KM/SEC « 


If the point (rocket), on departing beyond the boundary of the 
Earth's sphere of action, has sufficient residual geocentric 
velocity Vso. roc/#? its resultant heliocentric velocity is deter- 


mined from the théory of velocity additions: 


Yroc/o  ‘l.ci/o + "bo.roc/+* (9.1) 


The positive direction of velocities always is associated 
with the true direction of the Earth's orbital motion, which is 
taken as positive. 


Problem 9.1.. Determine the velocity at which a rocket must 
be launched from the Earth's surface |v,_| for the rocket to stop 


la 
at the boundary of the sphere of action (v = 0). How 


bo.roc/t 
much time will this flight take? Characterize the further helio- 
centric motion of the rocket. 


Solution. Let us write the interval of area of the launch 
point of the Earth's surface and at the boudary of the sphere of 


action: 
RZ - a Oe! _ o -2Hs | 7 
a nr a "bo.roc/+ Pod (942) 
whence when "bo.roc/t = 0 and » = 930,000 km, we determine the 
constant of energy h: R=-=H=- 2-398 00 0,857 km“/sec* < Q, 


p 930000 


such that the trajectory of the geocentric arrival at the boundary 
is an ellipse (stopping at a finite distance is not possible for {102 
parabolic motion). The velocity of the geocentric launch is 


Iva [PR RE = Vi + of =V60,857+ (,107=11,15' km/sec (9.3) 
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and it is close to the parabolic velocity Vrz = 11.19 km/sec. 


The flight time to the boundary should be determined from 
Kepler's equation. Let us find, first, from the integral of area 


; 2 hl 
a: Vi, =v = (7k - l/a), whence 


_ Rk _ 6370+ 348600 _ 
O =p ~R; OF = F3GBGO0 - 83570 (11,15)" 470,200 km. 


The eccentricity of the ellipse e -1- Se = 1- S259, 20,086. From the 


My 
O 


formula | 
0 = a(i-ecs&) = 930,000 km 


we-get cos E = = 0.992, sin E = 0.126, E = 170°44t = 3.015, and 
e sin E = 0.124, so that Kepler's equation becomes 


s/a 3/2 
= ax - so 4,702 +305 
T Tie fF e sln £) 2a x | 


x (3.015 - 0.124) = 1.476'°10° sec = , 
= 410 hours = 17.08 days. 


Let us characterize the further, heliocentric motion of a 
rocket. Within the frame of reference of the above assumptions 
it can be stated that a rocket, losing geocentric velocity, 
"always" remains at the boundary of the Earth's sphere of action 
and will move together with the sphere, that is, together with 
the Earth, around the Sun with circular velocity Voi/o 7 29.78 


km/sec. Actually, this motion will originate along a trajectory 
close to the Earth's orbit, and with a velocity along this trajec- 
tory that will be determined by the distance 149°106 km + 930,000 
km to the oun. 


Problem 9.2. <A rocket is launched from the Earth's surface 


with the second escape velocity Vir = 11.19 km/sec and will move 


along a parabolic trajectory intersecting the Earth's sphere of 
action. Determine with what residual parabolic velocity the rocket 
arrives at the boundary of the Earth's sphere of action and how 
much time will the flight take to reach the boundary? Characterize 
the heliocentric motion of the rocket after departing from the 
Earth's sphere of action. 


Solution. Let us determine the residual parabolic geocentric 
velocity at the boundary of the sphere of action for given initial /103 
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launch conditions: 
= | ef 2bhe y-s98 600 
Ino.roc/t! = Wi par! "Vp" geggor = 0-927 km/sec, 


The flight time to the boundary of the sphere of action can be 
found from the formula for the flight time along the hyperbola 
from problem 6.15: 


t(p)=(3¥Re) V2p-P (pp), : 


from which when p = 929,900 km, p = ar = 2R, = 12,740 km we get 


t = (37398600 )' ¥2-929900 -12740 (929900 - 12740) 


= 0.243*10° sec = 67.60 hours = 2.81 days. 


The character of the further, heliocentric motion of the 
rocket is determined by the direction v_. » which in this 
ci.roc/t 
problem is not fixed. In any case, the theorem of [velocity] 
addition (9.1) enables us to obtain either the increase in the 
heliocentric velocity within the limits 29.78 < Vnoc/o § 29.78 5+ 


+ 0.93 (29.78 < Yroc/e * 30.71) km/sec or.else its decrease within 
the limits 29.78 - 0.93 < Vnoc/e < 29.78 (28.85 < Vnoc/o * 29.78 


km/sec. In the first case the rocket will arrive at the elliptical 
trajectory whose semi-major axis is greater than the axis of the 
Earth's orbit, where the point of departure from the Earth's sphere 
of action. becomes the perihelion of the elliptical orbit, and the 
aphelion is on the other side of the Sun. When Viroc/e = 30./1 


km/sec, the orbital plane of the rocket coincides with the Earth's 
orbital plane. In the second case, the semi-major axis is smaller 
than the radius of the Earth's orbit, the departure point becomes 

the aphelion, and the coplanar orbit is formed when Vnoc/e = 28.85 


km/sec. The direction of motion in both cases coincides with the 
direction of Earth motion (direct orbits). We note that this 
consideration is known to be valid for any residual velocities 

in the interval 0 < ’bo.roc/t! < 0.93. 


The problem of the addition of heliocentric velocities at the 
boundary of the ~ Earth sphere of action and the determination of 


the trajectory of the rockets heliocentric motion is identical to 

the problem of the motion of a rocket moving along an initial 

circular orbit under the influence of the additional velocity 

impulse u or arbitrary direction and magnitude (see Chapter Eight). 

In this case the role of impulse u is played by residual velocity /104 


V . 
bo.roc/t 


Problem. 9.3. With what minimum geocentric velocity must a 
rocket be launched from the Earth's surface for its heliocentric 
trajectory to be tangent at its aphelion to the orbit of Mars? 
Solve this same problem for the case of tangency to the orbit of 
Venus. 


Solution. In the case of the heliocentric trajectory being 
tangent to the orbit of Mars, the minimum launch velocity hen 
is the velocity at which the direction coinciding to the residual 
= s a . C) o ? 
velocity “bo. roc/s coincides.with the direction of the Earth's 
circular velocity so that “bo .roc/& > 


theorem of [velocity] addition (9.1) we get the algebraic sum 
(the orbits are coplanar). Reducing this problem to the problem 
of a positive tangential velocity impulse u =v > QO, let 


O, and on the basis of the 


bo/roc/t 


us use the value u = = 2,95 km/sec (it is shown in 


“po.roc/t 
problem 8.5 that this velocity impulse is necessary and sufficient 
for the elliptical orbit of the rocket at its aphelion to be 
tangent to the orbit of Mars). Writing the integral of energy 
(9.2), we can determine the constant of energy for the geocentric 
trajectory: 


| 
ho = v2 . 28s. (2,95)°- SEE = 8 ,702-0,857 =7,845>0. 
bo.roc/¢ P 928 


Hence it follows that the geocentric trajectory of departure is a 
hyperbola. MIndicentally, this also follows from a comparison of 
the value of the residual velocity at the boundary of the sphere 
of action with the corresponding value of the local parabolic 
velocity (see problem 9.2): 


= 2.95 km/sec, v = 0.93 km/sec. 


“bo .roe/+ 1.par 


Let us determine the required hyperbolic launch velocity: 


=VRe+vp=v7.845+((1,t9)) = 11,53. km/sec. 
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It corresponds to the resultant heliocentric perihelion velocity 
at the boundary v, = Voc/o 7 29.78 + 2,95 = 32.73 km/sec, so 


that the heliocentric trajectory is an ellipse. 


In the case when the heliocentric trajectory is tangent To 
the orbit of Venus, the residual velocity corresponding to the 
minimum launch velocity lv, | must be directed so that there is /105 


— 


an algebraic subtraction of velocities and so that the resuitant 
aphelion velocity Ve = Ynoe/e = 29.78 = 2.53 = 27.25 km/sec, where 


u = - 2.53 km/sec is the negative tangential impulse 


“po.roc/+ 7 


necessary| and: sufficient for the elliptical trajectory to be 
tangent at its aphelion to the orbit of Venus. (see problem 8.5). 


The constant of energy corresponds to the geocentric trajectory 


C 9 ~ £. - . 
"bo. roc/$ _ 2h = 6.401 - 0.857 = 5.544 > 0, so that the 


geocentric trajectory of departure as before is a hyperbola. 
Actually, by comparing velocities it follows that 


h = 


y 
bo.roc/t 


2.53 km/sec, and = 0.93 km/sec. The minimum hyperbolic 


Iv) par | 
launch velocity here 


=~~h+v2 = 75,544 + (1,49) = 1.43)" km/sec, 


so that the heliocentric trajectory is an ellipse. 


Note that in both these cases the elliptical orbits are 
coplanar orbits of Earth, Mars, and Venus and have the same direc- 
tion of motion, coinciding with the direction of orbital motion of 
the Earth, that is, they are direct orbits. 


Problem 9.4. What must the minimum geocentric launch 
velocity of a rocket from the Earth's surface be |v,,| for the 


rocket to acquire the heliocentric velocity v needed to 


roc/o 
leave the solar system along a parabolic trajectory. (Problem 
of third escape velocity.) 


solution. The parabolic heliocentric velocity in the Earth's 


orbit is v, par 7 42,11 km/sec (see problem 5.6). Therefore, for 


the rocket to leave the solar system along a parabola it is 
necessary and sufficient for the direction of the geocentric 


escape velocity v to ensure that the quality v 


bo.roc/¢ roc/o 


= vy par is satisfied. Here the minimum launch velocity corre- 


Sponds to the case of the algebraic sum of velocities based on 
Eq. (9.1), which ensures the direct orientation of motion along 
a parabolic trajectory coplanar to the orbits of the Earth and 
the other planets (Fig. 26). 


Let us find the residual boundary velocity: 


= Vv = 42.11-29.78 = 12,33 km/sec, 


¥y 
bo.roc/t roc/® 


i Opium Jeune Let us calculate the constant /106 
4 of energy h from the integral 
of energy for the hyperbolic 


escape trajectory: 


~ fee _ 2_2°398600 —' 


Vv 
bo.roc/¢ 


—— «men sue 


151.18 km-/sec*! > OQ. 


Upexje ~UVnap/e ns 
To this value h there corresponds 


Fig. 26 the hyperbolic launch velocity 
Key: 1. Earth orbit _ 
” 2. Earth : Via! VA + ur = 015118 +49)? =) 
3. Parabola | 
/ “po.roce _ = 16.62 km/sec. 
5. OV 


roc/o  ‘par/e 
The minimum; velocity attained 
is called the third escape 


velocity v The launch of a rocket from the Earth's surface 


TIT’ 


with this velocity ensures realization of the heliocentric escape 
velocity. 


This problem on the escape of a rocket from the solar system 
along a direct parabolic orbit is identical to the problem of the 
application of a positive tangential impulse u = “po .roe/t = 12.33 
km/sec. This velocity is the upper limit to the velocities 
Vv with which the resultant heliocentric orbit will be 

bo.roc/t 


elliptical. By comparing this result with the result of problem 
9.2, we can conclude that with variation in ‘velocity, having any 


direction, within the limit 0 < bo. roe/t! < 12.33 km/sec, the 


heliocentric orbit will be an ellipse. 
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Based on the formulas in problems 9.1-9.4 we can solve the 
problem of the "falling" of a point for the Sun or the problem {107 
of an insertion along a parabola to the inverse orbit. In the 
problem of “falling"™ toward the Sun it is assumed that. the rocket 
at the boundary of the Earth's sphere of action has a zero helio- 
centric velocity, that is, that it is "at rest" relative to the 
sun. Under the effect of its gravity force the rocket begins to 
move toward the Sun rectilinearly (along a heliocentric radius- 


vector). H = — - 2 ~ 
) ere "bo .roc/# Y4/c 29.78 km/sec, lv, 4! 


~ 
ae 


+ ve = 31.81 km/sec, which corresponds to the recti- 


bo .roe/# II 
linear "falling" trajectory lying in the plane of the Earth's 
orbit. The theoretical time of motion can be determined from 
Kepler's’ third law, by examining rectilinear motion as the limiting 
case of motion along a severely elongated ellipse with semi-major 


axis ‘oa 2 150°10°/2 = 75.10° km, the period of motion along which 
is 130 days, so that the time of flight to the sun is 65 days. 


According to the problem of insertion into an inverse circular 
orbit v - 2,1 oye = ~ 59-56 km/sec, Ivi,| * 60.60 km/sec, 


where motion occurs along the Earth's orbit in the direction 
opposite to the true direction of its motion, such that in half 

a year the rocket will again encounter the Earth. According to 
the problem of escape from the solar system along a hyperbola for 
motion in the direction -opposite the motion of the Earth (inverse 
parabolic orbit), we have v = — 49,11- 


bo.wroc/+ — 


bo.roc/s ~ “1. par ~ Yilei 


29.78 = - 71.89 km/sec, |v = 72.80 km/sec: 


la! 
Problem 9.5. Determine the launch velocity Iv, | required 


for a rocket to be inserted into a circular heliocentric ‘orbit 
whose radius is equal to the radius of the Earth's orbit, and 
whose plane is perpendicular to the plane of the Earth's orbit. 


Solution. The problem has two solutions: flight to the 
"north pole" (i = 90°) and flight to the "south pole" (i = 270°). 
In both cases motion will occur with the same circular heliocentric 
velocity Vot/o = 29.78 km/sec as the motion of the Earth, so that 


for an equilateral vector triangle (Fig. 27), the formula estab- 
lished in problem 8.13 for the velocity impulse obtains: 


\ e of . ' 
- ") ; ut _ 1 
= 2 Uy, sin & =y¥2, 


Moo. roe/+! Iv. caso! = IY. par/e| 


= 42.11 km/sec. 


_ Fig. 27 
Key: 1. "bo.roe/+ 
Yvi/o 
3+ Vag 
4 Earth 
OR IGM, 
OF Roop” PAGE Io 


velocity v 


Fig. 27 shows the case of 


, flight to the "north pole," 


when the vector "vo .roe/t 
makes an angle of 135° with 
the vector of the Earth's 

t/o (for flight to 


the “south pole" this angle 

is 225°), In both cases the 
launch velocity at the Earth's 
surface is 


a, 2 2 


. ~ \ —_ 
Iv, ,! “bo .roc/t t Tr 


= 43.46 km/sec. 


™ 
re 
) 
Oo 
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CHAPTRR TERN 


TWO-BODY PROBLEM. 
THIRD KEPLER'TS LAW GENERALIZED 


In the previous sections we considered the limited two-body 
problem. In solving this problem it was assumed that the mass of 
a space object is sufficiently small so that its attraction of 
the central body does not affect the motion of the latter. However, 
in the case when natural celestial bodies interact, the central 
body under the influence of the other body executes some motion, 
which in turn is reflected in the motion of the former body. As 
a result, both bodies execute Keplerian motions relative to their 
common mass:center: (barycenter) with equal periods of revolution. 
In the problems given below we have examined some questions of 
the motion. of the "attracting" mass m in the case when the 
acceleration of the central attracting mass M under the influence 
of attraction by mass m cannot be neglected. 


Problem 10.1. Two free points with masses m and M (m <«< M) 
are moving under the influence of gravitational forces [mutual 
attraction forces]. Determine the law of motion of mass m rela- 
tive to mass M. 


solution. Let us consider the motion of points m and M in an /109 
absolute inertial coordinate system (Fig. 28). For the radius- 


* _ = - a— ~ 
vectors we” have = Py tf, dm _ 2 hn | a?r 
| | at? ~ qt? dt * 


. As a result of the action of the 

; attractive force, mass M will tend to 
shift relative to the origin of coordi- 
: nates in the direction r, and mass m 

| will tend to move in the direction 

| opposite to it. Here the equations 


| of absolute motion of both masses are;, 
d 


Dividing both parts of the equalities by the mass m, we get 


— 


i: ot a a oe d*hm 2 py | 
at? ~ fm Ss ate PM (10.1) 


The equation of motion of mass m relative to M becomes 


ar a? d*o, = 
“abe = “gph ~ “gph PUM +m) Fe (10.1") 
or 
| 
d?F ji : 
ae +F(M+m)2y<0. (10.1") 


By comparing the resulting equation (10.1") with the equation 
of absolute motion (1.5) for a "nonattracting"™ point from the 
restricted two-body problem, we conclude that in this case relative 
motion will occur. according to the same laws as absolute motion, 
but the gravitational parameter py = f (M +m). In other words, 
the gravitating satellite m will move about the central body M /110 
as if a "nonattracting" satellite m were moving around a central 
body with mass M + m. 


From the foregoing it follows that the formulas according to 
which in. Chaps. 149 we determined the various kinematic and 
dynamic characteristics of motion (velocity, time of motion, energy, 
and so on) jcan be used for determining these same characteristics 
of motion also in the general two-body problem. But in this case 
the gravitational parameter must not be wp = fM, but must have the 
new value p = f (M+mM). 


Problem 10.2. Two homogeneous spheres with radii R, and R 


will begin to move from a state of rest under the influence of 
forces of mutual attraction. Determine with which relative velocity 
Vn these spheres will collide if the initial distance between their 


2 


centers is L, and the mass are Mm. and M5 « 
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solution. Let us examine the case of rectilinear relative 
motion in the general two-body problem. As indicated in problem 
10.1, in this case relative motion can be replaced by the absolute 
motion of one (any) sphere in the field of attraction of the other 


(fixed) sphere with mass m4 + M., « Assuming that the spheres 


attract as material points, we will consider that, for example, 
the point with mass m. drops toward the point with mass ms + Ms » 


| 2 
where this dropping occurs from altitude L to altitude R + Ro> 
The force of attraction of mass m4 + mM. acting on point Me 
a 


is - fm, (m, + m,)r- (we associate the positive direction with 


the radius-vector of point m drawn from my + mM.) ; and the positive 


oO? 
work done by this attractive force in moving the mass m, over this 


section is fm (m, + m., ) (1/R, + R, 7 1/L). The theorem on change 


2 
in the kinetic energy of mass M5 when it is in absolute motion and 
has the initial velocity v, = 0 is of the form 


0 


; 2 
MeV f _ fo 
“Sr = fm, (m,+ m,) (ea, L ), | 


so that 


v.. = JF (1, +,)(gtg--). | (10.2) 


5 1 = Mae Ry = FR 


R, = 0, then again we get the formula for the absolute velocity /11i 
of dropping of the point onto the Earth's surface from altitude H: 


{ { af Ms Hi 
v= Y2PMs (a wr) ~ VAR, +H) Vag 


If in this formula we assume m, << m and 


+? 
@) 


that we derived in problem 1.9. Based on problems 1.8 and 1.9, 
it can be concluded that "receding" spheres must have this same 
velocity V., (as the initial velocity) in order to be separated by 


the given distance L. 
Problem 10.3. Two free points whose masses are m and M move 
under the influence of forces of mutual attraction. Determine the 


motion of the points relative to their common mass center (barycenter) 
C (Fig. 28). 
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Solution. The position of the barycenter C, lying along the 
line connecting the masses m and M, is defined by vector Pas Thus, 


(M + m) Po = Mey + mo. 


Tf we let Pu and ra stand for the vectors determining the 


positions of M and m relative to C, then Py = Pp, + Fs Pm =! 
= Pet ™m ’ so that we can write (+m) Po = p. + | 
+ Py) + mp +7) and 
MF, + mF, =O. (10.3) 


Since in this case the radius-vector of mass m_measured_ rela- 
tive to the radius-vector of mass M is of the form r = Mm = r_-r 


m  M’ 
we can write r,, =r -r, andr.=r+t+r. The successive sub- 
M m 14) M° 
stitution of these expressions into (10.3) gives M (r -r) + mr = 0, 
Mr +m (ry + Dy ) = 0, whence (M + m)r = Mr, (M + m)Fy = —- mr 
and 
F = am Fi =- “Mem r, 
(10.3') 
Thus, the orbits described by masses M and m around the /112 


common mass center C are similar to each other and similar to the 
orbit described by a one mass around the other. 


Let us look at the elliptical orbits of two bodies with masses 
Mandm shown in Fig. 29 a. For specificity, we assume M = em, 
which can correspond, for example, to the case of a binary star. 
Both bodies describe about their barycenter C, as about a focus, 
similar ellipses (with equal eccentricities), while continually 
remaining along a line drawn through the barycenter, on the opposite \ 
sides of it. The mass m describes an ellipse that is twice as large 
as the ellipse described by mass M. The relation u,,/M =a ym 


deriving from (10.3) corresponds to this motion. 


Let us evaluate this same phenomenon from the standpoint of 
an observer situated on a large star M for which it is fixed. 
Taking from Fig. 29 a the distances 'mM for each of the instants 
1, 2, ..-, 6 and plotting them in the corresponding direction, 
let us schematically represent the orbit of the star m relative 
to M (Fig. 29 b). Obviously, the major axis of the new orbit /113 
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must be- equal to the sum of the axes of the orbits of both stars 
in their barycentric motion. (Fig. 29 a), so that the similarity 
of. the three orbits is evident. In the general form, the relation 


&. . -. 
SP 5 Se = corresponds to this construction, 


Let us write out the 
equations of motion relative 
| to the mass center with refer- 
ence to the directions Py and 


re: 
mm. 
b gt 1 d* Fin a d*Fm pmy Te 
s mat fm Me “-s an 
SL | | 
Cy 
§ 
< 
_3 so that 
8 | | 
dF a*¥, . 
| Sf =-f Ms, “ata Fm =) ! 
i | 
Substituting in place of r 
Fig. 29 the expression (10.3') we 
Key: X|. Orbit of M obtained earlier, we get 
Y|. Orbit of m 
Z Orbit of m relative 
to M 7 ° 
d*Fm = £M®> Fm «sok Fst ok \ 
dt? (Mem rR? ate (Mem) re (10.4) 


Each of these equations has the same form as Eq. (10.1'), so that 
study of the relative motion of the points of both cases reduces 
to solving an equation of the form 


7 a (10.5) 


Thus, motion relative to the barycenter along the ellipses depicted 
in Fig. 29 a originates according to the same laws as absolute 
motion, but in this case the gravitational parameters are 


_ _fM°? _ _ fm 
Ha (M+m)y* ? bn ~ (Mem) ° (10.6) 


GE 
OF POOR QUALITY 


‘Here the motion of mass'm relative to mass M obeys the same law, 
but the parameter u = f.(M + m) (see Eq, (10.1")), 


Problem 10.4. What kind of function will obtain between the 
periods qT, of the revolutions of planets around the Sun and the 


semi-major axes a, of their elliptical orbits if the motion of the 


1 
Sun caused by the attraction of the corresponding planet is taken 

into account? Give the mathematical formulation of Kepler's third /114 
law, generalized. Consider the cases of motion of bodies relative 

to each other and relative to the barycenter., 


Solution. It is based on problem 10.Vin which it was established 
that the relative motion of a mass when one allows for mutual 
attraction obeys the same laws as absolute motion, but where the 
gravitational parameter p =f (M+m). Starting from the fact 
that for absolute motion the relation between the periods of 


revolution t, and the semi-major/axes of the planetary orbits Oy 


is of the form (see problem 6.5): 


ai _ fe _ fMe . a 
re gee tore = OR aE eo (10.7) 
it is of the form 


a3 _ f(Ma+mi) a ee Fe 
re awe ORV FC em): (10.8) 


for the relative motion of the i-th planet. For two planets we 
get, respectively, 


ai/T? _ f(M, +m) _ M, +m, 
ae/T,; F(M,+m,) Mag+m, * , (10.9) 
or 
ar Te A+ ma / Me 
as | T: {+m./M, - (10.10) 


The Eqs. (10.8)-(10.10) are given the name of Kepler's third 
law, revised or generalized, in distinction to Eq. (1.11) of 
absolute motion. 


Under barycentric motion in accordance with the equations of 
motion (10.4)-(10.6), for motions of small mass Mm. and large mass 


"Gee 
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Mo we can write, respectively, 


| mn. { &m, (M, +m,)* 
= ——i- = OT Se ; 
Treg on Bem, FM, (10,11) 


[im a2, (Metin, - 
Tn, * 2m Fo = aa [as Me tne - "| /115 


Based on the relations for the similar ellipses iy, [Mye Sg, Le 


(see problem 10.3, Fig. 29 a), it can be concluded that the periods 
of barycentric motions En and T are equal to each other. 

1 © 
Reducing the barycentric motion to the motion of a mass m; relative 
to the Sun M (Fig. 29 b), we can again replace the parameters 
Uy and Wa with the parameter p =f (Mo + ms); and the semi-axes 


i 


Q and a by the sum of the semi-axes a. = a + a 5, as the 
Mo 1 + Mo m4 


result of which we again get the formula 
_ (2s + @,,,) _ ai 
= ORY F(Me me) ~ ONV FTA Hm) (10.12) 


From the relation for similar ellipses 


Em, = Aste - OX ute + Qm; 
M me M + my, 


there follows 


M, tm, 


(ing Bes « (Se) (iH, 4m,)*=(- Se) (+m) 
© | 


The latter equation allows us to state that the period of relative 
motion coincides with the periods of barycentric motions: 


+ Problem 10.5. Determine the circular velocities of the Moon 
Cm) in its circular orbit (r = 384,400 km) and of a rocket 


(m = Q) executing circular motion in the same orbit. 


Solution. To determine the velocity of the rocket, let us 
use the formulas of the limited two-body problem, so that 


ps; 398600 | 
V =V-r- =/380400 = 1,018 km/sec, / 


Proc 


To find the Moonts velocity, we must adopt the formula from /116 
the general two-body problem and use the formula of rocket velo- 
city, replacing the parameter y= FM, in it by the parameter 
Hse 7F(Me +My) =pg the Knowing that the gravitational 


parameter of the Moon a 4900 km>/sec°, we have 
_ f/f PatPc _ ,/ 398600 +4900 | 
Ug = Pr 334400 = 1,029 km/sec ; 


which means that the Moon will move faster than the rocket. 
Actually, using Kepler's third law we can determine that the period 
of revolution of the rocket is greater than the period of revo- 
lution of the Moon: 


_ prs _ fr> oq yf r* 
Toe 7 eV pg = 2™VFH,* = 2"V ng AVE EMD ° 


By comparing the results we get 


T > 
POC 2 fan =Vi+ Ge = 1+ 316 ahy (10.13 
ro é 


« 


>T,-| 


therefore, TL, , 7 
Problem 10.6. The rocket will move in the circular Moon 

orbit, where at the initial moment it is at a point diameterally 

opposite the position of the Moon. How will the mutual disposi- } 


tion of the Moon and rocket change? 


Solution. As was shown in problem 10.5, the Moon will move 
along its orbit faster than the rocket and will overtake it. After 
a certain time t has ellapsed, the Moon will catch up to the rocket, 
that is, the rocket will fall behind it. Let us determine the time 
t by means of the formulas from problem 10.5, assuming that the 
circular motions of the Moon and the rocket are uniform. We can 
write the angular velocities as follows: 


ow 


7 1,025 -1,018 _ 0,007 -2 4 
— 2 = ¢ - 
We -@ sel“ = Saaa00 7584400 E29 segT 


_ 
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Obviously, the Moon will catch up to the rocket when its angular 
distance from the rocket, initially oy = Ts will reach $ = é7, 
so that 


a q 
q = eo fo NO 05 -10® 


Q- Win O-O5. 182 sec * 6 years, 


We can solve this problem by using Keplerts third law, /1LT 
generalized, by expressing the angular velocity in terms of the 
periods of revolution: 


Let us transform the expression within the parentheses by using 
Eq. (10.8): 


Te. Mey 4. (4-4. Me) 24 Me 
I-pisi-(tegt) =4-( at) =e ME 


Calculating the mean period of revolution of the Moon 


_ r3 _ (384 400) _ 
I =2h Wore e, = 239860074000 * 27,32 days, 


i 


we get 


3 
T 2 hae = 81,5 I= 84,5°27,32 = 2,235+10 . days, = 6 years. 


CHAPTER ELEVEN /117 


MISCELLANEOUS PROBLEMS 


Problem 11.1. Determine the useful work that must be done 
by @ rocket engine to lift a spacecraft with mass m to altitude 
H above the surface of a planet and at this altitude to impart 
circular and parabolic velocities to it. The weight of the space- 
craft at the planetary surface G; planetary radius R, and atmo- 
spheric drag can be neglected. 


Solution. The work done by the engine A is made up of the 
work expended in lifting the craft vertically to altitude H and 
the work needed for the initial-velocity vector to occupy a specific 
position in the plane of the local horizon. Thus, the total incre- 
ment in the kinetic energy of the craft mv¢/2 (we assume the initial 
velocity to be zero) is equal to the sum of the work done by the 
planetary gravity and the useful work done by the engine: 


mov 


3 =pm(gy ~ —a)tA ~~ atgeay the 


where »p is the gravitational parameter of the planet. Obviously, /118 
for the craft to be inserted into an orbit with circular velocity 


Vad Var | » the following work is necessary 
_ pom mH m(R+2H) 
Ag 7 WORRY * ROGHY* ERS AY + 


and for insertion with parabolic velocity v = /7H |, the 
par Rti 


following work must be done: 


A = Reet Rte RO 


par 


which does not depend on ascent to altitude H,. 
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If the weight of the craft on the planet G = mg, where g is 
the acceleration due to gravity at its surface, using the formula 


u = PM = gR- (see problem 1.6), we get 


_ GR(R+2H) _ 
oi 77 BR+A)” | and ‘var = GR. 


We can easily calculate that for a craft weighing 5 tons lifted 
directly from the Earth's surface with parabolic velocity, the 


desired work Adar is 31.85'109 kem. 


Problem 11.2, A spacecraft with mass m approaches a4 planet 
along a line extending through its center (along the radius-vector). 
At which altitude from the surface must the engine be fired so 
that the constant braking force mT it produces will ensure a soft 
landing (landing with zero velocity)? The velocity of the craft 
at the instant of engine firing. is Vo» the gravitational para- 


meter of the planet is u, and its radius is R. The attraction of 
other celestial bodies, atmospheric drag, and change in engine mass 
can be neglected. 


Solution. This motion occurs under the influence of the 
planetary gravitational force and the braking force of the engine, 
whose resultant of which is of the form 


The change in the kinetic energy of the craft is equal to the sum 
of these forces expended in moving it from the beginning of the 
liftoff (instant of engine firing) when ry = R + H to the soft 
landing (r = R) 


mV ; 
lan _ Lue = pm (t - +) + mT(R-10), 


By substituting v = 0 in this equality and dividing both sides 


lan 
by m, we get the equation 


ve yt (h)-R) R-r.). 
“2. RT, + TC 0) fits 


from which we can determine the altitude at which the engine is 


fired, by solving the corresponding quadratic equation for Po. 
The altitude is 


2 \2 
Haz! | +TR +t + e +TR+ a) yr bR. (11,1) 


Determine with which sign we must take the root. By denoting 
the radicand with A, let us write (11.1) as Bava - 


~aTH+TR-+, 


» and let us write the integral ‘of energy in the 


2 as 7 
form B=TH-h +H, 50 that Tr, - Ft YAS From 
this it follows that when T>o@tm|, the root must be taken with 
the "plus" sign, and when I< SA -- with the "minus" sign. 
| Me = [T- gefgy] H > Oe: ine it 4 
But since a2 R(IR+A) ! for a soft landing it is necessary 
} Be | 
that ! > RRFAY >? RIDE This means that in Eq, (11.1) 
. . J 


there must only be the "plus" sign in front of the root, 


Problem 11.3. Set up the equations of rectilinear motion of 
a point with mass m in the Earth's gravity field for the cases of 
ascent from the Earth's surface (r = R,) to altitude H and falling 
O 

to the Earth's surface from the same altitude without initial 
velocity. Determine the time of motion and calculate it for the 
cease H=R,. 

+ { 

O 


solution. The equations of rectilinear motion in the Earth's 
gravity field are of the form 


ae. fmmMs dd? #M Its Ri 
m ata eo e , We 77 ae 2 -£§ 2-2 || 
RA | 

In problem 1.8 the launch velocity Uy = eae needed 

to bring a point to altitude H was determined. It is equal to 


the landing velocity of the point when falling from the same 
height, so that the initial equations are of the form: 


° 2ghsH | 
for the ascent (t=0. Rss hea eV 7% 


for the descent ¢-0, ry=Rs +H, Fa. 
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gh 
in solving the second-order nonlinear equation Fa- Tp we /120 


{ we 


must use the substitution Fdr=Ffdr,' after which the equation 


takes on the form dt a- ZR ar, Integrating it, we get the 

integral of energy P= PRC, | where we find the constant 
c= (G+ BE) from the initial conditions, so that in both 

cases C = - £85. Let us determine the time of motion from the 


integral of area, by transforming it into a first-order equation 
with separable variables: 


r R+H r 
at=*Vargaacry  ** RV ag Vera 


where the "plus" sign corresponds to the ascent (dr > QO), and the 
"minus" sign corresponds to the descent (dr < 0) of the point. Thus, tihe 
ascent time is equal to the descent time: 


R+H 


r Ret 
t=— 2g . where I= aa ar. | 


~ R 


The integral I is taken by means of the substitution 


Ten Ns Bee Var Ar ~2sin¥gdg =(l-eos2g)dg, 


so That 
a % 
_I_ .\ ({-cos2g@)dg = 3-4, + SiM 4, COSGQq, 
aig = | (1-cos2g)de = 3-4 
4. 
where 


R ' 
singe Vase = 8 (F- Ge). sing, cosy=Re | 


= 3 ——uee SS 4 - ; 
’, arccos|/p E z 3 GFCcos po carcsin Yr . 


rls 


Finally, we have 


_ 4 /R+H R+H R-H 
te SH ( VR + S54 arccos AM) , (11.2) 


or ; 


t 5) | YRA+ (R+H) arcsln Re | 
| 


(11.2') 


We can easily see that Eq. (11.2t), as R + », is transformed into /121 
the formula of uniformly accelerated motion t <V2H/g,| corresponding 
to Gallileo's formula v=V2gH (see problem 1.8). 


Let us calculate the time of motion of a point for the alti- 
tude H = Ry by formula (11.2): 
O 


t =f Gt (i+ 3) = 2,57/ 827010 = 2065 sec = 34.5 min. 


Problem 11.4. Two points with masses m and M begin to move 
from a rest state under the influence of forces of mutual attrac- 
tion. Determine the time T by which the points will collide if 
the initial distance between them is L. 


Solution. As indicated in problem 10.1, the relative motion 
of the points in the common two-body problem can be replaced by 
the absolute motion of one of them, executed under the influence 
of the gravity force of a new mass M +m, and thus we can reduce 
this problem to the problem of the motion of a moving point relative 
to a fixed point. The latter problem was already solved by us for 
the limited two-body problem _(see problem 11.3). By substituting 
in Eq. (11.2), H = L, and gpR°o = »p = f (M +m), and by passing to 
the limit as R + », we find 


oh 
as 


(11.3) 


We note that this result can be obtained directly from the integral 


of energy -a_ SP oh, Actually, for the case of descent 
(7,=0, Pa= L) | when t = 0, we get AoE and ar. Yep - LT, 
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whence 


' L 4 (EE 
2 Tk gprs | 
4 0 


Making the substitution r=isin?a, | we again get 


x/a 
r=~5 J sin'g dy = Se. ! 


The hypothetical time of fall of the point to Earth adopted 
as the attracting center with gravitational parameter Wy from 


O 
altitude H = L = R,, calculated by this formulas, is t = 890 sec = 


+? 
oO 


= 14.8 min. 


Problem 11.5. <A heavy sphere moves along an imaginary recti- /122 
linear channel extending through the Earth's center. The force of 
attraction within the Earth is proportional to thé, distance to its 
center and is oriénted toward this center. Determine the velocity 
at which the sphere transits the Earth's center and the velocity 
at which departs from the surface for the cases Vy = O and I%!>0, 


and also the time of motion in both cases. 


Solution. Writing the condition of equality of attractive 
force and gravity at the surface of the Earth, we get 


Fe-kimr, |Fl, =k'mRy=mg, and k* = g/R,, on the basis 
O 
of which.the exact law of action of the force will be of the form 
-_ mor 6 , _ . 2 
Pao m ae Fe a TEE F., The elementary work done by this force is 


SA =f .dF <= _ mg rar 
Rs ’ 


moves within the globe is A=- 32 xfhar = - 2g (7 fi 7 "reg de 
beg 


» While the total work done as the sphere 


Determine the change in the kinetic energy of the point as it 
moves toward the center and from the center toward the surface: 


mu: mu? max ay 4 


a) 
me moe _ m7 { 
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where v5 is the velocity of departure at the surface, always equal 
tO Vo: The velocity of transiting the center is 4%eVuj+gk, =VurtagR, . 


For the case of motion without initial velocity, it is x =V9R, @! 


= 7.90 km/sec. The equality Yo = Vi = O means that the sphere 


at once will fall back and its further motion will be an undamped 


harmonic oscillation. If Vo = Vy ? Q, this motion will occur in 


the "channel" whose length would be apparently increased by 2H, 
where H is the altitude at which the sphere is lifted above the 
Farth. The law of oscillation_in this case will be more complex, 
Since the force of attraction F*-p,F¢/r®* operates outside the 


Farth (see problem 11.3). 


The equation of motion of the sphere within the Earth is a 
homogeneous second-order linear equation (the equation of free /123 
harmonic oscillation): 2K? » or H4k?y agi, where 

| mre-K mr . 


keVg/R,. Its general solution, as we know, is of the form #F-C,coskt+ | 


+C, stn kt ! so that r=-€,k sinkt + C, k coskt. 


Let’ us look at four cases of motion. 


I. Falling to the center without initial velocity (v.20, =O): 


Ty=R,, r,=u,=0, C.=R,, C,=0, 7, > R, cos/F ¢ , 
| as TM 


: | _—— 
ty Vp arceos0 = 2y/H | When t-9 rR, Cet,, r 20. | 
It is of interest to compare the calculated ty = 21.1 min 
with the hypothetical time of fall of the point to Earth, 
attracted to the fixed center. It i8, as was shown in problem 
11.4, 14.8 min. This difference is due to the fact that for the 
law of attraction F ~ 1/r¢ the acceleration of a point with this 
same section of the trajectory will be greater than for the law 
of attraction F ~ r. 


II. Falling to the Earth's center with initial velocity 


Vv > QO: . 
| 5 | r=Re r= d, =~/v,J)<0, C, =; , C, 0 He = -luo|/Be <0, 


Tks cosy-f- | v9| fet sin ff t = cosy/Z- t (Rs ~ lool y/H ty PE t), 
ty Bt aretg YOR, t, (v,=0) = e, | When | ‘ =D yay, 


| 
| 
€=t,, 7, 20. 
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TII. Ascent to Earth's surface, where the sphere halts 


(v, = 0): 
r,=0, T= U, = gR, >0O, C,=0, C, = 2 aVgR, y - ; 
. x 
rT, =, siny-f t, ts - 58 arcsin | ak ot,. 
y 
When t-0 rr, =0, t=tg , y=Rz- /124 


IV. Ascent to the Earth's surface, where at the residual 
velocity |v,| > 0: 


To°0y Fam Ve -VgRa rh, C=O, Cy ag, +d? YR 
T= VOR, +02 ye sin ft » by -/(% arcsin VIR: = 


=/ 6 : 
g arcty Lote, te (%=0)- FY at, -4,. 


When | t=0 rp =O, t=-t,, ry =Rs. 


Problem 11.6. A point with mass m moves under the influence 


of the central force Fa-m(4, + 3,),° where wu > Oandv are certain 


constants. Determine the trajectory of the point, considering 


that the incremental force - my/r> can be either a force of 
attraction or a force of repulsion (depending on the sign). 

Set up the equation of conservation of energy of the point, using 
Binet's formulas. 


Solution. Based on Binet's second formula (see Chapter Four) 
in which u = 1/r is the Binet variable and c = r*$ is the constant 
of areas, we can write 


a7u _ F \ { je Ju 
gat tu =-aorge= Gat ps) ig tar to | 
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from whence we obtain the inhomogeneous second-order linear equa- 
tion in the Binet variable 


vu = + 
agi + (1- ¢ u= ars, (11.4) 


from which when v = 0 there follows Eq. (4.11), given in problem 
4.10. The solution of this equation can be represented in a dif- 
ferent form, depending on the sign of the quantity 1 - v/c?. {125 


) : 
i. For the case when 1-73 >0 (¥<e%) | we seek the solution 
of the form ' 


e7(4 - = 


Pa 2 | A (c*-¥) | 
usa 1- <3 (g-€) - fy - [+4 x ; 


| 


x COS /1- %(¢-6)| ‘ 


from whence it is clear that the equation of the trajectory 
Fae is the conic section 


r= . 
~ “Pecos K(@e-€) °? (11.5) 


that is moving, rotating about a focus, where P=—a—~, Keopl-a, 


and e = A) (e and e are arbitrary constants). Obviously, 
| 
when v = 0, Eq. (11.5) is converted into the equation of the 


Pixed conic section r= pS secey? obtained in problem 4~ 10 


=~ 


as the solution to Eq. (4.11) only for the attractive ‘force fu dey 


7 or om 


that is, when motion is strictly periodic. 


In this case, when the incremental force F * l/r? is present, 
the motion of a point described by Eq. (11.5) loses its periodicity, 
that is, when ¢ is replaced with 21 and the radius-vector has an 
initial direction, its value differs from its initial value. Its 


initial value is restored when the angle ¢ is changed by the 
amount oR -an(i-2)" that is, when the radius-vector is 
rotated by an angle that is somewhat larger than ea in 
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the case v > Q (the incremental force is an attractive force), and 
by an: angle somewhat smaller than 27 in the case v < 0 (the 
incremental force is a force of repulsion). When e < 1, the 
motion of the point in both cases is along an ellipse con- 
tinuously rotating in its plane in the direction of motion of the 
point v > 0 or in the opposite direction, when v < 0 (positive or 
negative precession Of orbit, where the angle of rotation of the 


Tocal axis of the precessing ellipse during the time of one com- 


. ‘ ~1/ 
plete revolution is a = 2%[(I - wey "-4], IPF v/e* is sufficiently /126 


small that its square can be neglected, we can write the approxi- 
mate formula thusly: 


az an[(i+$-%)-i]- BS. (11.6) 


The trajectory of motion 
has the form of a rosette (Fig. 30) 
continuously filling with lobes the 
circular part of the plane bounded 
by the circles 


Pain 2 E(1-e)= Zee 


and 


riasa(irejere, 
Pig. 30 
where p, a, and e are characteristics of the ellipse. 
TI. In the case when 1-2-0 ( Yec*). Eq. (11.4) 
becomes 
o = £; = const . | 


f 


Its ‘double integration in » allows us to establish that ‘the 
trajectory is a fairly complicated spirallike curve (precessing' 
spiral) originating from the circle r = Py (d¢ = 0), 


III. In the case when !-2e<0 (%>c*) | the trajectory is a 


spirallike curve with an infinite number of revolutions, originating 
from the circumscribed circle r = tax (r > 0, 6 > ~) 

IV. If p = 0, that is, if the force of attraction F ~ L/r* 
disappears, Eq. (11.4) becomes uniform, and one of its solutions 
(when v > c%) can be the equation of the logarithmic spiral 


per) er? , The solution of the equation of this type for the 
3 


force F ~- 1/r~ was obtained by us in problems 3.9 and 4.2. 


The force fields of these types are encountered, for example, 
in the theory of motion of microparticles ‘however the case of s3 f/127 
elliptical motion has an interesting interpretation also for 
planetary motions. To show this, let us set up the equation of 
conservation of energy of a point moving under the effect of the 
force indicated in the conditions of the problem, utilizing Binet's 
formulas. 


Since the potential corresponding to the given force is of 
the form U =m (yr7! + 1/2 +vr-72), the law of conservation of 


energy T-U = h becomes $mv?-mprt -4£ mor en, , 


whence we have 


; | 
yea 2h +2mu + Vu, 


. @ 
Comparing this expression with Binet's first formula ven e*[( $e) | 


+ u®] we get the equation of energy 
du\* a » 2 
(s5) rui(i--3)-ugh- Sh, | (11.7) 


This equation, defining the trajectory of a particle moving according 
to Newtonian mechanics under the combined influence of the forces 


Fe 1/r* and F ~ 1/r has an analogy in relativistic mechanics. 
It coincides in form with the equation of motion of a particle 
moving according to relativistic mechanics only under the influence 


of the force of attraction F . re. Here we have in mind that the 
mass Of a particle will vary in accordance with the Lorenz formula 


mem “(1 3 ve a 


» where m, is the rest mass, and Vo is the velocity 


0 
of light. 
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By writing out the equation of conservation of energy, instead 
of the kinetic energy we must introduce the intrinsic energy of the 


particle e -mv/playing the same role in relativistic mechanics as 


the kinetic energy in Newtonian mechanics. 


Multiplying Binet's first formula by m°, we get 


Let us set up for the numerator the formula derived from the 
Lorenz formula after identity transformations: 


24,3 3 4.2 yey-4 273 
mv" =m G- 4, - 22 Ue 2 2 
oe 3) (mév pet mv) x 


The quantity move can be expressed in terms of the total 


energy of a particle: E «= @-mpr'=«mv -mpu, |! 
so that mtu? « vit (E +mpu) and -miyt=v2 (£+mpu)*-miv. 


By substituting this expression into Binet's formula, we get the 
equation of energy 


fdu\®, 2 2fp _ £%* mus 
Gig) + @' (- dor) - & salar - Saat (11.8) 
coinciding in form with Eq. (11.6) if m = m, = const. Here the 


0 
role of coefficient v is played by uo /ve, The coefficient u can 


be replaced by Eu/mavs and so on. 


This analogy of formulas enables us to give a qualitative and 
quantitative estimation of the relativistic effect of displacement 
of Mercury's perihelion. 


The change in the mass of Mercury for motion along an orbit 
proves to be more significant than for other planets, since on the 
one hand its orbital velocity (mean) is higher than for these 
pianets. On the other hand, its orbit has a fairly high eccen- 


tricity e = 0.206 so that the velocity change in the motion from 
perihelion to aphelion also proves to be greater than for other 
planets. Because of this, finding the influence of the relativistic 
dependence of mass on velocity on the nature of the motion of 
Mercury, manifested in the precession of the orbit and the displace- 
ment of the perihelion, is easier than for other planets. 


The analogy of equations found above enables us to use Eq. 
(11.6) in calculating the angle a of displacement of the perihelion 
under the influence of the velocity dependence of mass. After the 
coefficient v has been replaced by ue/wes Eq. (11.6) becomes 


oh =npe/e*r ve, 


We know that Vo = 3°10° km/sec. We can determine the constant of 


areas c for Mercury's orbit from these velocities and distances: 


r= 46*10° km, V_ = 57.8 km/sec, r= 70°10° km, Vo = 38.0 km/sec, 


so that ¢c = rpove= Privy = 2.66°102 km°/sec°. The displacement of 


the perihelion in one orbit of Mercury around the Sun (T = 88 days) 


3,44 (1327 -10°)3 - - 
is “ =[37495)* (2,66-107)* 8,62:10" rey - which corresponds to the 
secular displacement a = 7'%e2. This value of a is about one- f129 


sixth of the observed displacement 42"9 per century, whose value 
is due also to other relativistic effects. 


This effect of the displacement of the planetary perihelion 
can be justified also from the standpoint of the law of conserva 
tion of the kinetic energy of a point in a central force field 
(law of areas). Actually, since the planetary velocity of the 
perihelion Vv is at a maximum, the planetary mass m, is also 


larger than the mass at the aphelion mo: Due to this law of areas, 
mv r= mvir. is not satisfied for a point with constant mass, and 
excess kinetic moment appears at the perihelion: 


AK = M, Vv, one m,Us re = UT (m.,,~- M,) = v, Py (m,-™,)= | ! 


139 


140 


However, the orbital planetary motion is impossible without 
satisfying the law of areas. Therefore this excess moment at the 
perihelion AK is compensated by the incremental kinetic moment AK'! 
arising due to the rotation of the orbit around the sun in the ' 
direction of planetary motion, that is, it is caused by the posi- 
tive precession of the orbit. This incremental moment AK' will 
be larger at the aphelion than at the perihelion because roo Pr. 
It is 


16 , | ; 
BK = My Ue MU = (Mar2 -m,7r?)w » myw(ri-ri), | 


where w is the angular velocity of orbital revolution satisfying 


the condition -%a de! By equating AK to AK', we get 
- _. je . 
_Ugrx  va-vi _ 57,8°46+108 _(57.8)°- (38,0)" ! 


WO -Dor * ne r2 ~~ 2(3-407)2  (r0-1D")* —(46 108)? ~ 


- 1.02°1072" seen, 


[ 


Since a century contains t = 3.15°107 sec, the rotation of 
Mercury's orbit in a century in are seconds (1 rad = 2.06°10°2") /130 


is a= w"t = (1,02-10°")(2,06 -10°)(3,15 10%) 3 742, 


Problem 11.7. <A spacecraft with mass m moves in a central 
attractive field under the influence of the force fo ~ HEE pe! 


in a circular orbit with radius Dy: On the craft acts the engine 
thrust yo SEF? (a continuously-acting variable radial thrust). 


Determine the nature of the change in the craft's orbit as a func- 
tion of the change in the magnitude and sign of the constant 
coefficient a. (Problem of flight with radial thrust.) 


Solution. After firing the engine, the craft is continuously 
under the influence of two forces directed along the same straight 
line (along the radius-vector), where when a < 0 both forces are 
directed toward the force center, and when a > O -- to opposite 
sides. The resultant of these forces is 


The replacement of (1 - a)u by the new gravitational constant 
u' Signifies "replacing" the true gravitational u-field by the 
imaginary u'-field. Also valid for the new field are Kepler's 
laws, since no fundamental changes occur in the equations of motion. 
This replacement can be interpreted as a "change" in the mass of the 
attracting center by the amount aM so that the relation of the 
gravitational parameters is as follows: : 

w=fM, ple(i-a)p = ({-a)/M =F (M-aM)= -AH", | 


Let us look at several particular cases. 


I. a < 0 (radial thrust is directed toward the attracting 
center). At the instant of engine firing the u-field is replaced 
by the u'-field, where nu! = (1 - a)u > up. The motion of the craft 
along the former circular orbit with circular velocity vay~p/r, | 


now proves to be impossible since the theoretical value of the 
circular velocity for a given Po in the y'-field is v'=Vp/r, >0,| 


In other words, the actual velocity v < v' proves to be elliptical 
with respect to the u'-field and the craft passes into an ellip- 
tical orbit situated within the initial circular orbit (Fig. 31), 
where the engine firing point Mo proves to be the orbital apocenter. 


II. a = O (thrust is absent)./131 
In this case no changes in craft 
motion occur. 


III. O<a< 0.5 (the 
radial thrust is directed from 
the center). At the moment the 
u-field is replaced by the u-field, 
when un! = (1 - a)u < uy, the actual 
circular velocity y.ypjr| proves 


to be larger than the theoretical 
circular velocity vaVp7e | and is 


elliptical (v > v') with respect 

to the yw'-field. The craft passes 
into an elliptical orbit lying 
outside the circular orbit, and 

Fig. 31 the point M, becomes its pericenter. 


IV. a= 0.5. In this case 
the actual velocity ve¥e/. is parabolic with respect to the : 


O5<a«<? 


u't-field, since the theoretical value of the circular velocity 
ul = JE = / eh = oa » whence 2-v¥2 , and the craft leaves the 
u'-field along a parabola. 

V. a= 1 (one field "damps" another). The gravitational 


field is entirely eliminated and the craft moves rectilinearly 
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and uniformly with velocity v > 0 in a gravity-free space, along 
a tangent to the initial circular orbit. The circular velocity 


with respect to the p'-field here is absent: i er) 
0 


When 0.5 < a < I, we get hyperbolic orbits with respect to the 
u'-field (v > viv2). 


VI. a> I. The effect of superimposing the fields is such 
as to give the effect of the central field repelling the craft. 
The circular velocity v' in the ur'-field does not exist at all 
(it is imaginary). The second branch of a hyperbola not containing 
the focus 0 at which the central mass is situated is the trajectory 
of motion. 


We note that the equivalent of forces F’ = F + T, just like 
the force F, is a potential force, as the result of which the 
actual radial thrust T is a potential force: 


~~ 
= 
W 
NO 


_pm ~f{-appm _ jt re Oh pt rr y. 2 BTR 


POR Up Sm Be r poe (Ue ae 
| 


The thrust of a plane mirror solar sail positioned perpen- \ 
dicular to the solar rays can serve as an example of a radial 
thrust varying under the law T . 1/r*, since the force of solar 
radiation pressure varying according to the law F . 1/r? can be 
viewed as a thrust force (we have in mind heliocentric motion). 


Problem 11.8. Prove that the motion of a spacecraft can 
follow a logarithmic spiral r = rye*$ if, in addition to':the 


attractive force of the central body F v 1/r“%, the craft is 
continuously acted on by a variable thrust R vv 1/r* tangent to 
the trajectory. Determine the velocity along the spiral. Give 
the energy characteristic of the thrust R. Find the time of 
motion along the spiral. (Problem with tangential thrust.) 


solution. Since the point moves in a logarithmic spiral 


r=r er? (4 = ctg 6) (Fig. 32) under the effect of a single central 


@) 
force (see problem 4.2), this problem can only be the force F v 1/r3 
(see also problems 3.9 and 11.6). However, it can be shown that 
this motion is possible under the influence of the attractive 

force F v 1/r2 (force of Newtonian attraction) if to it we add 

the tangential thrust F ~ 1/r2. Let us write out the equations 


of motion in polar coordinates: m(#-r@*)=F,+,, and m(ré +27°@)-R, 


whence the radial and transversal craft accelerations are as follows: 


w.-#-rg?=-+, +4 cose, we = 7G +2rg-K sing. (11.9) 


By differentiating the equation of the spiral, we find /133 
r=7T@ cty p and F = re 2 ety” p+ rg ctg p From the second equation of 
(11.9) we have Beats a(ré+2rg ctgp). By substituting this ratio into 
the first equation of (11.9) we get: 


re (1 +ctg?a) =4; gst, sintp | 


whence 


2G =— 3B sin’p, G=-% @'ctgp =- 3-5, stnp cosp. | 


If from the first equation of (11.9) we determine R/m and 
replace the derivatives F, @° and @ by their expressions, we 


find the acceleration imparted to the craft by the tangential 
thrust, therefore, we find the thrust R: 


m = aap (F-1@* + $s) = £y (1- sin’p + costp - $ cosip) = Lapse / 
R=pmcosp/2r’. 


Thus, for the motion to occur along a logarithmic spiral, 
the tangential thrust R must be a force of repulsion proportional 


TO 1/r. These conditions are satisfied, for example, by the 
pressure of solar radiation so that the theoretical possibility 

of a‘ craft with a solar sail moving in outer space is evident 
when the spiral angle 8 is equal to the angle of sail positioning, 
that is, equal to the angle of incidence of the radiation at the 
sail surface y. Taking the motion to be heliocentric (yp = Uods 


let us determine the craft velocity along a spiral, by replacing 
the derivative r in the formulas by the radial and transversal 


velocities: ; 
v= Fr =r ctg p =|/-Fe- COSBy =r -/Be sing, 
so that 


-Vor vos - fF. (11.10) 
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This means that at each point of the spiral the craft velocity, 
always oriented along the tangent to this spiral, is equal to the 
local circular heliocentric velocity at the given point and 
decreases with growing separation of the craft from the Sun. 

Using this formula we can get the energy characteristics of the 
thrust R by writing the change in the craft's kinetic energy as 

it passes from one point of the spiral to the other as equal to 
the sum of work done by the solar gravity and the work done by the 
thrust: 


ye ab ant es \ 
or 
Ee -a)= me -) Age | 


so that the work done by the thrust 4: =-zmp.(L-4)iis positive 
as the craft goes farther and farther from the Sun (Yr, > rs) (in 


contrast to the work done by the gravity) and is negative with 
increasing proximity to the Sun (Pr, < ry). In the case when the 


craft goes away from the Sun, the kinetic energy drops off. Just 
like the velocity, the work done by the thrust A_ does not depend 
R 


on the shape of the spiral (on the spiral angle 6), that is, does 
not depend on the number of orbits that the craft must make around 
the Sun to pass from one point to another. 


We can easily see that the tangential thrust R is not a 
potential force. We can be convinced of this by writing out the 
projections of the resultant F' = F + R: 


r r3 


: _He moos’p _ Me™ (,_ 1, ' 
fr.-Hem, be = fer (1 Zc py 


Fy = et cosp sing. ! 


Actually, the force Ft is not a potential force since it is 
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impossible to select the force function U' satisfying simultaneously 


. . . f! 2 ay’ Frat. ov 
the conditions of potentiality ‘r * @ | and er 84! . So the 


force R is not a potential force as well. We note that the resul- 


tant F' is also not a central force, as the result of which neither 


the integral of areas nor the Binet formulas are satisfied in this 
problem. 


To determine the time of flight along the spiral it is suffi- 
cient to integrate any of the relations for the velocity projec- 


tions, for example, se = Ve sing. . Introducing the value of r 
from the equation of a spiral into this relation, we get the 
3ArAG/2 o . 
equation em d@ = PFE sinp at . ; from whence we have 
P a(r, er?) _ ar? /135 


“Sryvp, snp 3, cosp 


Thus, to pass from ry to D,, takes the time 


Veo cusp” (11,11) 


which depends on the spiral angle pcocpse §).. If the angle 
B is smail, the spiral is weakly generated, and the time of motion 


t = 2 (7,4 — 7,98) 
3VHe 


corresponds to the case of radial motion gp = 0 (see problem 11.7). 
if the angle @ is large, then the spiral turns sharply and the time 
of motion is long. When £8 = 1/2, we again get the initial circular 
orbit for which it is not possible to determine the time of motion 
based on Eq. (11.11). In this case we must use the formulas in 
Chapter Six. 


is short, where the minimum possible time 


We should note that the explicit or implicit (in terms of r) 
dependence of the motion characteristics on the spiral angle £6 is 
Getermined by their dependence on the parameters of the sail pro- 
ducing the thrust, since in this case we can assume that each 8B 
corresponds to fixed sail parameters. For example, we assume that 
some given "sail factor" of the craft corresponds to 8 = 88° 
(cos 8 = 0.03) (strict correspondence is possible only for the 
case when the spiral angie 8 is equal to the angle y of incidence 
of the rays at the plane surface). Based on Eq. (11.11) we can 
calculate the time of motion of the craft from Earth orbit to Mars 
orbit: 


- 108)22_ 650° &\W2 - 
, --2[(2,28-10°)" = (1,50 -40°)""] 0.98°108 


3 y1327-+108 + 0,03 ~ sec = 1130 days = 3.1 years. 


Here the radius-vector of the craft (heliocentric) will turn by 
the angle 


L45 


Or ve 


fn 3tVite | . : 


8 8 
-_* 30,9840 °71327 107 -0,03! 2 yo _ yz? 
~ 30,03 tm 2(4,50-10%)*/ °. - 


2 
< = Setgp 


8 


that is, in the time of motion from. the Earth orbit to Mars orbit /136 
the craft will make about half a revolution along the spiral 
around the Sun. 


Problem 11.9. Determine the distance r from Earth's center 
to points situated along the same line as the centers of the Earth 
and Moon at which the force of attraction by the Earth of a mass m 
is equal to the force of attraction of this same mass by the Moon. 
Take the Earth-Moon distance to be d = 384,400 km, and the mass 


ratio M,/M,=,1:81.5. (Problem of sphere of attraction.) 


Solution. The equality of forces of attraction of mass m by 
the Earth and by the Moon 


Fel =tF.l = fms fimMe M, Me 


r’ (d-r)jee ore * =r | 
enables us to write out the quadratic equation that has the 
following roots: 
(M, -M,)r?-2M; ra +M,a@*=0, 
d /M_M. 
"2 Wom, (Ms + MM, ) = 


— | 
ae 

JIN 
mM 


Fig. 32 
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Thus, there are at least two equi-attracted points situated along 
the lines of centers of attraction (collinear points) for which — 
(Fig. 33) we have 


a 
fe *TaynM, * 80d. 
.- — _4_«. 


so that the point C is between Earth and Moon at a distance ro oF 
= 346,000 km from the Earth's center and d - roe 38,400 km from 


the Moon's center, while point D lies "beyond the Moon" at the 


distance ry = 430,500 km from Earth's center and Py - d = 46,100 


km from Moon's center. Point D is 7/00 km farther from the Moon's 
center than is the point C. 


A mass m placed at point C experiences equal attraction by 
both Earth and Moon, where these forces are directed in different /137 
directions so that the mass, devoid of velocity, stays at point C. 
A mass placed at point D will also experience equal attraction, 
but the forces in this case are directed in the same direction so 
that the mass, devoid of velocity, will fall toward the Moon under 
the influence of the resultant of these attractive forces. The 
force of attraction at point D is smaller than at point C: 


- 
El tt (1 MITT 56. 
| - Me/ Ms ° 


The conditions of pairwise equality of the attractive forces 
of Earth and Moon enable us to write 


Be od -y, “4s BD ty -a gf Mg 
AC — rT, VMs * AD ~ ro OCU My 


from whence it follows that points C and D divide up the distance 
between Earth's center and Moon's center AB in the ratio YPM; | 


internally and externally. Using the techniques of elementary 
geometry, we can prove that if a sphere is constructed on the 
segment CD as a diameter, the sphere will be the geometrical locus 
of points for each of which the force of Barth attraction is equal 
to the force of Moon attraction. The value of this force differs 
from the value at the neighboring force owing to the change in 
distance (for the same mass m). We call this sphere the sphere 

of attraction of the Moon relative to the Earth. Within this 
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sphere the force of attraction of mass m by the Moon is always 
greater than the force of attraction by the Earth, The radius 
of the sphere of attraction op and the distance of its center E 
from the Moon's center BE = A can be calculated using the formulas 


p=3eD=3(%-%) =! 42,300 km, A=7+p-a@=4740 km, 


from which it follows that the center of the sphere of attraction 
lies "beyond the Moon" by a distance of 3100 km from its surface 
(Re =| 1740 km). 

The concept of the "sphere of attraction of a lesser mass 
relative to a greater" applies to any point masses. To determine 
the parameters of the sphere if we know only the masses and the 


distances between the masses, we can use the following formulas 
that are valid for the Earth-Moon case: 


r,=a-(p-b), %=d+(p+d), | 


d /138 
= +D+Ac= a I —_ 
ao pth= V+¥M./M, ° a+? A= MMs” | 
| 
V li.le 
A = _aMe/Ms_ _o,79 kin = AVM /Ms ( ) 
= { ~Me/Ms “2? P 1-M./Mz ° 


By writing out these formulas for the Sun-Earth case, we can 
determine the parameters of the sphere of attraction of the Earth: 
relative to the Sun: M, = 6°10°! as, My = 271993 Gs, d = 150°10° 
e) 
km, p = 259,500 km, and A = 450 km. Thus, the center of the Earth's 
sphere of attraction relative to the Sun lies within the globe. It 
is interesting to note that the Moon is outside the Earth's sphere 
of attraction relative to the Sun, as a consequence of which the 
Sun attracts the Moon more strongly than does the Earth. Actually, 
let us set up the ratio of the magnitude of the attractive forces 
acting on the Moon: 


lFol =} F st = (M6/Mg) (Tee/Poc) *) 
By adopting the mean distance from Moon to Sun as 150°10° km, we get 


_ = 2 -4955 5.84407 \2 | 
VEL AL = —esioe— (savior) = 2.18 5 


so that the Sun attracts the Moon roughly twice as strongly as does 
the Earth. The stable geocentric motion of the Moon is due to the 
corresponding geocentric velocities and accelerations. 


Using Eqs. (11.12) we can also show that any point K of a 
sphere constructed on diameter CD exhibits the above-indicated 
property of equality of attractive forces. To see this, it is 
sufficient to prove the following equality (see Fig. 33): 


KB = Be Me - | 
‘A const AC =V 4M; =a. ; 


By denoting a?=oM/Ms,| we can rewrite Eqs, (11.12) as 


aa” dat ad+A=d+ gd a | 


P= yaa" faa? *T-aF 


On the other hand, based on the theorem of cosines for the triangles 


KBE and AKE, we have 


(KB)* =A3+ p’-2 pA cosa = (Sx) (a*+ a" 20° 008 a), 


(KA)* = (d+A)° + p- 2(d+A) pcosat = (Hn) (1+ u“*- 2acosa), | 


from whence there follows (KB/KA)* = 0* 


is what we set out to prove. 


» that is, KB/KA = o, which 


Problem 11.10. Determine at which velocity a missile must 
be launched from the Earth's surface aimed at the Moon for it to 
reach the point of equal attraction closest to the Earth and to 
remain in equilibrium at it. 


Solution. As indicated in problem 11.9, the collinear point 
of equal attractions C situated along the line of the centers of 
the Earth and Moon (see Fig. 33) is the point of equal attractions 
closest to the Earth. The distance from the Earth's center to this 
point r, = 346,000 km, and from the Moon's center -- d - ri = 


= 38,400 km. In this problem was also shown that if some mass (a 
missile or a spacecraft) reaches point C with zero velocity rela- 
tive to Earth and Moon, this mass will remain at this point in 
equilibrium by virtue of the above presented equality of attrac- 
tive forces. 


When a missile moves from the Earth's surface to point C, a 
change occurs in its kinetic energy equal to the sum of work done 
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by the attractive forces of Earth and Moon: 


where 


The equality of the missile at point C means that it arrives there 
at zero velocity, so that the launch velocity is 


», «arn, CE = Lela, aa) | 


Rs a-r% ad-Rs 
| 


Thus, knowing that fM, = 398,600 km/sec“, Ry = 6370 kn, /140 


O O 
d = 384,400 kn, ri = 346,000 km, and d -r,. = 38,400 km, let us 


find Vo = 11.04 km/sec. 


C 


Note that the formula g = fm, /R is no longer valid for the 
Oo Oo 
motion of a point in the field of attraction of two centers. 
Assuming the resultant force of attraction of Earth and Moon at 
the surface of the Earth to be equal to the weight of a body, we 
get 


=o M M = 
~mgr" = fru |-g8- + cae | rs 


whence 


aa M 
g* flat ~ ay] 


When M, =0, this formula is converted into the formula eRe = {[M 


which is valid for one attractive center, 


Similarly, we determine the falling velocity of mass m onto 
the Moon if this mass is at the collinear point of equivalent 
attractions D without initial velocity. In this case 


2 2 
TU — IRV y = A- "+ A.- 
2 ee) Fa Fe? 


= { f { ! 
Ao FIM s (aa a) 20 Age mM -gte)>9, 


i 


so that 


s 
e _— 


When 2%! = 1740 km, r, = 430,500 km, and d = 384,400 km, we 


D 
get v = 2.38 km/sec. Obviously, this then must be the launch 
velocity from the surface of the Moon required for the rocket to 
reach the collinear point of equal attractions D, after which it 
falls to the surface of the Moon. 


Problem 1lil.il. <A spacecraft moves along an elliptical trajec- 
tory around an attractive center M with gravitational parameter 
Uy The craft passes near some attractive mass m. Assuming that 


the perturbing mass is in the orbital plane of the craft, determine 
the instantaneous change in the constant of energy dh/dt and the 
semi-major axis da/dt of the craft orbit. 


Solution. The orbital energy of a craft moving in the field /141 
of attraction of the central mass M.is determined by the integral 


a a BBay BE ing 
of energy hom vt SR rn aa We will call the variation in 


h and a for constant r the instantaneous change in these quantities: 


adh dv aa Da*v dv 


geal Saal —  —— 


ade ~ °° at* dt ~~ Pa ° ae? 


where dv/dt is the instantaneous velocity change. 


In this problem (see Fig. 34) this velocity change is defined 
thusly. Suppose the force |F| = u,/o? is a force of attraction 


of the perturbing mass m (m << M) acting on a unit-mass craft at 
the distance p from it. Under the influence of this force 
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tangential and normal craft 
accelerations are induced 


Of these accelerations, only the 
tangential acceleration Ww, can 


alter the orbital energy, so that 


Fig. 34 


2 Hm Y. cos@, : 


p* 


— 
—_ 


cos§, ——-= 


t pe ps 


ah 
at 


From the resulting formulas: it is clear that the orbital 
energy h increases or decreases, depending on the direction of 
craft motion, that is, on the change in the angle 6. The change 
in energy dh/dt depends only on the gravitational parameter of 
the perturbing mass and the distance p at which the craft is 


separated from it. The value of cos 6 must be assumed negative 
if the perturbing force causes orbital rotation in a direction 
opposite to the direction of the orbital motion of the craft. 

In this case the amount of energy (h < 0) and the semi-major axis 
are reduced. It is precisely this case that is shown in Fig. 34. 


Problem 11.12. When an artificial Earth satellite is moving 
in the atmosphere, scattering (dissipation) of the total energy 
of the AES is observed, obeying the law dh/dt = - 2%, where ¢@ is 
the positive scattering function. Determine the nature of the 


evolution of the elliptical orbit of the AES (h=Bh<0).! 


Write out the equation of the energy balance and establish the 
nature of the change in the orbital kinetic and potential energies 
of the AES. 


Solution. Based on the law of scattering, the total energy 
h acted on by atmospheric drag decreases. From the formula h = 
= - um/2a it becomes obvious that the decrease in the negative h 
(elliptical orbit) is accompanied by a decrease in the total semi- 
major axis of the orbit a, that is, the orbit sinks into the denser 
atmospheric layers. 


We can easily show that during this descent of the orbit, 
the apogee sinks faster than the perigee. Actually, let us look 
at one revolution of the orbit: apogee - perigee - apogee. At 


~ 
}H 
a one 
NO 


| 


the perigee the atmospheric density for a given revolution is the 
maximum, therefore the atmospheric drag appearing in the drag 
formula will also be at a maximum. In addition, the drag is a 
function of AES velocity, and its velocity at the perigee is at a 
maximum so that it is obvious that the effect of atmospheric drag 
will be at a maximum precisely at the perigee. As a result, the 
AES partially loses velocity at the perigee, resulting in the 
former apogee altitude no longer being attained. A similar effect 
can be observed for any point on the orbit, but its intensity 
decreases as the satellite moves from the perigee to the apogee 
(the density decreases, since the altitude of the AES increases, 
while the velocity diminishes according to the law of areas). At 
the apogee the deceleration effect will be the smallest, while a 
slight loss in velocity will not substantially influence the 
position of the next perigee. Thus, the apogee sinks faster than 
the perigee, resulting in the elliptical orbit being converted 
into a circular orbit with simultaneous reduction in a and r: 

e > 0, and a > en i 


In spite of the loss of orbital velocity in each revolution, 
overall as the orbit sinks from revolution to revolution the mean /143 
orbital velocity will increase. Actually, let us determine the 
limit of the mean AES velocity (see problem 5.19): 


dim Vay =e dim (¥, + U,) =ylim Vis (F -+) + 
+/ BC ~ a) =P Paes 


that is, the mean velocity tends to the value of the circular 
velocity, which as we know increases as the orbit sinks. Thus, 
the so-called "satellite paradox" can take place, namely, 

that during motion in the atmosphere the satellite experiences: 
acceleration in the direction of its motion. The approximate 
nature of the velocity change for motion through the atmosphere 
is shown in Fig. 35, where n is the number of revolution. 


The increase in the mean orbital velocity entails an increase 
in the kinetic energy T = mv-/2. By writing out further the energy 


balance equation h=T-+-U-= T+t+V= mv°/2 - um/r and substituting 
this equation. into the expression of the law of energy scattering, 
we get 
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Hence it follows that if there 

1s @ Simultaneous increase in 

T due to an increase in the 
velocity v and a decrease in the 
total energy h resulting from the 
scattering of energy, the poten- 
tial energy of the AES V = - U 
must fall off faster than the 
kinetic energy T increases. 


It should be noted that the 
actual slowing down of AES in 
the atmosphere can be viewed as 
the result of applying a con- 
tinuous series of infinitely 
small impulses, where the drag 
is directed along the tangent to this trajectory (tangential 
impulse) and is proportional to the instantaneous velocity. _After 
applying each impulse, we get a new orbit with new constant h. 

The new ellipse is smaller than the former, and the tangent to it 
at the point of impulse application even has one focus coinciding 
with the focus of the initial ellipse. 


The principles presented above can be illustrated by the 
following formulas. Suppose that at the perigee the’ velocity 


P . ? ° 
receives some increment 42% ..so that ¥, =U,+ AY. The sign of 
the increment is not known in advance. From the integral of 
energy when r= const, we have 

~72 _ 2h 
ho =v -S--4£, Ah =2v, Av, = a Aa, 
yAt+WTSr. 
a r_§ a, Aa = Ale 


whence 


AT,=2Aa= 2(2u,4x%,)a/p - 40, av, 0°/p.. 
_ (11.13) 


This means that the increments Ar, and 4r have the same signs, /145 


since the reduction in the velocity at the perigee necessarily 
entails a decrease in ree And this means the descent of the orbit. 


Let us establish a relationship between the velocity incre- 
ments at the perigee and at the apogee. Let us write out the 
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integral of energy for the apogee and the increment Ah for the 
variable re: 


2 
hoev- tk --£ 


a Pra 


‘Ah = 20, AU, + 2h ar, = +; haat, Sie, 


on 


so that 


{ 2 
Ponce (sae - Fy)an = gllng (rida) ane 


~ sar (ry + 2rnr) Ar. 


Substituting in this expression in place of Ar > its value 
from Eq. (11.13), we find 


r 


Av,= — = AY, [ (y+ 7a], | (11.14) 


From this formula it follows that the increments AV and Av 


have. different signs. In addition, the orbit tends to a circular 


orbit, that is, sed and el Therefore the impulse 
relation AU, =-5A Ue obtains. 


Problem 11.13. Determine the local circular and local para- 
bolic velocities of a point with mass m over the surface of the 
Earth at the distance Lo from its center if the motion occurs in 


the equatorial plane of the Earth. The Earth is taken as a com- 
pressed ellipsoid of revolution (spheroid). 


Solution. The potential of a compressed Earth taken as an 
ellipsoid of revolution, to second-order terms relative to the /146 
geometrical compression of the Barth, is of the form 


ye MEaly+t5 i(2e ‘V(- -3sin'g)+ dy 0 (SE V(1- “10 sinty +32 Fan, 
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where J 2 1624'107°, and D = 6°10 
containing the geometric compression a in the first and second 
powers; ¢ is the geometrical latitude of the points; and Ry, is the 
equatorial radius of the Earth. fe) 


are dimensionless coefficients 


Let us write out the formula of the potential for $9 = 0 
(equatorial plane): 


Um TE lt +S a(R) + ae 0 (48)']. | (11.15) 


We know that the acceleration w of a point moving in a field 
of arbitrary potential force F in a system of spherical coordinates 
consists of three components: 


w,. “w,7r® = fe Res W = We ¥ a ¢° ’ ®. * wa 
F, yo ss BU —_ tL OO entite eo ==. 8 
aabh’, F ro} Fe ~ 7 oe ? Fy reose on "| 


Suppose the point executes uniform circular motion in the equatorial 
plane of the Earth taken as an ellipsoid of revolution. The force 
of gravity F® (¢ = 0) in this case, as is true of the force of 
gravity of a spherical Earth, is a central force (see Fig. 36 a, 

b, and c), that is, FY = FO = 0 and FO = FY. Then the radial 
acceleration is simultaneously a normal acceleration, and the 
"Longitudinal" component of acceleration corresponds to the tan- 
gential acceleration which is always zero for uniform circular 


motion. Thus, 


2 | o 0 . ! 
w.e-w =-a tl. ($=) =const <0, w,=0, w,= w,= Sting, ; 
J 


{ 
mm mH 


whence we can obtain the circular velocity 


, = 
v= - 2 (= id oe const >0. 


° ‘ o,po _ ay 
The gravitational force Fe2F = = 


from Eq. (11.15): /147 


corresponding to it is found 
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so that the circular velocity is defined by the formula 


(11.17) 


To calculate the local parabolic velocity of a point, let us 
use the theorem of the change in the kinetic energy of a point as 
it moves “at infinity": 


2 
mus ““par. 4 


2 Z) ron” 


where A. ~ 68 the work done by the gravitational force po in 
0? 
moving the point in the equatorial plane from a given surface of 


the potential level "at infinity", equal to ,y /148 


- U*(o0) - “U'r)<0, | ——— 
; o., - Ts 
In this case v.. = 0, Vieoy= 0, so that Voar =y2 UP (Te) « Using 


Eq. (11.15) (where the distance r = Py)s let us determine the 


L537 
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parabolic velocity: 


(11.18) 
A comparison of the formulas obtained for Vag and Voar shows 
that the familiar formula Yoar *¥2 Vg5 used by us in solving the 


problems in this collection corresponds only to the first term of 
&a series in the expansion of the potential of the compressed Earth, 
that is, it is valid only for the spherical model of the Earth. °. 
Using Eqs. (11.17) and (11.18), let us write out the corrections 


to the values vr = 7.90 km/sec and Viz = 11.19 km/see for points 
on the equator (Po = Ri): 
O 


v, Be (147+ 5D) = 7,90 71 +1624 10% + $40" - 7,90 ¥'1,001627, 
& 

v = 2he (i+ 49+ 5D) =11,19171+542-10 “+ 0,3-10 = 11,197 1,000542. 
$ 


From these corrections it follows that the influence of the non- 
sphericity of the Earth for Vrq is smaller than for Vy: As a 


point ascends above the Earth's surface this influence diminishes, 
which follows from the velocity formulas. 


In passing from the equatorial plane (49 # 0), the force of 
attraction of a compressed ellipsoid of revolution, owing to the 
disruption of symmetry, ceases to be a central force (Fig. 36 da) 

Here two components of the attractive [gravitational] force appear -- 
-- radial 


| 
\ 


2 
Fio= oe = ~ 7H | { +7 (32) (1-3sin“g ) + 0(S2) (1-10 sin’ +3 sint a 


and "latitudinal" 


F< +H ~ “ES - ad ($2) sin i cose + 4p(28)(mo sinecosy - 60sing cosy). 


The “longitudinal" component B for the model of an ellipsoid 


of revolution is always equal to zero since F does not depend on }. 
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Corresponding to these force components is the variable radial 


E; { eu . 
acceleration W.. =sf =a (+= ro and the nonzero acceleration 
F, - _] gu. ‘ . 
Wy = Fe mrs (SE) ron eee? so that circular motion proves 


to be impossible for an ellipsoid of revolution when ¢$ # 0, 


In this problem we have limited ourselves to examining the 
model of an ellipsoid of revolution, or a spheroid, taking the 
corresponding expansion to an accuracy of second-order terms. But 
if in the expansion (11.15) we take series terms that enable us not 
only to allow for the higher orders of compression for the model of 
a spheroid but also terms characterizing the triaxiality of the 
Earth (nonsphericity of the equator), the oblateness of the 
northern hemisphere, and so on, it becomes obvious that since 
in actuality the Earth's symmetry relative to the equatorial plane 
does not exist, neither does the strictly circular motion of an 
AES exist even in the equatorial plane. So for actual AES motions 
we can speak only of near-circular orbits. 


Problem 11.14. A point is moving under the effect of the 


a-mf{Ht. 
central attractive force * m (Es + a). where yp > 0. and v > O are 


certain constants. Set up an equation of motion of the point using 
Binet's formula. 


Solution. The physical significance of the problem posed 
can be established by using the formulas from problem 11.13 where 
Eq. (11.16) was derived for the force of attraction in the equa- 
torial plane. From Eq. (11.16) it follows that this force, cal- 
culated to the first degree of the geometrical compression of the 
Earth, is of the form 


ea mht [1+ 9(4)], 70. 


This law of action of the force coincides with the law given in 
the conditions of the problem, if we take 


pespy>O, | and =p, IR,>0., 


{ 


Thus, the first summand can be interpreted es the force of attrac- /150 
tion of a spherical Earth, and the second summand can be inter- 

preted as the incremental force of attraction of the equatorial 

access of the Earth's mass. 
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Using Binet's second formula (see Chapter Four), let us write 


dtu ; FF ¢h  Y\V td oP Jur? 
fa ths male teow a 


Qe 


where u = 1/r is Binet's variable, and ¢=7"¢! is the constant of 


| 


areas... Hence 


du -=-‘Y pe | 
age teu tun or | (11.19) 


(compare with Eqs. (11.4) from problem 11.6, and (4.11) from prob- 
lem 4.10). This is a second-order nonlinear inhomogeneous equation 


of the type yl oye dytcy+d=0 (a +0), , which is insolvable in elementary 

functions. It can, for example, be reduced to an equation of the 
oe , 

type y'* + gaytZby*rcy +edy+c=0 (a =0), that is solvable in 


elliptical functions. However, from the physical point of view 
this equation is the equation of the perturbed motion of a point 


‘ain a central field of attraction. In celestial mechanics, solutions 


of equations of perturbed motion are set up by using methods of 
successive approximations or methods of the variation of arbitrary 
constants. If these methods are applied to the solution of this 
equation (11.19), one must take as the first approximation the 
solution of the equation of unperturbed motion derived from (11.19) 
when v = Q. This equation corresponds to the motion of a point 

in the field of attraction of a spherical Earth and its solution 

is an unperturbed Keplerian ellipse. 


Problem 11.15. Determine the maximum possible period of 
revolution of an artificial Earth satellite T. What is the maximum 
duration of residence of this AES in the Earth's umbra? 


Solution. We wiil assume the surface of the Earth's sphere 
of action to be the boundary of the domain of existence of the AES. 
Taking as the radius of the sphere p = 929,900 km (see Chapter Nine), 
which is 2.41 times greater than the distance of Earth to Moon 


re =}384,400 km, let us express the period of revolution of this 
AES in terms of the known period of revolution of the Moon around 
the Earth /« =| 27.3 days (see problem 10.7), using Kepler's third 
law: : 


T=T, (£)*- (2,41) °"T, = 3,72°27,.3 = 104.6 days. 


Obviously, when r > o, the body will be converted into an f/151 
artificial planet so that the period found can be assumed To be 
approximately of the possible periods. 


Let-.us determine the duration of residence of an AES in the 
Rarthts umbra. The length of the conical shadow of the Earth 
OC ee (Fig, 37) can be found from the congruence of the triangles 


O'K'C' and OKC, By setting O'C = K'C = a@+x,|, we can write 


Ro [Rem +Xyx- Taking R, = 6,96°10? km, R, = 6370 km, and 
° O 
a = 149.6*10° km, we get Bs tO 12496 +10" + : 
6 ; °. 
ins _ 2 : 
from whence xe te = 46 38> 0" km, 


2. Earth 
3. Orbit of AES 


The half-aperture angle of the cone of the Earth's umbra can 
be determined from AOKC: 


-. & _ Rg _— 6,37°10> _ a _ ont 
sin ar = gy = -agsige” 79.0046, = =0 16 | 


We can readily see that aa the length of the cone of the 
Earth's umbra is ywe 1.38°10°% km, an object at a distance p = 


930,000 km = 0.93°106 km from the Earth's center can fall within 
this cone, for a specific orbital inclination. The time spent in 
the umbra for an AES orbit situated in the plane of the ecliptic 
will be at a maximum when the orbital plane will lie in the plane 
of symmetry of the umbra cone (this orbital inclination is equal 
to the inclination of the ecliptic to the equator i = 23°.5). 
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The time of residence of this AES in the umbra can be deter- 
mined as follows. Let us find the length of the segment AB 
(chords of circular segment) from the congruence of triangles OBD 


and OKC, It is 


AB L 2-P , ~~ 
2Re ~ =~ cos $ ‘\ {152 
whence when cos St | we have 


We determine the corresponding central angle: 
: AB _ 4, 
sin £ = 6° Stgsa0r a 0 = 90,0022 , 
£ 20°08’, 6 =0°s6! = 0,0047. 


The proportion of the above-found period of revolution of 
the AES corresponding to its residence in the Earth's umbra is 


T | 0047 _ : 
+ =f. = = 2 = 0,00075, |! 
whence the time of residence in the umbra (maximum possible) is 


T= 0,00075 7 = 0,00075 -401,6 = 0.076 day = 1 hour 50 min, 
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